CN110646848A - Seismic interpretation method for fault - Google Patents
Seismic interpretation method for fault Download PDFInfo
- Publication number
- CN110646848A CN110646848A CN201910977084.7A CN201910977084A CN110646848A CN 110646848 A CN110646848 A CN 110646848A CN 201910977084 A CN201910977084 A CN 201910977084A CN 110646848 A CN110646848 A CN 110646848A
- Authority
- CN
- China
- Prior art keywords
- fault
- coherent
- section
- seismic
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001427 coherent effect Effects 0.000 claims description 38
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000012986 modification Methods 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 238000003708 edge detection Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000011179 visual inspection Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 abstract 8
- 230000005540 biological transmission Effects 0.000 abstract 2
- 238000007670 refining Methods 0.000 abstract 2
- 238000007599 discharging Methods 0.000 abstract 1
- 238000004513 sizing Methods 0.000 abstract 1
- 239000003245 coal Substances 0.000 description 7
- 238000005065 mining Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the technical field of oil exploration, in particular to a fault earthquake interpretation method which comprises a base, a material supporting disc, a machine frame, a roller, a transmission system, a gas transmission pipeline and a control panel, the control panel is provided with an operating button for controlling each part of the thermal refining machine to work, the upper end surface of the base is detachably provided with a material supporting disc, two machine frames are fixedly arranged on the upper end surface of the base, the two machine frames are respectively arranged on two sides of the material supporting disc, two ends of the rollers are respectively and rotatably arranged on the two frames through bearings, the number of the rollers is two, the two rollers are arranged in parallel, through the arrangement of the positive and negative rotation motors, the screw and the charging basket, the auxiliary materials are sprinkled onto the sizing material from the discharging port of the charging basket moving left and right, therefore, auxiliary materials are uniformly added, the rubber material can be accelerated to reach an expected mixing state, and the heat refining efficiency of the rubber material is improved.
Description
Technical Field
The invention relates to the technical field of petroleum exploration, in particular to a fault seismic interpretation method.
Background
China is a country which takes coal as a main energy source, the requirement on coal mine exploration is increasingly improved along with the requirement of mechanized coal mining of mines in recent years, faults are non-negligible disastrous geological abnormal bodies in the coal mining process, the existence of the faults breaks the continuity of coal seams, influences the stability of surrounding rocks of a roadway, hinders mechanized coal mining, can form good water guide channels or gas gathering places, becomes potential risks of coal seam roof and floor water inrush and gas inrush in the mine construction and production processes, different design schemes and pretreatment measures can be adopted according to the size of the faults, so that only the correct identification and fault interpretation are carried out to determine the production state and the property of the faults, and reliable geological basis can be provided for the optimization design, the reasonable arrangement of mining areas and working faces,
with the continuous improvement of seismic exploration technology, the requirement on seismic data interpretation precision is higher and higher, interpretation work of large faults is basically completed, the production capacity and benefit of a fully mechanized mining face depend on the pinpointing degree of small faults to a great extent, so the interpretation work of small faults is more and more important.
Disclosure of Invention
In order to solve the technical problem, the invention provides a fault seismic interpretation method.
The technical problem to be solved by the invention is realized by adopting the following technical scheme:
a method of seismic interpretation of a fault, the method comprising the steps of:
(1) fault section interpretation: the method is characterized in that a color variable density section is adopted to identify small faults, some small faults may have distorted wave crests but broken wave troughs in a section picture, the faults are on a conventional section, because people are used to pay attention to characteristic information of the wave crests and often ignore wave trough information, the situation is difficult to judge whether the situation is caused by fault layers or lithological change, and on the color variable density section, the wave crests and the wave troughs are filled with different colors;
(2) tomographic three-dimensional interpretation: the essence of the coherent body technology is that the abnormal phenomena of outstanding correlation and irrelevance between traces are calculated by utilizing seismic information, and a corresponding and effective calculation method is selected from a three-dimensional data body to realize the conversion of the coherent data body, so that the fracture characteristics are shown and the plane distribution of a crack development zone is predicted;
(3) the darker area on the coherent body slice represents lower coherence degree, the brighter area represents higher coherence degree, because the reflected wave continuity at the fault is poor, the fault can be explained on the coherent body slice according to the dark line or curve-shaped characteristic, at the initial stage of the explanation, time slice explanation can be made at intervals of large time-point, because the large fault extends far in the longitudinal direction and has reaction in the range of 10-100ms, the large fault can be firstly explained and then gradually encrypted until one horizontal time slice is formed every 1ms, so that after the explanation of the small fault is recognized on the slice, the fault can project a breakpoint at the corresponding position on the cross section, and the main measuring line, the contact measuring line and various arbitrary cross sections can be used for checking and correcting the explained breakpoint on the cross section, the combination of the breakpoints is carried out, so that small faults which are not easy to directly observe on the section are solved;
(4) using all the information of the data to automatically track the layer group to form three-dimensional data of the horizon surface, and adopting automatic superposition processing of surface slices or combination of human-computer interaction surface slices to form t according to the multi-solution degree of the horizon0Time, three-dimensional visual inspection, projection to vertical section inspection and modification to form t0Color pixel structure diagram, for t0Time or t0The color pixel map carries out other conventional interpretation work, and the target fault track lines on the seismic main section and the contact section are corrected and interpreted.
Furthermore, in the conversion process of the coherent data volume, an estimated value of three-dimensional seismic correlation can be obtained by calculating the local waveform similarity in the longitudinal direction and the transverse direction, seismic traces in a small range cut by a fault plane usually have different characteristics from adjacent traces, so that the local trace-to-trace correlation is suddenly changed, the correlation value on each grid point is calculated along a time slice, a profile with a low correlation value along the fault can be obtained, the process is repeated for a series of time slices, the profiles become the sections through the extraction of the three-dimensional correlation attribute body, and the three-dimensional reflection amplitude data volume can be converted into the three-dimensional correlation coefficient data volume.
Further, the coherent body algorithm has C1,C2,C3A total of 3 algorithms, C1The correlation between the transverse measuring line and the adjacent longitudinal measuring line of each channel is calculated by a coherent algorithm, two results are mixed and standardized by energy, and the essence is the expansion of the two correlation algorithms; c2The coherent algorithm is to carry out coherent calculation on any multi-channel seismic data based on a similar coherent algorithm; c1The coherent algorithm is implemented by means of a covariance matrix, since C2Coherent algorithms can greatly suppress noise when the effective signal is greater than the average of the noise, so C3The coherent algorithm has higher horizontal resolution and vertical resolution than the other two coherent algorithms in fault identification and edge detection.
Compared with the prior art, the invention has the beneficial effects that: the invention relates to a fault seismic interpretation method, a fault section interpretation method and a fault three-dimensional interpretation method, which utilize a plurality of groups of data to correct and interpret target fault track lines on a seismic main section and a contact section, can use main survey lines, contact survey lines and various arbitrary sections on the sections to check and correct explained breakpoints, and carry out the combination of the breakpoints, thereby interpreting and interpreting small faults which are not easy to be directly observed on the sections.
Detailed Description
The following will clearly and completely describe the technical solutions in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
The invention provides a technical scheme that: a method of seismic interpretation of a fault, the method comprising the steps of:
(1) fault section interpretation: the method is characterized in that a color variable density section is adopted to identify small faults, some small faults may have distorted wave crests but broken wave troughs in a section picture, the faults are on a conventional section, because people are used to pay attention to characteristic information of the wave crests and often ignore wave trough information, the situation is difficult to judge whether the situation is caused by fault layers or lithological change, and on the color variable density section, the wave crests and the wave troughs are filled with different colors;
(2) tomographic three-dimensional interpretation: the essence of the coherent body technology is that the abnormal phenomena of outstanding correlation and irrelevance between traces are calculated by utilizing seismic information, and a corresponding and effective calculation method is selected from a three-dimensional data body to realize the conversion of the coherent data body, so that the fracture characteristics are shown and the plane distribution of a crack development zone is predicted;
(3) the darker area on the coherent body slice represents lower coherence degree, the brighter area represents higher coherence degree, because the reflected wave continuity at the fault is poor, the fault can be explained on the coherent body slice according to the dark line or curve-shaped characteristic, at the initial stage of the explanation, time slice explanation can be made at intervals of large time-point, because the large fault extends far in the longitudinal direction and has reaction in the range of 10-100ms, the large fault can be firstly explained and then gradually encrypted until one horizontal time slice is formed every 1ms, so that after the explanation of the small fault is recognized on the slice, the fault can project a breakpoint at the corresponding position on the cross section, and the main measuring line, the contact measuring line and various arbitrary cross sections can be used for checking and correcting the explained breakpoint on the cross section, the combination of the breakpoints is carried out, so that small faults which are not easy to directly observe on the section are solved;
(4) using all the information of the data to automatically track the layer group to form three-dimensional data of the horizon surface, and adopting automatic superposition processing of surface slices or combination of human-computer interaction surface slices to form t according to the multi-solution degree of the horizon0Time, three-dimensional visual inspection, projection to vertical section inspection and modification to form t0Color pixel structure diagram, for t0Time or t0The color pixel map carries out other conventional interpretation work, and the target fault track lines on the seismic main section and the contact section are corrected and interpreted.
Furthermore, in the conversion process of the coherent data volume, an estimated value of three-dimensional seismic correlation can be obtained by calculating the local waveform similarity in the longitudinal direction and the transverse direction, seismic traces in a small range cut by a fault plane usually have different characteristics from adjacent traces, so that the local trace-to-trace correlation is suddenly changed, the correlation value on each grid point is calculated along a time slice, a profile with a low correlation value along the fault can be obtained, the process is repeated for a series of time slices, the profiles become the sections through the extraction of the three-dimensional correlation attribute body, and the three-dimensional reflection amplitude data volume can be converted into the three-dimensional correlation coefficient data volume.
Further, the coherent body algorithm has C1,C2,C3A total of 3 algorithms, C1The correlation between the transverse measuring line and the adjacent longitudinal measuring line of each channel is calculated by a coherent algorithm, two results are mixed and standardized by energy, and the essence is the expansion of the two correlation algorithms; c2The coherent algorithm is to carry out coherent calculation on any multi-channel seismic data based on a similar coherent algorithm; c1The coherent algorithm is implemented by means of a covariance matrix, since C2Coherent algorithms can greatly suppress noise when the effective signal is greater than the average of the noise, so C3Coherent algorithmCompared with other two coherent algorithms, the method has higher horizontal resolution and vertical resolution in fault identification and edge detection.
While the invention has been described above with reference to an embodiment, various modifications may be made and equivalents may be substituted for elements thereof without departing from the scope of the invention. In particular, the various features of the embodiments disclosed herein may be used in any combination, provided that there is no structural conflict, and the combinations are not exhaustively described in this specification merely for the sake of brevity and conservation of resources. Therefore, it is intended that the invention not be limited to the particular embodiments disclosed, but that the invention will include all embodiments falling within the scope of the appended claims.
Claims (3)
1. A method for seismic interpretation of a fault, the method comprising the steps of:
(1) fault section interpretation: the method is characterized in that a color variable density section is adopted to identify small faults, some small faults may have distorted wave crests but broken wave troughs in a section map, the fault is on a conventional section, people are used to pay attention to wave crest characteristic information and often ignore wave trough information, so that the situation is difficult to judge whether the situation is caused by the fault or lithological change, and the wave crests and the wave troughs are filled with different colors in the color variable density section, because human eyes are sensitive to the identification of the colors compared with the tiny change of waveforms, when the information of the wave crests and the wave troughs is displayed in front of an interpreter in a color form, the identification of the interpreter is easier, and the information can be used as a basis for identifying the small faults;
(2) tomographic three-dimensional interpretation: the essence of the coherent body technology is that the abnormal phenomena of outstanding and irrelevant relativity between traces are calculated by utilizing seismic information, and a corresponding and effective calculation method is selected from a three-dimensional data body to realize the conversion of the coherent data body, so that the fracture characteristics are shown and the plane distribution of a crack development zone is predicted;
(3) the darker area on the coherent body slice represents lower coherence degree, the brighter area represents higher coherence degree, because the reflected wave continuity at the fault is poor, the fault can be explained on the coherent body slice according to the dark line or curve-shaped characteristic, at the initial stage of the explanation, time slice explanation can be made at intervals of large time and small time, because the large fault extends far in the longitudinal direction and reacts in the range of 10-100ms, the large fault can be firstly explained and then gradually encrypted until one horizontal time slice is obtained every 1ms, so that after the explanation of the small fault is recognized on the slice, the fault can project a breakpoint at the corresponding position on the section, and the main measuring line, the contact measuring line and various arbitrary sections can be used for checking and correcting the explained breakpoint on the section, the combination of the breakpoints is carried out, so that small faults which are not easy to directly observe on the section are solved;
(4) using all the information of the data to automatically track the layer group to form three-dimensional data of the horizon surface, and adopting automatic superposition processing of surface slices or combination of human-computer interaction surface slices to form t according to the multi-solution degree of the horizon0Time, three-dimensional visual inspection, projection to vertical section inspection and modification to form t0Color pixel structure diagram, for t0Time or t0The color pixel map carries out other conventional interpretation work, and the target fault track lines on the seismic main section and the contact section are corrected and interpreted.
2. A method of seismic interpretation of faults according to claim 1, wherein: in the conversion process of the coherent data volume, an estimated value of three-dimensional seismic correlation can be obtained by calculating the waveform similarity of longitudinal and transverse parts, seismic channels in a small range cut by a fault plane usually have different characteristics from adjacent channels, so that the sudden change of the correlation between the local channels is caused, the correlation value on each grid point is calculated along a time slice, a profile with a low correlation value along the fault can be obtained, the process is repeated for a series of time slices, the profiles become the sections and are extracted through the three-dimensional correlation attribute body, and the three-dimensional reflection amplitude data volume can be converted into the three-dimensional correlation coefficient data volume.
3. A method of seismic interpretation of faults according to claim 2, wherein: the coherent body algorithm has C1,C2,C3A total of 3 algorithms, C1The correlation between the horizontal line and the adjacent vertical line of each channel is calculated by a coherent algorithm, two results are mixed and normalized by energy, and the essence is the expansion of the two correlation algorithms; c2The coherent algorithm is to carry out coherent calculation on any multi-channel seismic data based on a similar coherent algorithm; c1The coherent algorithm is implemented by means of a covariance matrix, since C2Coherent algorithms can greatly suppress noise when the effective signal is larger than the average value of the noise, so C3The coherent algorithm has higher horizontal resolution and vertical resolution than the other two coherent algorithms in fault identification and edge detection.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910977084.7A CN110646848A (en) | 2019-10-15 | 2019-10-15 | Seismic interpretation method for fault |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910977084.7A CN110646848A (en) | 2019-10-15 | 2019-10-15 | Seismic interpretation method for fault |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110646848A true CN110646848A (en) | 2020-01-03 |
Family
ID=68994065
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910977084.7A Pending CN110646848A (en) | 2019-10-15 | 2019-10-15 | Seismic interpretation method for fault |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110646848A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104597485A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Micro-fault detecting method and fault detecting device |
| CN107703542A (en) * | 2017-08-31 | 2018-02-16 | 中国石油天然气集团公司 | A kind of determination method and apparatus of seismic horizon |
| CN107765301A (en) * | 2017-10-13 | 2018-03-06 | 中国煤炭地质总局地球物理勘探研究院 | The method for quickly identifying and device of coal seam craven fault |
| CN109425899A (en) * | 2017-09-01 | 2019-03-05 | 中国石油天然气股份有限公司 | A kind of prediction technique and device of the distribution of carbonate rock fault belt |
| CN109557580A (en) * | 2017-09-27 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic interpretation method of tomography |
| US10338252B1 (en) * | 2009-06-01 | 2019-07-02 | Emerson Paradigm Holding Llc | Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface |
-
2019
- 2019-10-15 CN CN201910977084.7A patent/CN110646848A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10338252B1 (en) * | 2009-06-01 | 2019-07-02 | Emerson Paradigm Holding Llc | Systems and processes for building multiple equiprobable coherent geometrical models of the subsurface |
| CN104597485A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Micro-fault detecting method and fault detecting device |
| CN107703542A (en) * | 2017-08-31 | 2018-02-16 | 中国石油天然气集团公司 | A kind of determination method and apparatus of seismic horizon |
| CN109425899A (en) * | 2017-09-01 | 2019-03-05 | 中国石油天然气股份有限公司 | A kind of prediction technique and device of the distribution of carbonate rock fault belt |
| CN109557580A (en) * | 2017-09-27 | 2019-04-02 | 中国石油化工股份有限公司 | A kind of seismic interpretation method of tomography |
| CN107765301A (en) * | 2017-10-13 | 2018-03-06 | 中国煤炭地质总局地球物理勘探研究院 | The method for quickly identifying and device of coal seam craven fault |
Non-Patent Citations (2)
| Title |
|---|
| 杨瑞召等: "煤田三维地震采区小断层解释方法及应用 ", 《矿业工程研究》 * |
| 汤祥武等: "三维地震勘探小断层的解释方法及应用 ", 《煤矿开采》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cheng et al. | A zoning model for coal mining-induced strata movement based on microseismic monitoring | |
| CN101770038B (en) | Intelligent positioning method of mine microquake sources | |
| CN102385070B (en) | Radio wave perspective captive test (CT) method for overlong working surface | |
| Zhou et al. | Formation process and mechanism of seepage channels around grout curtain from microseismic monitoring: a case study of Zhangmatun iron mine, China | |
| CN105334548A (en) | Geological forecasting method for tunnel construction in karst area | |
| Liu et al. | Activation characteristics analysis on concealed fault in the excavating coal roadway based on microseismic monitoring technique | |
| CN105700010B (en) | Colliery joint earthquake holographic forecast method | |
| CN104766242A (en) | Method for evaluating dangerousness of water inrush from coal floor | |
| CN103726844A (en) | Automatic coal mining method based on coal mining work surface | |
| Roach et al. | Initial 4D seismic results after CO2 injection start-up at the Aquistore storage site | |
| CN110609335A (en) | Multi-means-based residual mining area complex condition detection method | |
| CN105891896A (en) | Feature information recognition and analysis method for underground mined area | |
| Ikuma | Maintenance of the undersea section of the Seikan Tunnel | |
| Jian et al. | Determining areas in an inclined coal seam floor prone to water-inrush by micro-seismic monitoring | |
| Zhang et al. | Vibration events in underground heading face and useful index for rock burst monitoring | |
| CN103344987A (en) | Method for explaining active fault moving strength and growing directions through artificial earthquake | |
| CN110646848A (en) | Seismic interpretation method for fault | |
| CN107784437B (en) | Stress concentration-based coal roadway driving face outburst danger identification method | |
| Tu et al. | Research situation and prospect of fully mechanized mining technology in thick coal seams in China | |
| Hu et al. | Imaging mining hazards within coalbeds using prestack wave equation migration of in-seam seismic survey data: A feasibility study with synthetic data | |
| CN104297799B (en) | High-speed layer top interface location determines method and depth determination method | |
| Zhang et al. | A case study of floor failure characteristics under fully mechanised caving mining conditions in extra-thick coal seams | |
| CN102221707A (en) | Seismic attribute extracting method for geologic structure containing over-thrust fault | |
| Wang et al. | Advanced prediction for multiple disaster sources of laneway under complicated geological conditions | |
| CN106680892A (en) | Exploration improvement method in mineral mining process |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200103 |