CN115718327B - Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium - Google Patents
Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN115718327B CN115718327B CN202310036668.0A CN202310036668A CN115718327B CN 115718327 B CN115718327 B CN 115718327B CN 202310036668 A CN202310036668 A CN 202310036668A CN 115718327 B CN115718327 B CN 115718327B
- Authority
- CN
- China
- Prior art keywords
- rock
- integrity
- coefficient
- function
- rock mass
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000010606 normalization Methods 0.000 claims abstract description 42
- 238000012545 processing Methods 0.000 claims abstract description 19
- 230000001131 transforming effect Effects 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 48
- 230000008569 process Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000009412 basement excavation Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 238000002759 z-score normalization Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 5
- 238000009826 distribution Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 229910052742 iron Inorganic materials 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010224 classification analysis Methods 0.000 description 1
- 238000007621 cluster analysis 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
- 238000013519 translation Methods 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 mass integrity coefficient acquisition method, a rock mass integrity coefficient acquisition system, electronic equipment and a storage medium. The method utilizes all resistivity data in the acquired tunnel section to carry out primary standardization processing according to a standardization function; and transforming the standardized function in a form of firstly translating and then scaling so as to ensure that the maximum value of the transformed function is 10 times of the minimum value, performing re-normalization processing, converting the re-normalization function into an initial value of the rock integrity coefficient, and determining an adjustment coefficient of the integrity profile to obtain a final rock integrity coefficient. When the method is used for forming the resistivity calculation integrity coefficient calculation formula, the overall resistivity data can be considered, and the difficulty of data interpretation is simplified.
Description
Technical Field
The invention relates to the technical field of tunnel exploration, in particular to a rock mass integrity coefficient acquisition method, a system, electronic equipment and a storage medium.
Background
In the process of iron tunnel investigation, in order to find out lithology limit, structure and other bad geological conditions in the range of tunnel body and evaluate the integrity of the rock mass of the tunnel body, a geophysical prospecting method is often adopted for detection; the current method for tunnel geophysical prospecting mainly comprises an earthquake method and an electromagnetic method, wherein the earthquake method mainly obtains wave velocity parameters of a rock mass, and the electromagnetic method mainly obtains resistivity parameters of the rock mass. The adoption of wave velocity parameters for rock mass grading or rock mass integrity evaluation has achieved a number of productive research results and has been incorporated into the relevant survey design specifications. At present, rock mass grading or rock mass integrity evaluation by adopting resistivity parameters is in a research stage, and a plurality of problems still exist in a practical stage. The rock mass integrity of a deeply buried tunnel is studied by adopting the electromagnetic resistivity of the earth in the middle-iron four-hospital Seey and the like (invention patent CN 102495430B), the true resistivity of the rock mass is obtained by utilizing the weighted average of multiple points of 20 meters near the buried depth elevation of the tunnel, and the rock resistivity and the weak medium resistivity in the rock mass are obtained by an outcrop small quadrupole method, so that the rock mass integrity coefficient is obtained. China iron design Xu Anchun and the like research a method for judging the grade of a tunnel rock mass based on magnetotelluric sounding resistivity (application number CN 110968840A), after a plurality of normal distributions are adopted, the resistivity intervals are relatively concentrated, 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 rock integrity by resistivity still have a plurality of limitations in use: firstly, the calculated rock 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 global property of rock evaluation is lacking; secondly, the adopted calculation formula has more parameters, the difficulty of acquiring the parameters is high, and the extreme value easily appears in the calculation process; and thirdly, the magnitude difference of the resistivity profile is larger, which is not beneficial to classification and cluster analysis.
Therefore, the application provides a rock mass integrity coefficient acquisition method, a rock mass integrity coefficient acquisition system, electronic equipment and a storage medium, and the problems are solved.
Disclosure of Invention
Therefore, the invention aims to provide a rock mass integrity coefficient acquisition method, a system, electronic equipment and a storage medium, which are used for carrying out numerical conversion after standardized processing, so that the converted numerical values are distributed in the range of [0,1], and when a resistivity calculation integrity coefficient calculation formula is formed, global resistivity data can be considered, and the difficulty of data processing is simplified.
In order to achieve the above object, the method for obtaining the integrity coefficient of the rock mass of the present invention comprises the following steps:
s1, acquiring all resistivity data rho in a tunnel section, and performing primary standardization processing according to a standardization function, wherein the primary standardization function is represented by f (rho);
s2, transforming the standardized function in a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, performing re-normalization on the resistivity data subjected to the first normalization, and recording the re-normalization function as f g (ρ)；
S3, normalizing the function f again g Conversion of (ρ) to initial value K 'of rock integrity factor' r According to the determined adjustment coefficient, the rock mass is completedAnd (5) adjusting the initial value of the integrity coefficient to obtain the final rock integrity coefficient.
Further preferably, in S1, the first normalization process is to perform normalization by using a Z-score normalization method after removing extreme values and abnormal values from all resistivity data.
Further preferably, the reconverted normalization function f g (ρ) is related to the first normalization function f (ρ) according to the following formula:
f g (ρ)＝(f(ρ)+a)*b；
wherein ,
f(ρ) max representing the maximum value after the first normalization treatment; f (ρ) min Representing the minimum value after the first normalization process.
Further preferably, in S3, the function f g Conversion of (ρ) to initial value K 'of rock integrity factor' r The following formula is adopted:
wherein ,fg (ρ) is a re-normalization function, f g (ρ) max The maximum value after the re-normalization treatment is shown.
Further preferably, in S3, the value of the adjustment coefficient λ is 0.75, and the adjustment coefficient is obtained by performing fitting correction on the rock integrity coefficient by using existing resistivity profile data and rock grading data revealed by excavation.
Further preferably, the final rock mass integrity factor is represented by K using the formula r ＝λ*K′ r Wherein Kr is the final rock mass integrity factor, K' r For the initial value of the rock mass integrity coefficient, λ=0.75 represents the adjustment coefficient.
Further preferably, the method also comprises S4, utilizing the rock mass integrity coefficient K r And (5) grading rock mass:
when the rock mass integrity coefficient is >0.75, the rock mass grade belongs to grade i;
when the rock integrity coefficient is more than 0.60 and less than or equal to 0.75, the rock grade belongs to grade II;
when the rock integrity coefficient is more than 0.40 and less than or equal to 0.60, the rock grade belongs to grade III;
when the rock integrity coefficient is more than 0.20 and less than or equal to 0.40, the rock grade belongs to IV grade;
when the rock integrity coefficient is less than or equal to 0.20, the rock grade belongs to the V grade.
The invention also provides a rock mass integrity coefficient calculation system, which is used for implementing the rock mass integrity coefficient acquisition 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 the standardization function;
the data processing module is used for transforming the standardized function in a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, re-normalizing the resistivity data after the first normalization, and recording the normalized function after the transformation as f g (ρ); will function f g (ρ) converting into an initial value of a rock mass integrity coefficient;
the data fitting module is used for carrying out fitting correction on the rock integrity coefficient according to the existing resistivity profile data and 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.
The present invention also provides an electronic device including: a memory and a processor; the memory is used for storing at least one group of instruction sets; the processor is used for calling and executing the instruction set in the memory, and executing the rock mass integrity coefficient acquisition method through executing the instruction set.
The invention also provides a storage medium storing a computer program of the rock mass integrity coefficient acquisition method, the computer program being executed by a processor to implement the rock mass integrity coefficient acquisition method.
Compared with the prior art, the rock mass integrity coefficient acquisition method disclosed by the application has the following advantages:
when the resistivity is obtained, all resistivity data in the tunnel section are adopted, the global property of the rock mass is evaluated, powerful data support is carried out, and due to the fact that the data discrete degree is large and extreme value characteristics exist, on the other hand, when the normalization processing is carried out, the two normalization processing processes are adopted, the range of the data is adjusted through the first normalization processing, the transformed data relationship can be suitable for the logarithmic relationship during the second normalization processing, so that the global resistivity data can be converted by utilizing the accordant logarithmic relationship, the global resistivity data can be considered, the difficulty of data interpretation is simplified, and the clustering analysis is facilitated and the calculation process is reduced.
Drawings
FIG. 1 is a flow chart of a rock mass integrity factor acquisition method of the present invention;
FIG. 2 is a schematic representation of all resistivity data within a tunnel profile acquired in accordance with the present application;
FIG. 3 is a graph of the first normalized resistivity value distribution characteristics of the present application;
FIG. 4 is a re-normalization function f in the present application g (ρ) correspondence with the integrity coefficient;
FIG. 5 is K in the present application r Corresponding relation with resistivity;
FIG. 6 is a schematic cross-sectional view of a tunnel integrity factor in accordance with the present invention;
fig. 7 is a cross-sectional view of a certain tunnel rock mass grade in the present invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and the detailed description.
As shown in fig. 1, a method for obtaining a rock integrity coefficient according to an embodiment of the present invention includes the following steps:
s1, acquiring all resistivity data rho in a tunnel section, and performing primary standardization processing according to a standardization function, wherein the primary standardization function is represented by f (rho);
s2, transforming the standardized function in a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, performing re-normalization on the resistivity data subjected to the first normalization, and recording the re-normalization function as f g (ρ)；
S3, normalizing the function f again g Conversion of (ρ) to initial value K 'of rock integrity factor' r 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.
As shown in fig. 2 and 3, since the degree of data dispersion is large and there is a characteristic of an extreme value, in S1, the first normalization process is to perform normalization processing by using a Z-score normalization method after removing the extreme value and the abnormal value from all the resistivity data. Assuming that the resistivity is a variable ρ, the normalization is performed using the Z-score method,
in the formula
ρ: profile resistivity;
mu: average value of resistivity
Sigma: standard deviation of resistivity.
This way of processing will give the data a characteristic that the average value of the data must be 0 and the standard deviation must be 1. The value of f (ρ) function greater than 3 is assigned 3 and the value less than-3 is assigned-3.
Converting f (ρ) function, translating the whole, and scaling to convert f g The maximum value of (p) is 10 times the minimum value.
The maximum value and the minimum value are 10 times, so that f is convenient for subsequent calculation g (ρ) min ＝1，f g (ρ) max =10. Assuming a translation amount a, a scaling amount b is:
f g (ρ) min ＝(f(ρ) min +a)*b＝1
f g (ρ) max ＝(f(ρ) max +a)*b＝10
the deduction is as follows:
f g (ρ)＝(f(ρ)+a)*b
in the formula ：
as shown in fig. 4, the distribution rule of the function value and the integrity coefficient after conversion approximately accords with the logarithmic relationship.
f(ρ) max Representing the maximum value after the first normalization treatment; f (ρ) min Representing the minimum value after the first normalization process.
Will function f g Conversion of (ρ) to initial value K 'of rock integrity factor' r The following formula is adopted:
wherein ,fg (ρ) is a function value when the resistivity is ρ after the reconversion, f g (ρ) max The maximum value after the re-normalization treatment is shown.
In S3, the value of the adjustment coefficient λ is 0.75, and the adjustment coefficient is obtained by performing fitting correction on the rock integrity coefficient by using existing resistivity profile data and rock grading data revealed by excavation. As shown in FIG. 5, according to K r And correcting the corresponding relation with the resistivity. In order to accommodate the differences between different regions, the adjustment coefficient may be adjusted in a small range, where λ=0.75±β, where β is an adaptation parameter adjusted according to lithology components of different regions.
The final rock integrity factor is represented by K using the formula r ＝λ*K′ r Wherein Kr is the final rock mass integrity factor, K' r For the first time of rock integrity coefficientThe value, λ=0.75, represents the adjustment coefficient.
As shown in fig. 6 and 7, the method also comprises S4, and utilizes the rock mass integrity coefficient K r And (5) grading rock mass:
when the rock mass integrity coefficient is >0.75, the rock mass grade belongs to grade i;
when the rock integrity coefficient is more than 0.60 and less than or equal to 0.75, the rock grade belongs to grade II;
when the rock integrity coefficient is more than 0.40 and less than or equal to 0.60, the rock grade belongs to grade III;
when the rock integrity coefficient is more than 0.20 and less than or equal to 0.40, the rock grade belongs to IV grade;
when the rock integrity coefficient is less than or equal to 0.20, the rock grade belongs to the V grade.
The invention also provides a rock mass integrity coefficient calculation system, which is used for implementing the rock mass integrity coefficient acquisition 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 the standardization function;
the data processing module is used for transforming the standardized function in a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, re-normalizing the resistivity data after the first normalization, and recording the normalized function after the transformation as f g (ρ); will function f g (ρ) converting into an initial value of a rock mass integrity coefficient;
the data fitting module is used for carrying out fitting correction on the rock integrity coefficient according to the existing resistivity profile data and 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.
The present invention also provides an electronic device including: a memory and a processor; the memory is used for storing at least one group of instruction sets; the processor is used for calling and executing the instruction set in the memory, and executing the rock mass integrity coefficient acquisition method through executing the instruction set.
The invention also provides a storage medium storing a computer program of the rock mass integrity coefficient acquisition method, the computer program being executed by a processor to implement the rock mass integrity coefficient acquisition method.
It is apparent that the above examples are given by way of illustration only and are not limiting of the embodiments. Other variations or modifications of the above teachings will be apparent to those of ordinary skill in the art. It is not necessary here nor is it exhaustive of all embodiments. While still being apparent from variations or modifications that may be made by those skilled in the art are within the scope of the invention.
Claims (9)
1. The rock mass integrity coefficient acquisition method is characterized by comprising the following steps of:
s1, acquiring all resistivity data rho in a tunnel section, and performing primary standardization processing according to a standardization function, wherein the primary standardization function is represented by f (rho);
s2, transforming the first normalization function f (ρ) according to a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, and performing re-normalization on the resistivity data subjected to the first normalization to obtain a re-normalization function denoted as f g (ρ)；
The re-normalization function f g (ρ) is related to the first normalization function f (ρ) according to the following formula:
f g (ρ)＝(f(ρ)+a)*b；
wherein ,
f(ρ) max representing the maximum value after the first normalization treatment; f (ρ) min Representing the minimum value after the first normalization treatment;
s3, normalizing the function f again g Conversion of (ρ) to initial value K 'of rock integrity factor' r Based on the determined adjustment coefficient, the rock integrity systemAnd (5) adjusting the initial value of the number to obtain a final rock mass integrity coefficient.
2. The method according to claim 1, wherein in S1, the first normalization process is a normalization process performed by using a Z-score normalization method after extremum and outlier are removed from all resistivity data.
3. The rock mass integrity factor obtaining method according to claim 1, wherein in S3, a function f is applied g Conversion of (ρ) to initial value K 'of rock integrity factor' r The following formula is adopted:
wherein ,fg (ρ) is a re-normalization function, f g (ρ) max The maximum value after the re-normalization treatment is shown.
4. The method for obtaining the rock mass integrity coefficient according to claim 1, wherein in S3, the value of the adjustment coefficient λ is 0.75, and the adjustment coefficient is obtained by fitting and correcting the rock mass integrity coefficient by using existing resistivity profile data and rock mass grading data revealed by excavation.
5. The method of claim 4, wherein the final rock integrity factor is represented by K using the formula r ＝λ*K′ r wherein ,Kr For the final rock mass integrity factor, K' r For the initial value of the rock mass integrity coefficient, λ=0.75 represents the adjustment coefficient.
6. The method of claim 5, further comprising S4 using a rock integrity factor K r Grading rock massDividing:
when the rock mass integrity coefficient is >0.75, the rock mass grade belongs to grade i;
when the rock integrity coefficient is more than 0.60 and less than or equal to 0.75, the rock grade belongs to grade II;
when the rock integrity coefficient is more than 0.40 and less than or equal to 0.60, the rock grade belongs to grade III;
when the rock integrity coefficient is more than 0.20 and less than or equal to 0.40, the rock grade belongs to IV grade;
when the rock integrity coefficient is less than or equal to 0.20, the rock grade belongs to the V grade.
7. A rock mass integrity factor acquisition system for carrying out the rock mass integrity factor acquisition method of any of the preceding claims 1-6, 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, wherein the primary standardization function is represented by f (rho);
the data processing module is used for transforming the standardized function in a form of firstly translating and then zooming so that the maximum value of the transformed function is 10 times of the minimum value, re-normalizing the resistivity data after the primary normalization, and recording the re-normalized function as f g (ρ); will function f g (ρ) converting into an initial value of a rock mass integrity coefficient;
the data fitting module is used for carrying out fitting correction on the rock integrity coefficient according to the existing resistivity profile data and 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.
8. An electronic device, comprising: a memory and a processor; the memory is used for storing at least one group of instruction sets; the processor is configured to invoke and execute the instruction set in the memory, and execute the rock mass integrity coefficient acquisition method according to any one of claims 1-6 by executing the instruction set.
9. A storage medium storing a computer program of the rock mass integrity coefficient acquisition method of any one of claims 1-6, the computer program being executed by a processor to implement the rock mass integrity coefficient acquisition method of any one of claims 1-6.
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 CN115718327A (en) | 2023-02-28 |
| CN115718327B true 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) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AP2016009458A0 (en) * | 2014-02-28 | 2016-09-30 | Barnard Andries Jacobus | Grouted rock support testing apparatus and method |
| CN112819240B (en) * | 2021-02-19 | 2021-11-26 | 北京科技大学 | Method for predicting shale oil yield based on physical constraint LSTM model |
| CN113010942B (en) * | 2021-02-25 | 2022-10-25 | 中国铁路设计集团有限公司 | Tunnel excavation risk early warning and surrounding rock grading evaluation method |
-
2023
- 2023-01-10 CN CN202310036668.0A patent/CN115718327B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN115718327A (en) | 2023-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Marwan et al. | Cross recurrence plot based synchronization of time series | |
| CN110865357A (en) | Laser radar echo signal noise reduction method based on parameter optimization VMD | |
| CN111597981B (en) | Magnetotelluric signal denoising method and system based on improved multi-scale diffusion entropy | |
| CN108230365B (en) | SAR image change detection method based on multi-source difference image content fusion | |
| CN117368920B (en) | D-insar-based coal mining area subsidence monitoring method and system | |
| CN114595732A (en) | Radar radiation source sorting method based on depth clustering | |
| CN118010849B (en) | Expressway bridge and culvert damage detection method and system | |
| CN114429151A (en) | Magnetotelluric signal identification and reconstruction method and system based on depth residual error network | |
| CN116068619A (en) | Self-adaptive multi-order frequency dispersion surface wave pressing method, device and equipment | |
| CN115718327B (en) | Rock mass integrity coefficient acquisition method, system, electronic equipment and storage medium | |
| CN106548031A (en) | A kind of Identification of Modal Parameter | |
| CN114091538B (en) | Intelligent noise reduction method for discrimination loss convolutional neural network based on signal characteristics | |
| CN113514743A (en) | Construction method of GIS partial discharge pattern recognition system based on multi-dimensional features | |
| CN114488293A (en) | High-resolution inversion method based on sensitive logging curve | |
| Ding et al. | Automatic scaling of F2-layer parameters from ionograms based on the empirical orthogonal function (EOF) analysis of ionospheric electron density | |
| CN116484184A (en) | Method and device for enhancing partial discharge defect sample of power equipment | |
| CN116343051A (en) | Geological environment monitoring method and system based on remote sensing image | |
| CN109029735B (en) | Earth surface emissivity calculation method and device | |
| CN105093312A (en) | Seismic relative wave impedance prediction method and device based on frequency domain multi-order differentiation | |
| CN115293090A (en) | Multi-harmonic source responsibility quantification method based on reconstruction data processing | |
| CN112630839B (en) | Logging curve standardization method and system | |
| CN110794469B (en) | Gravity inversion method based on minimum geological feature unit constraint | |
| Qiu et al. | Iteration empirical mode decomposition method for filling the missing data of GNSS position time series | |
| CN112684493B (en) | Solid tide response removing method based on improved variation modal decomposition | |
| CN118013301B (en) | BIM-based highway bridge construction information digital management 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 |