CN115903032A - Underground fluid storage space and migration channel detection method and device - Google Patents
Underground fluid storage space and migration channel detection method and device Download PDFInfo
- Publication number
- CN115903032A CN115903032A CN202211307019.1A CN202211307019A CN115903032A CN 115903032 A CN115903032 A CN 115903032A CN 202211307019 A CN202211307019 A CN 202211307019A CN 115903032 A CN115903032 A CN 115903032A
- Authority
- CN
- China
- Prior art keywords
- matrix
- data
- storage space
- fluid storage
- migration
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 32
- 238000013508 migration Methods 0.000 title claims abstract description 32
- 239000012530 fluid Substances 0.000 title claims abstract description 24
- 238000001514 detection method Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 10
- 238000003384 imaging method Methods 0.000 claims abstract description 10
- 238000005516 engineering process Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 36
- 230000006870 function Effects 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 6
- 238000012216 screening Methods 0.000 claims description 4
- 230000002159 abnormal effect Effects 0.000 abstract description 5
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 4
- 239000011707 mineral Substances 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000007789 sealing Methods 0.000 abstract description 2
- 238000010438 heat treatment Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the technical field of detection, and relates to a method and a device for detecting an underground fluid storage space and a migration channel, wherein common offset seismic wave field data of a region to be processed are collected firstly; processing the common offset seismic wave field data to obtain scattered wave data; then, performing homing imaging on scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, and realizing the detection of the underground fluid storage space and the migration channel; the method takes the kinematic characteristic difference of the wave field and the dynamic characteristic of the wave field into consideration, so that the wave field information can be extracted more quickly and conveniently, the matrix decomposition technology is applied, the scattered wave data of discontinuous geologic bodies can be accurately acquired, the calculation efficiency and the accuracy are high, the method is beneficial to positioning underground abnormal geologic bodies, and the method is convenient for geothermal heating, mineral resources and CO 2 And items such as sealing and the like are safely and effectively promoted.
Description
The technical field is as follows:
the invention belongs to the technical field of detection, and relates to a method and a device for detecting an underground fluid storage space and a migration channel.
Background art:
underground discontinuous geological anomalous bodies have important influence on the development and utilization of geothermy, mineral resources and the like. The discontinuous geological abnormal body can be used as a geothermal storage space heat conduction channel for heat conduction, and is of great importance for geothermal resource utilization and well drilling position selection. The underground discontinuous fractures serving as storage spaces and migration channels of mineral resources such as coal bed gas, petroleum and natural gas play a vital role in selecting drilling positions and exploring, developing and utilizing oil and gas resources. The scattering data is the wave field response of the underground geological abnormal body, can be used for accurately positioning the non-uniform discontinuous body and provides stronger underground space illumination, and the premise of using the scattering data is to extract the scattering data from the wave field data and then carry out high-precision imaging to accurately position and image the geological abnormal body.
In the existing method, a plane wave decomposition method is usually adopted to obtain scattering data, but the method only considers the kinematic characteristic difference of a wave field, ignores the dynamic characteristic of the wave field and is not beneficial to the extraction of wave field information.
The invention content is as follows:
the invention aims to overcome the defects in the prior art, and provides a method and a device for detecting an underground fluid storage space and a migration channel.
In order to achieve the purpose, the specific process of detecting the underground fluid storage space and the migration channel comprises the following steps:
(1) Acquiring common offset distance seismic wave field data of a region to be processed;
(2) Processing the common offset seismic wave field data to obtain scattered wave data;
(3) And performing homing imaging on the scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, so as to realize detection of the underground fluid storage space and the migration channel.
As a further technical scheme of the invention, the specific process of acquiring the scattered wave data in the step (2) is as follows:
(21) Constructing frequency space domain seismic data by utilizing Fourier transform on the common offset seismic wave data acquired in the step (1), and converting the frequency space domain seismic data into a Hankel matrix;
(22) Selecting a plurality of rows and columns from a Hankel matrix according to a matrix row-column contribution value weight function to construct a decomposition submatrix;
(23) And approximating the low-rank information of the Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component, and acquiring scattered wave data.
As a further technical scheme of the invention, the Hankle matrix obtained in the step (21) is as follows:
wherein d is i For the ith channel of data at a certain frequency, m and n are the rows and columns of the matrix.
As a further technical solution of the present invention, the matrix row and column contribution weight functions in step (22) are respectively:
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row and column contribution weight function.
As a further technical solution of the present invention, the approximate low rank component shown in step (23)Scattered wave data->
As a further technical scheme of the invention, the shift algorithm adopted in the step (3) is a Kirchoff shift algorithm.
Compared with the prior art, the method simultaneously considers the kinematic characteristic difference of the wave field and the dynamic characteristic of the wave field, so that the wave field information is extracted more quickly and conveniently, the matrix decomposition technology is used, the scattered wave data of discontinuous geologic bodies can be accurately obtained, the calculation efficiency and the accuracy are high, the method is favorable for positioning underground abnormal geologic bodies, and is convenient for geothermal heat, mineral resources and CO 2 And items such as sealing and the like are safely and effectively promoted.
Description of the drawings:
FIG. 1 is a block diagram of the working process of the present invention.
Fig. 2 is a block diagram of a process for acquiring scattered wave data according to the present invention.
FIG. 3 is a block diagram of the underground fluid storage space and migration channel detection apparatus according to the present invention.
FIG. 4 is a raw seismic stack data record of common-offset seismic wavefield data of example 1 of the present invention.
FIG. 5 is a graph of scattered wave data extracted in embodiment 1 of the present invention.
FIG. 6 is a diagram showing the result of the scattered wave data imaging in example 1 of the present invention.
The specific implementation mode is as follows:
to make the objects, technical solutions and advantages of the embodiments of the present invention clearer and more complete, the technical solutions of the present invention will be described below with reference to the accompanying drawings, it is obvious that the described embodiments are some, but not all embodiments of the present invention, and all other embodiments obtained by a person skilled in the art without making creative efforts based on the embodiments of the present invention belong to the protection scope of the present invention.
Example 1:
as shown in fig. 1 and fig. 2, the present embodiment provides a method for detecting an underground fluid storage space and a migration channel, which includes the following specific processes:
(1) Determining a region to be processed, and acquiring common offset seismic wave field data of the region to be processed, as shown in figure 1;
(2) Constructing frequency-space domain seismic data by utilizing Fourier transform according to the common offset seismic wave field data acquired in the step (1), and converting the frequency-space domain seismic data into a Hankel matrix H:
wherein d is i The ith channel of data under a certain frequency, and m and n are rows and columns of a matrix;
(3) Selecting a plurality of rows and columns from a Hankel matrix to construct a decomposition submatrix according to a matrix row and column contribution value weight function, wherein the matrix row and column contribution weight function is respectively as follows:
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method comprises the following steps of obtaining three matrixes A, U and B through matrix row and column contribution weight functions as a generalized inverse matrix;
(4) According to the submatrix obtained in the step (3), utilizing a matrix decomposition technology to approximate Hankle matrix low-rank information to obtain approximate low-rank componentsThereby acquiring scattered wave data->As shown in fig. 5;
(5) And (3) performing homing imaging on the scattered wave data obtained in the step (4) by using a Kirchhoff migration algorithm, and obtaining a fine positioning result of the underground fluid storage space and the migration channel as shown in fig. 6 to realize detection of the underground fluid storage space and the migration channel.
Example 2:
as shown in fig. 3, the present embodiment provides an underground fluid storage space and migration channel detection apparatus, including:
the data acquisition module is used for acquiring common offset seismic wave field data of the area to be detected;
the data processing module is used for processing the common offset seismic wave field data to obtain scattered wave data;
and the imaging module is used for carrying out offset imaging on the scattered wave data to obtain a fine detection result of the underground fluid storage space and the migration channel.
Specifically, the data processing module comprises a Hankel matrix conversion unit, a sub-matrix decomposition unit and a scattered wave data acquisition unit, the Hankel matrix conversion unit constructs frequency space domain seismic data from the collected common offset seismic wave field data by utilizing Fourier transform and converts the frequency space domain seismic data into a Hankel matrix H, and the sub-matrix decomposition unit selects a plurality of rows and columns from the Hankel matrix to construct a decomposition sub-matrix and construct a target function by utilizing a matrix row-column contribution value weight function; the scattered wave data acquisition unit approximates the low-rank information of a Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component, so that the scattered wave data is acquired.
More specifically, the Hankel matrix H obtained by the Hankel matrix conversion unit is:
wherein d is i The ith channel of data under a certain frequency, and m and n are rows and columns of a matrix;
matrix row and column contribution weight functions adopted by the sub-matrix decomposition unit are respectively as follows:
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and the matrix U is formed by intersecting elements of the A and the B, namely:
U=A + HB +
wherein, (.) + The matrix is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row-column contribution weight function;
approximate low rank component obtained by scattered wave acquisition unitThereby obtaining scattered wave data
In the embodiment, the common offset seismic wave field data of the area to be processed is obtained through a data acquisition module; the data processing module constructs frequency-space domain seismic data by utilizing Fourier transform and converts the frequency-space domain seismic data into a Hankel matrix; selecting a plurality of rows and columns from the Hankel matrix according to the matrix row and column contribution value weight function to construct a decomposition submatrix; then, approximating Hankle matrix low-rank information by using a matrix decomposition technology to obtain scattering wave field data; and finally, the imaging module performs homing imaging on the scattered wave data through a migration algorithm to obtain detection results of the underground fluid storage space and the migration channel.
The computer program product of the method and the device for finely detecting the underground fluid storage space and the migration channel provided by the embodiment of the invention comprises a computer readable storage medium storing program codes, wherein instructions included in the program codes can be used for executing the method described in the previous embodiment, specific implementation can refer to the method embodiment, and detailed description is omitted, and algorithms not described in detail are all general technologies in the field.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk, and various media capable of storing program codes.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: those skilled in the art can still make modifications or changes to the embodiments described in the foregoing embodiments, or make equivalent substitutions for some features, within the scope of the disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. A method for detecting underground fluid storage space and migration channels is characterized by comprising the following specific processes:
(1) Acquiring common offset distance seismic wave field data of a region to be processed;
(2) Processing the common offset seismic wave field data to obtain scattered wave data;
(3) And performing homing imaging on the scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, so as to realize detection of the underground fluid storage space and the migration channel.
2. The method for detecting the underground fluid storage space and the migration channel according to claim 1, wherein the step (2) of acquiring the scattered wave data comprises the following specific steps:
(21) Constructing frequency space domain seismic data by utilizing Fourier transform on the common offset seismic wave data acquired in the step (1), and converting the frequency space domain seismic data into a Hankel matrix;
(22) Selecting a plurality of rows and columns from the Hankel matrix according to a matrix row and column contribution value weight function to construct a decomposition submatrix;
(23) And approximating the low-rank information of the Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component and acquiring scattered wave data.
4. A method according to claim 3, wherein the matrix row and column contribution weighting functions of step (22) are:
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row and column contribution weight function.
6. The method for detecting underground fluid storage space and migration passage according to claim 1, wherein the migration algorithm used in step (3) is Kirchhoff migration algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211307019.1A CN115903032B (en) | 2022-10-25 | 2022-10-25 | Underground fluid storage space and migration channel detection method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211307019.1A CN115903032B (en) | 2022-10-25 | 2022-10-25 | Underground fluid storage space and migration channel detection method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115903032A true CN115903032A (en) | 2023-04-04 |
CN115903032B CN115903032B (en) | 2023-12-15 |
Family
ID=86470065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211307019.1A Active CN115903032B (en) | 2022-10-25 | 2022-10-25 | Underground fluid storage space and migration channel detection method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115903032B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018129844A1 (en) * | 2017-01-10 | 2018-07-19 | 中国科学院地质与地球物理研究所 | Seismic diffracted wave separation method and device |
CN111999770A (en) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | Precise beam offset imaging method and system for converting TTI medium into PS wave |
CN112684504A (en) * | 2020-12-03 | 2021-04-20 | 中国矿业大学 | Urban underground cavity rapid detection method based on total scattering model |
CN113945967A (en) * | 2021-10-14 | 2022-01-18 | 中国矿业大学(北京) | Diffraction wave separation method and device |
CN113945969A (en) * | 2021-10-19 | 2022-01-18 | 中国矿业大学(北京) | Discontinuous geologic body detection method and device and electronic equipment |
CN114924315A (en) * | 2022-05-10 | 2022-08-19 | 北京世纪金道石油技术开发有限公司 | Polymorphic seismic exploration method and system |
-
2022
- 2022-10-25 CN CN202211307019.1A patent/CN115903032B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018129844A1 (en) * | 2017-01-10 | 2018-07-19 | 中国科学院地质与地球物理研究所 | Seismic diffracted wave separation method and device |
CN111999770A (en) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | Precise beam offset imaging method and system for converting TTI medium into PS wave |
CN112684504A (en) * | 2020-12-03 | 2021-04-20 | 中国矿业大学 | Urban underground cavity rapid detection method based on total scattering model |
CN113945967A (en) * | 2021-10-14 | 2022-01-18 | 中国矿业大学(北京) | Diffraction wave separation method and device |
CN113945969A (en) * | 2021-10-19 | 2022-01-18 | 中国矿业大学(北京) | Discontinuous geologic body detection method and device and electronic equipment |
CN114924315A (en) * | 2022-05-10 | 2022-08-19 | 北京世纪金道石油技术开发有限公司 | Polymorphic seismic exploration method and system |
Non-Patent Citations (3)
Title |
---|
史恒;刘洪;高星;丁仁伟;王伟;: "散射波叠前时间偏移成像", 石油地球物理勘探, no. 04 * |
李纯阳,周官群,章俊: "散射波成像技术在地铁空洞探测中的应用", 安徽理工大学学报( 自然科学版), vol. 38, no. 3, pages 40 - 44 * |
王瑞贞;白旭明;王金宽;赵利庆;韩力;邱文平;: "朝克乌拉凹陷火成岩覆盖区地震勘探方法研究", 石油物探, no. 03 * |
Also Published As
Publication number | Publication date |
---|---|
CN115903032B (en) | 2023-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Simmons et al. | LLNL‐G3Dv3: Global P wave tomography model for improved regional and teleseismic travel time prediction | |
Afonso et al. | 3‐D multi‐observable probabilistic inversion for the compositional and thermal structure of the lithosphere and upper mantle. II: General methodology and resolution analysis | |
Xiao et al. | Siamese earthquake transformer: A pair‐input deep‐learning model for earthquake detection and phase picking on a seismic array | |
Yang et al. | Dual-mode FPGA implementation of target and anomaly detection algorithms for real-time hyperspectral imaging | |
Hu et al. | Using deep learning to derive shear‐wave velocity models from surface‐wave dispersion data | |
Zhang et al. | Extracting dispersion curves from ambient noise correlations using deep learning | |
US10436925B2 (en) | Method for processing marine seismic data | |
Zhang et al. | Moho depth variations from receiver function imaging in the northeastern North China Craton and its tectonic implications | |
Yuan et al. | A probabilistic shear wave velocity model of the crust in the central West Australian craton constrained by transdimensional inversion of ambient noise dispersion | |
Kumar et al. | Dynamic time warping as an alternative to windowed cross correlation in seismological applications | |
CN112114362A (en) | Method for reconstructing space-time field of underground shallow layer explosion | |
Sopher | Converting scanned images of seismic reflection data into SEG-Y format | |
Chevrot et al. | Optimized discrete wavelet transforms in the cubed sphere with the lifting scheme—implications for global finite-frequency tomography | |
Li | Recent advances in earthquake monitoring II: Emergence of next-generation intelligent systems | |
Wang et al. | Denoising with weak signal preservation by group-sparsity transform learning | |
Guevara et al. | A data-driven workflow for predicting horizontal well production using vertical well logs | |
López et al. | Spectral gap-based seismic survey design | |
Ojeda G et al. | Spatio‐temporal entropy analysis of the magnetic field to help magnetic cloud characterization | |
AlSalmi et al. | Automated seismic semantic segmentation using attention U-Net | |
CN115903032A (en) | Underground fluid storage space and migration channel detection method and device | |
Gu et al. | A modified excitation amplitude imaging condition for prestack reverse time migration | |
Liu et al. | Line survey joint denoising via low-rank minimization | |
Song et al. | Domain adaptation in automatic picking of phase velocity dispersions based on deep learning | |
CN114185095B (en) | Method for suppressing multiple waves of three-dimensional plane wave domain seismic data | |
CN113484906B (en) | High-resolution energy field reconstruction method based on low-frequency energy spectrum data driving |
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 |