CN115718327A - Rock integrity coefficient calculation method and system, electronic equipment and storage medium - Google Patents
Rock integrity coefficient calculation method and system, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN115718327A CN115718327A CN202310036668.0A CN202310036668A CN115718327A CN 115718327 A CN115718327 A CN 115718327A CN 202310036668 A CN202310036668 A CN 202310036668A CN 115718327 A CN115718327 A CN 115718327A
- Authority
- CN
- China
- Prior art keywords
- integrity
- rock mass
- rock
- coefficient
- function
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 133
- 238000004364 calculation method Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000013519 translation Methods 0.000 claims abstract description 7
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 41
- 238000010606 normalization Methods 0.000 claims description 25
- 238000009412 basement excavation Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000002759 z-score normalization Methods 0.000 claims description 3
- 230000002159 abnormal effect Effects 0.000 claims description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The invention discloses a rock integrity coefficient calculation method, a rock integrity coefficient calculation system, electronic equipment and a storage medium. According to the method, all resistivity data in a tunnel section are acquired, and primary standardization processing is carried out according to a standardization function; and transforming the standardized function according to a mode of translation and scaling, then carrying out standardization processing again after the maximum value of the transformed function is 10 times of the minimum value, converting the standardized function again into an initial value of the integrity coefficient of the rock mass, and determining an adjustment coefficient of the integrity profile to obtain the final integrity coefficient of the rock mass. When the resistivity calculation integrity coefficient calculation formula is formed, the method not only can give consideration to the overall resistivity data, but also simplifies the difficulty of data interpretation.
Description
Technical Field
The invention relates to the technical field of tunnel exploration, in particular to a rock integrity coefficient calculation method, a rock integrity coefficient calculation system, electronic equipment and a storage medium.
Background
In the iron tunnel exploration process, in order to find lithological boundaries, structures and other unfavorable geological conditions in the tunnel body range and evaluate the integrity of the tunnel body rock mass, a geophysical prospecting method is often adopted for detection; at present, methods for tunnel geophysical prospecting mainly comprise an earthquake method and an electromagnetism method, wherein the earthquake method mainly obtains wave velocity parameters of a rock body, and the electromagnetism method mainly obtains resistivity parameters of the rock body. The wave velocity parameter is adopted to carry out rock mass grading or rock mass integrity evaluation, so that a plurality of fruitful research results are obtained, and relevant investigation and design specifications are included. At present, the resistivity parameter is adopted to evaluate the rock mass grading or the rock mass integrity in the research stage, and a plurality of problems still exist in the practical stage. The method comprises the steps of carrying out research on the integrity of a rock body of a deep-buried tunnel by adopting earth-electric magnetic resistivity (patent CN 102495430B), obtaining the true resistivity of the rock body by utilizing a weighted average value of multiple points 20 meters near the buried depth elevation of the tunnel, and obtaining the resistivity of the rock and the resistivity of soft and weak media in the rock body by a small outcrop quadrupole method, so as to obtain the integrity coefficient of the rock body. The method for determining the grade of the tunnel rock mass based on the magnetotelluric sounding resistivity (application number CN 110968840A) is researched by China iron, namely, xiguangchun and the like, after multiple times of normal distribution is adopted, the resistivity intervals are relatively centralized, and then the corresponding relation between the grade of the tunnel rock mass and the resistivity intervals is established. The two representative methods for evaluating the integrity of the rock mass by adopting the resistivity still have a plurality of limitations in use: firstly, the calculated rock mass integrity coefficient only selects the resistivity data near the tunnel body, so that the resistivity information outside the tunnel body is easy to ignore, and the globality of rock mass evaluation is lacked; secondly, the adopted calculation formula has more parameters, the difficulty in obtaining the parameters is high, and extreme values are easy to occur in the calculation process; thirdly, the numerical magnitude difference of the resistivity profile is large, and classification and cluster analysis are not facilitated.
Therefore, the application provides a rock integrity coefficient calculation method, and the problems are solved.
Disclosure of Invention
Therefore, the invention aims to provide a rock integrity coefficient calculation method, system electronic equipment and a storage medium, which are used for carrying out numerical conversion after standardization processing so that the converted numerical values are distributed in the range of [0,1], and can not only give consideration to global resistivity data but also simplify the difficulty of data processing when forming a resistivity calculation integrity coefficient calculation formula.
In order to achieve the purpose, the rock integrity coefficient calculation method comprises the following steps:
s1, acquiring all resistivity data in a tunnel section
Carrying out primary standardization treatment according to a standardization function;
s2, converting the normalized function according to a mode of translation and scaling so that the maximum value of the converted function is 10 times of the minimum value, carrying out re-normalization processing on the resistivity data subjected to the first normalization processing, and recording the re-normalized function as
；
S3, normalizing the function again
Converted into initial value of rock mass integrity coefficient
And adjusting the initial value of the rock integrity coefficient according to the determined adjustment coefficient to obtain the final rock integrity coefficient.
Further preferably, in S1, the first normalization process is performed by performing normalization using a Z-score normalization method after removing extrema and outliers from all resistivity data.
Further preferably, the re-transformed normalization function
Establishing and first normalizing function according to the following formula
The relationship of (1):
wherein ,
，
；
represents the maximum value after the re-normalization process;
the minimum value after the re-normalization process is shown.
Further preferably, in S3, the function is divided
Converted into initial value of rock mass integrity coefficient
The following formula is adopted:
wherein ,
after being transformed again, the resistivity is
The function value of the time is calculated,
after representing the re-normalizationA maximum value;
the minimum value after the re-normalization process is shown.
Further preferably, in S3, the adjustment coefficient
The value of (2) is 0.75, and the adjustment coefficient is obtained by fitting and correcting the integrity coefficient of the rock mass by adopting the existing resistivity profile data and the rock mass grading data revealed by excavation.
Further preferably, the final rock mass integrity factor is expressed by the following formula
wherein ,
in order to obtain the final rock mass integrity factor,
is the initial value of the integrity coefficient of the rock mass,
indicating the adjustment factor.
Further preferably, the method also comprises S4 and utilizes the integrity coefficient of the rock mass
Grading rock mass:
when the integrity coefficient of the rock mass is greater than 0.75, the rock mass grade belongs to the I grade;
when the integrity coefficient of the rock mass is more than 0.60 and less than or equal to 0.75, the rock mass grade belongs to II grade;
when the integrity coefficient of the rock mass is more than 0.40 and less than or equal to 0.60, the rock mass grade belongs to grade III;
when the integrity coefficient of the rock mass is more than 0.20 and less than or equal to 0.40, the rock mass grade belongs to IV grade;
and when the integrity coefficient of the rock mass is less than or equal to 0.20, the rock mass grade belongs to the V grade.
The invention also provides a rock mass integrity coefficient calculation system for implementing the rock mass integrity coefficient calculation method, which comprises the following steps:
the data acquisition module is used for acquiring all resistivity data in the tunnel section and carrying out primary standardization processing according to a standardization function;
the data processing module is used for transforming the standardized function according to a mode of firstly translating and then scaling so that the maximum value of the transformed function is 10 times of the minimum value, carrying out standardization processing on the resistivity data subjected to the first standardization processing again, and recording the transformed standardized function as
(ii) a A function of
Converting the initial value of the integrity coefficient of the rock mass;
and the data fitting module is used for fitting and correcting the rock integrity coefficient according to the existing resistivity profile data and the rock mass grading data revealed by excavation, determining the adjustment coefficient of the integrity profile, and adjusting the initial value of the initial rock integrity coefficient to obtain the final rock integrity coefficient.
The present invention also provides an electronic device, comprising: a memory and a processor; the memory for storing at least one set of instructions; the processor is used for calling and executing the instruction set in the memory, and executing the rock integrity coefficient calculation method by executing the instruction set.
The invention also provides a storage medium, which stores a computer program of the rock integrity coefficient calculation method, wherein the computer program is executed by a processor to realize the rock integrity coefficient calculation method.
Compared with the prior art, the rock mass integrity coefficient calculation method has the following advantages:
when the resistivity is obtained, all resistivity data in a tunnel section are adopted, the global property of a rock mass is evaluated, powerful data support is carried out, the data discrete degree is large, the extreme value characteristic exists, on the other hand, two times of standardization processing processes are adopted in standardization processing, the range of the data is adjusted through the first time of standardization processing, the transformed data relation can be suitable for the logarithmic relation in the second time of standardization processing, the global resistivity data can be converted by utilizing the logarithmic relation which accords with, the global resistivity data can be considered, the difficulty of data interpretation is simplified, the clustering analysis is facilitated, and the calculation process is reduced.
Drawings
FIG. 1 is a flow chart of a method for calculating the integrity coefficient of a rock mass according to the invention;
FIG. 2 is a graphical representation of overall resistivity data within a tunnel section acquired by the present application;
FIG. 3 is a graph of the first normalized resistivity value distribution characteristic of the present application;
FIG. 4 is a re-normalization function of the present application
Corresponding to the integrity coefficient;
FIG. 5 is a drawing of the present application
Corresponding to the resistivity;
FIG. 6 is a schematic cross-sectional view of a tunnel integrity factor according to the present invention;
FIG. 7 is a cross-sectional view of a rock mass grade of a tunnel according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and the detailed description.
As shown in fig. 1, an embodiment of the invention provides a method for calculating a rock integrity coefficient, which includes the following steps:
s1, acquiring all resistivity data in a tunnel section
Carrying out primary standardization treatment according to a standardization function;
s2, converting the normalized function according to a mode of translation and scaling so that the maximum value of the converted function is 10 times of the minimum value, carrying out re-normalization processing on the resistivity data subjected to the first normalization processing, and recording the re-normalized function as
；
S3, normalizing the function again
Converted into initial value of rock mass integrity coefficient
And adjusting the initial value of the rock integrity coefficient according to the determined adjustment coefficient to obtain the final rock integrity coefficient.
As shown in fig. 2 and 3, since the data dispersion degree is large and the extreme value characteristic exists, in S1, the first normalization process is to perform normalization process by using a Z-score normalization method after removing extreme values and abnormal values from all resistivity data. Assuming resistivity as a variable
The Z-score method is adopted for standardization,
in the formula
This process gives the data a characteristic that the mean value of the data is always 0 and the standard deviation is always 1. Will be provided with
The value of the function greater than 3 is assigned a value of 3 and the value less than-3 is assigned a value of-3.
Will be provided with
The function is converted, and the whole translation is firstly carried out and then the scaling is carried out, so that the converted function
The maximum is 10 times the minimum.
The maximum value and the minimum value are 10 times of the relation, so that the subsequent calculation is convenient
，
. Assuming that the translation amount a and the scaling amount b are as follows:
derived to
in the formula ：
as shown in fig. 4, the distribution rule of the converted function value and the integrity coefficient approximately conforms to a logarithmic relationship.
A function of
Converted into initial value of rock mass integrity coefficient
The following formula is adopted:
wherein ,
after being transformed again, the resistivity is
The function value of the time is calculated,
represents the maximum value after the re-normalization process;
the minimum value after the re-normalization process is shown.
In S3, the adjustment coefficient
The adjustment coefficient is obtained by fitting and correcting the integrity coefficient of the rock mass by adopting the existing resistivity profile data and the rock mass grading data revealed by excavation. According to K, as shown in FIG. 5 r And correcting the corresponding relation with the resistivity. It should be noted that if the adjustment coefficient can be adjusted in a small range to adapt to the difference of different regions, in this case
, wherein ,
is an adaptive parameter adjusted according to lithologic compositions of different regions.
The final rock integrity coefficient is expressed by the following formula
wherein ,
in order to obtain the final integrity factor of the rock mass,
is the initial value of the integrity coefficient of the rock mass,
indicating the adjustment factor.
As shown in fig. 6 and 7, the method further comprises S4, utilizing the integrity coefficient of the rock mass
Grading rock mass:
when the integrity coefficient of the rock mass is greater than 0.75, the rock mass grade belongs to grade I;
when the integrity coefficient of the rock mass is more than 0.60 and less than or equal to 0.75, the rock mass grade belongs to II grade;
when the integrity coefficient of the rock mass is more than 0.40 and less than or equal to 0.60, the rock mass grade belongs to grade III;
when the integrity coefficient of the rock mass is more than 0.20 and less than or equal to 0.40, the rock mass grade belongs to IV grade;
and when the integrity coefficient of the rock mass is less than or equal to 0.20, the rock mass grade belongs to the V grade.
The invention also provides a rock integrity coefficient calculation system, which is used for implementing the rock integrity coefficient calculation method and comprises the following steps:
the data acquisition module is used for acquiring all resistivity data in the tunnel section and carrying out primary standardization processing according to a standardization function;
the data processing module is used for transforming the standardized function according to a mode of translation and then scaling so that the maximum value of the transformed function is 10 times of the minimum value, the resistivity data subjected to the first standardization processing is subjected to the second standardization processing, and the transformed standardized function is recorded as
(ii) a Will function
Converting the initial value of the rock mass integrity coefficient;
and the data fitting module is used for fitting and correcting the rock integrity coefficient according to the existing resistivity profile data and the rock mass grading data revealed by excavation, determining the adjustment coefficient of the integrity profile, and adjusting the initial value of the initial rock integrity coefficient to obtain the final rock integrity coefficient.
The present invention also provides an electronic device comprising: a memory and a processor; the memory for storing at least one set of instructions; the processor is used for calling and executing the instruction set in the memory, and executing the rock integrity coefficient calculation method by executing the instruction set.
The invention also provides a storage medium, which stores a computer program of the rock integrity factor calculation method, wherein the computer program is executed by a processor to realize the rock integrity factor calculation method.
It should be understood that the above examples are only for clarity of illustration and are not intended to limit the embodiments. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications therefrom are within the scope of the invention.
Claims (10)
1. A rock integrity coefficient calculation method is characterized by comprising the following steps:
s1, acquiring all resistivity data in a tunnel section
Carrying out primary standardization treatment according to a standardization function;
s2, converting the normalized function according to a mode of translation and scaling so that the maximum value of the converted function is 10 times of the minimum value, carrying out re-normalization processing on the resistivity data subjected to the first normalization processing, and recording the re-normalized function as
；
S3, normalizing the function again
Converted into initial value of rock mass integrity coefficient
And adjusting the initial value of the rock mass integrity coefficient according to the determined adjustment coefficient to obtain the final rock mass integrity coefficient.
2. The method for calculating the rock integrity coefficient according to claim 1, wherein in S1, the first normalization process is performed by using a Z-score normalization method after removing extreme values and abnormal values from all resistivity data.
3. A method of calculating a rock mass integrity factor according to claim 1 wherein the re-normalisation function
Establishing and first normalizing function according to the following formula
The relationship of (c):
wherein ,
，
；
represents the maximum value after the re-normalization process;
represents the minimum value after the re-normalization process.
4. A method of calculating a rock mass integrity factor in accordance with claim 1, wherein in S3 the function is applied
Converted into initial value of rock mass integrity coefficient
The following formula is adopted:
wherein ,
after being transformed again, the resistivity is
The function value of the time is calculated,
represents the maximum value after the re-normalization process;
the minimum value after the re-normalization process is shown.
5. The method of calculating a rock integrity coefficient of claim 1 wherein in S3 the adjustment coefficient
The value of (2) is 0.75, and the adjustment coefficient is obtained by fitting and correcting the integrity coefficient of the rock mass by adopting the existing resistivity profile data and the rock mass grading data revealed by excavation.
6. The method of calculating a rock mass integrity coefficient of claim 5 wherein the final rock mass integrity coefficient is expressed by the following formula
wherein ,
in order to obtain the final rock mass integrity factor,
is the initial value of the integrity coefficient of the rock mass,
indicating the adjustment factor.
7. The method for calculating the integrity coefficient of a rock mass according to claim 1, further comprising S4, using the integrity coefficient of a rock mass
Grading rock mass:
when the integrity coefficient of the rock mass is greater than 0.75, the rock mass grade belongs to grade I;
when the integrity coefficient of the rock mass is more than 0.60 and less than or equal to 0.75, the rock mass grade belongs to II grade;
when the integrity coefficient of the rock mass is more than 0.40 and less than or equal to 0.60, the rock mass grade belongs to grade III;
when the integrity coefficient of the rock mass is more than 0.20 and less than or equal to 0.40, the rock mass grade belongs to IV grade;
and when the integrity coefficient of the rock mass is less than or equal to 0.20, the rock mass grade belongs to the V grade.
8. A rock integrity factor calculation system for carrying out the rock integrity factor calculation method of any one of claims 1 to 7, comprising:
the data acquisition module is used for acquiring all resistivity data in the tunnel section and carrying out primary standardization processing according to a standardization function;
the data processing module is used for transforming the standardized function according to a mode of firstly translating and then scaling so that the maximum value of the transformed function is 10 times of the minimum value, carrying out standardization processing on the resistivity data subjected to the first standardization processing again, and recording the transformed standardized function as
(ii) a A function of
Converting the initial value of the integrity coefficient of the rock mass;
and the data fitting module is used for fitting and correcting the rock integrity coefficient according to the existing resistivity profile data and the rock grading data revealed by excavation, determining the adjustment coefficient of the integrity profile, and adjusting the initial value of the initial rock integrity coefficient to obtain the final rock integrity coefficient.
9. An electronic device, comprising: a memory and a processor; the memory for storing at least one set of instructions; the processor is used for calling and executing the instruction set in the memory, and the rock integrity coefficient calculation method according to any one of claims 1 to 7 is executed by executing the instruction set.
10. A storage medium storing a computer program of the method of calculating a rock integrity factor of any one of claims 1 to 7, the computer program being executable by a processor to implement the method of calculating a rock integrity factor of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310036668.0A CN115718327B (en) | 2023-01-10 | 2023-01-10 | Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310036668.0A CN115718327B (en) | 2023-01-10 | 2023-01-10 | Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115718327A true CN115718327A (en) | 2023-02-28 |
| CN115718327B CN115718327B (en) | 2023-05-02 |
Family
ID=85257978
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310036668.0A Active CN115718327B (en) | 2023-01-10 | 2023-01-10 | Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115718327B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015128831A1 (en) * | 2014-02-28 | 2015-09-03 | BARNARD, Andries Jacobus | Grouted rock support testing apparatus and method |
| CN112819240A (en) * | 2021-02-19 | 2021-05-18 | 北京科技大学 | Method for predicting shale oil yield based on physical constraint LSTM model |
| CN113010942A (en) * | 2021-02-25 | 2021-06-22 | 中国铁路设计集团有限公司 | Tunnel excavation risk early warning and surrounding rock grading evaluation method |
-
2023
- 2023-01-10 CN CN202310036668.0A patent/CN115718327B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015128831A1 (en) * | 2014-02-28 | 2015-09-03 | BARNARD, Andries Jacobus | Grouted rock support testing apparatus and method |
| CN112819240A (en) * | 2021-02-19 | 2021-05-18 | 北京科技大学 | Method for predicting shale oil yield based on physical constraint LSTM model |
| CN113010942A (en) * | 2021-02-25 | 2021-06-22 | 中国铁路设计集团有限公司 | Tunnel excavation risk early warning and surrounding rock grading evaluation method |
Non-Patent Citations (2)
| Title |
|---|
| JIACHUANG WANG: "Research on Rock Mass Quality Classification Based on An Improved Rough Set Cloud Model" * |
| 化希瑞: "基于岩体电性参数的深埋隧道围岩完整性判识方法" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115718327B (en) | 2023-05-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110361778B (en) | Seismic data reconstruction method based on generation countermeasure network | |
| CN107688201B (en) | RBM-based seismic prestack signal clustering method | |
| CN110687610B (en) | Gravity and magnetic data correlation analysis-based field source positioning and attribute identification method | |
| CN114595732B (en) | Radar radiation source sorting method based on depth clustering | |
| CN108230365B (en) | SAR image change detection method based on multi-source difference image content fusion | |
| CN106548031A (en) | A kind of Identification of Modal Parameter | |
| CN110490061B (en) | Uncertainty modeling and measuring method for remote sensing image characteristics | |
| CN116432988A (en) | Intelligent management method, medium and equipment for valve production process data | |
| CN116304549A (en) | Wavelet threshold denoising method for tunnel health monitoring data | |
| CN114091538B (en) | Intelligent noise reduction method for discrimination loss convolutional neural network based on signal characteristics | |
| Li et al. | Magnetotelluric data denoising method combining two deep-learning-based models | |
| CN115718327A (en) | Rock integrity coefficient calculation method and system, electronic equipment and storage medium | |
| CN113514743A (en) | Construction method of GIS partial discharge pattern recognition system based on multi-dimensional features | |
| CN116310851B (en) | Remote sensing image change detection method | |
| CN116776181A (en) | Terminal side load identification method, medium and system based on improved fuzzy clustering | |
| CN116629686A (en) | Method and device for evaluating enterprise energy consumption data | |
| CN116484184A (en) | Method and device for enhancing partial discharge defect sample of power equipment | |
| CN115204277A (en) | Asymptotic optimal lithology prediction method based on multi-threshold quantitative output observation | |
| Ng et al. | Reconstructing ice‐flow fields from streamlined subglacial bedforms: A kriging approach | |
| CN110794469B (en) | Gravity inversion method based on minimum geological feature unit constraint | |
| CN115374687A (en) | Numerical-shape combined intelligent diagnosis method for working conditions of oil well | |
| CN113033908B (en) | Optimization method for submarine resource detection station | |
| CN104614768A (en) | Linear and nonlinear combined seismic wavelet phase correction method | |
| CN113837163B (en) | Tunnel monitoring method and system based on three-dimensional ground penetrating radar and storage medium | |
| CN112183473B (en) | Geological curved surface visual semantic feature extraction 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 |