CN108535779A - A kind of optimization method for seismic exploration data - Google Patents
A kind of optimization method for seismic exploration data Download PDFInfo
- Publication number
- CN108535779A CN108535779A CN201810295757.6A CN201810295757A CN108535779A CN 108535779 A CN108535779 A CN 108535779A CN 201810295757 A CN201810295757 A CN 201810295757A CN 108535779 A CN108535779 A CN 108535779A
- Authority
- CN
- China
- Prior art keywords
- data
- time migration
- crp
- numeric field
- stack time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 230000003068 static effect Effects 0.000 claims abstract description 46
- 238000012937 correction Methods 0.000 claims abstract description 40
- 238000013508 migration Methods 0.000 claims abstract description 38
- 230000005012 migration Effects 0.000 claims abstract description 38
- 241001269238 Data Species 0.000 claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000010276 construction Methods 0.000 claims description 8
- 230000002547 anomalous effect Effects 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012795 verification Methods 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
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 present invention relates to a kind of optimization method for seismic exploration data, step includes carrying out fine processing to original seismic exploration data and obtaining CDP numeric field datas；1. and CRP numeric field datas pre-stack time migration processing is carried out to CDP numeric field datas and obtains pre-stack time migration data；Remaining static correction and pre-stack time migration processing are carried out to CDP numeric field data sequences, obtain pre-stack time migration data 2.；From pre-stack time migration data 1. and 2. in choose interface, 1. 2. the pre-stack time migration data more than the interface are matched into splicing with positioned at interface pre-stack time migration data below, obtain final seismic exploration data.The optimization method for seismic exploration data of the present invention can eliminate stratum depths hidden anomaly body to the distortion caused by back wave structural configuration, and then substantially improve the precision of seismic migration imaging.
Description
Technical field
The invention belongs to seismic exploration technique fields, and in particular to a kind of optimization method for seismic exploration data.
Background technology
Seismic prospecting is the difference using underground medium elasticity and density, by observing and analyzing the earth to artificial excitation
The geophysical exploration method of nature of ground and form is heuristically descended in the response of seismic wave.In the method for many seismic prospectings, fold
Preceding time migration technology has become to improve the important means of data imaging precision, can make data using prestack time migration technique
Correct playback, keeps the contact relation on tomography, stratum apparent, to make underground complicated structure obtain more accurately being imaged.
But for the hidden anomalies body such as igneous lithofacies existing for earth formation deep (such as 1000m or less), due to its presence
The cross directional variations of speed and thickness etc. so that seismic wave required hourage when passing through changes, to make
The back wave that underlying strata is reflected below igneous lithofacies all by (i.e. the time difference) the phenomenon that time of occurrence delay or time advance, is compeled
These back wave forms are made to generate distortion.Such time difference generation is directly related with hidden anomaly body, on same place over the ground
The influence at lower different depth interface is the same, and belongs to the static time difference, in order to distinguish other static time differences in this application by
Referred to as " the remaining static state time difference ", but utilize static correction known near surface static correction and earth's surface-consistent residual static corrections etc.
Method can not extract the remaining static time difference at all, and then can not eliminate the hidden anomaly body to caused by back wave structural configuration
Distortion, influences the precision of seismic migration imaging.
Invention content
In order to solve the problems, such as that above-mentioned all or part, the present invention provide a kind of optimization method for seismic exploration data,
It can eliminate stratum depths hidden anomaly body to the distortion caused by back wave structural configuration, and then substantially improve seismic migration
The precision of imaging.
The present invention provides a kind of optimization method for seismic exploration data, step includes：Step 1, to primitively
Seismic exploration data carry out fine processing, obtain CDP numeric field datas, wherein the fine processing includes at least near surface static correction and ground
Table consistency residual static correction is handled；Step 2, pre-stack time migration processing is carried out to the CDP numeric field datas, obtains prestack time
1. and CRP numeric field datas offset data；Step 3, CDP numeric field datas sequence is carried out at remaining static correction and pre-stack time migration
2. reason obtains pre-stack time migration data；Step 4, from the pre-stack time migration data 1. with pre-stack time migration data 2.
There is the stratum where hidden anomaly body in middle acquisition, and chooses interface from the stratum where the hidden anomaly body, by position
Pre-stack time migration data more than the interface 1. with positioned at the interface pre-stack time migration number below
According to 2. matching splicing, final seismic exploration data is obtained.Wherein, the remaining static correction includes the following steps：Step 3.1,
From the seismic channel data that offset distance is zero in each CRP trace gathers is extracted in the CRP numeric field datas, to set up total zero-offset road
Collect data；Step 3.2 is thrown oneself on the ground under at least two be total in zero-offset trace gather data below identification hidden anomaly body
Layer, and each CRP is therefrom obtained in the CRP numeric field datas in aforementioned at least two underlying strata caused by hidden anomaly body
The remnants static time differences；The remnants of each CRP static time differences are divided into two parts by step 3.3, and respectively as correspondence
Then the shot correction amount and geophone station correcting value of seismic channel use the shot correction amount and detection point calibration to the CDP
Numeric field data carries out static correction.
Further, the step 3.2 includes：Step 3.2.1, it is identified in the zero-offset trace gather data altogether hidden
At least two underlying stratas below anomalous body, and based on seismic interpretation technique middle level position pick-up method pickup aforementioned at least two
The layer position data of underlying strata；Step 3.2.2, obtain the intrinsic construction data of each layer from each layer position data, and by its
It is removed from each layer position data and obtains the remnants static time differences of each CRP in the CRP numeric field datas.
Further, in step 3.2.2, each layer is obtained from each layer position data with exponential smoothing or trend analysis method
Intrinsic construction data.
Further, in step 3.3, recycle earth's surface-consistent remaining using shot correction amount and geophone station correcting value
The method of static correction carries out static correction to CRP numeric field datas.
Further, the interface is adjacent with the bottom of hidden anomaly body or virtual level that is being in contact position interface.
The optimization method for seismic exploration data of the present invention can eliminate stratum depths hidden anomaly body to back wave
Distortion caused by structural configuration avoids occurring false construction in seismic profile, substantially improves the precision of seismic migration imaging.This
Outside, the optimization method for seismic exploration data of the invention is in addition to being suitable for there is transverse gage and speed to become stratum inside
Except the influence of the igneous lithofacies of change is eliminated, it is also applied for the shadow to existing low-velocity zone under mountain front near surface
Sound is eliminated, and is applied also for the Velocity anomalies that are generated by natural gas stack effect below seabed in marine seismic prospectiong
Influence is eliminated.
Description of the drawings
It, below will be to specific in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Embodiment or attached drawing needed to be used in the description of the prior art are briefly described.In all the appended drawings, similar element
Or part is generally identified by similar reference numeral.In attached drawing, each element or part might not be drawn according to actual ratio.
Fig. 1 is the flow chart of the optimization method for seismic exploration data of the embodiment of the present invention.
Fig. 2 is the flow chart of the remnants static corrections in the optimization method；
Fig. 3 is the migrated section partial enlarged view handled without above-mentioned optimization method；
Fig. 4 is the migrated section partial enlarged view for having been subjected to above-mentioned optimization method processing.
Specific implementation mode
The embodiment of technical solution of the present invention is described in detail below in conjunction with attached drawing.Following embodiment is only used for
Clearly illustrate technical scheme of the present invention, therefore be only used as example, and the protection model of the present invention cannot be limited with this
It encloses.
Fig. 1 shows the flow chart of the optimization method for seismic exploration data of the embodiment of the present invention.As shown in Figure 1,
The optimization method for seismic exploration data of the embodiment of the present invention specifically includes：Step S1, to original seismic exploration data into
Row fine processing obtains CDP numeric field datas, and wherein the fine processing includes at least near surface static correction and earth's surface-consistent residue is quiet
Correction；1. and the domains CRP number step S2 carries out pre-stack time migration processing to CDP numeric field datas, obtains pre-stack time migration data
According to；Step S3 carries out remaining static correction to CDP numeric field data sequences and pre-stack time migration is handled, obtains pre-stack time migration number
According to 2.；Step S4, from pre-stack time migration data, 1. 2. there is hidden anomaly body place in middle obtain with pre-stack time migration data
Stratum, and interface is chosen from the stratum where hidden anomaly body, by the pre-stack time migration more than the interface
1. 2. data match splicing with positioned at interface pre-stack time migration data below, obtain final seismic exploration data.
In fact, basic using static correcting method known near surface static correction and earth's surface-consistent residual static corrections etc.
It can not extract the remaining static time difference, but present inventor after long felt by having found below hidden anomaly body extremely
All there is identical change, the i.e. feature of multilayer bit integrity when the back wave approach hidden anomaly body of few two underlying stratas,
Also it can only be embodied in total zero-offset trace gather data, it is possible to effectively be sought in total zero-offset trace gather data residual
The remaining static state time difference.As shown in Fig. 2, the remnants static corrections include the following steps：Step S3.1, it is extracted from CRP numeric field datas each
The seismic channel data that offset distance is zero in CRP trace gathers, to set up total zero-offset trace gather data；Step S3.2, from total zero bias
It moves away from least two underlying stratas identified in trace gather data below hidden anomaly body, and therefrom obtains each in CRP numeric field datas
Remnants of the CRP in aforementioned at least two underlying strata caused by the hidden anomaly body static time differences；Step 3.3, by each CRP
The remnants static time differences be divided into two parts, and respectively as the shot correction amount and geophone station correcting value of corresponding seismic channel, so
Shot correction amount and detection point calibration is used to carry out static correction to CRP numeric field datas afterwards.It should be noted that step 3 belonging to above-mentioned is removed
Include except abovementioned steps 3.1 to step 3.3 further include each sub-steps for executing pre-stack time migration and handling, but
Because pre-stack time migration processing belongs to conventional technical means well known to those skilled in the art, so for the sake of saving length not
It repeats again.
It is hidden that the optimization method for seismic exploration data of the embodiment of the present invention can eliminate stratum depths (1000m or less)
Anomalous body is covered to the distortion caused by back wave structural configuration, is avoided occurring false construction in seismic profile, can be substantially improved
The precision of seismic migration imaging, refers to Fig. 3 and Fig. 4.In addition, the optimization method for seismic exploration data of the present invention is in addition to suitable
Except there is the influence of the igneous lithofacies of transverse gage and velocity variations to eliminate stratum inside, it is also applied for mountain
The preceding influence with existing low-velocity zone under near surface is eliminated, and is applied also for in marine seismic prospectiong below seabed
The influence of the Velocity anomalies generated by natural gas stack effect is eliminated.
The step S3.2 of the present embodiment can specifically be executed according to following sub-step：Step S3.2.1, in total zero-offset
At least two underlying stratas below hidden anomaly body are identified in trace gather data, and based on horizon picking side in seismic interpretation technique
Method picks up the layer position data of aforementioned at least two underlying strata；Step S3.2.2, by exponential smoothing, trend analysis method or other often
Rule extractive technique obtains the intrinsic construction data of each layer from the data of each layer position, and it is removed from each layer of position data and is obtained
To the remnants static time differences of each CRP in CRP numeric field datas.
In the present embodiment, in step S3.3, recycle earth's surface consistent with geophone station correcting value using shot correction amount
Property residual static correction method to CRP numeric field datas carry out static correction.By experimental verification, in addition to quiet using earth's surface-consistent residue
Except correction, the method that other static corrections can also be used, but utilize earth's surface-consistent residual static correction embodiment remnants quiet
Correction can more precisely, preferably eliminate stratum depths hidden anomaly body at least two layers of underlying strata back wave structure below
The distortion of appearance state.
In fact, the interface that step 3 is related to can be selected as by hidden anomaly body or without its virtual level position
Interface.But in the present embodiment, interface is preferably (± 10ms, or ± 10 be multiplied by ground adjacent with the bottom of hidden anomaly body
Total of the seismic wave in the speed v of underground propagation) or the virtual level position interface that is in contact, this can ensure final seismic prospecting number
For being imaged the data of true hidden anomaly body and without remaining static correction in, it is ensured that Velocity anomalies can be by more precisely
Imaging.
In conclusion can to eliminate stratum depths hidden for the optimization method for seismic exploration data of the embodiment of the present invention
Anomalous body avoids occurring false construction in seismic profile, it is inclined to substantially improve earthquake to the distortion caused by back wave structural configuration
Move into the precision of picture.
Finally it should be noted that：The above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；To the greatest extent
Present invention has been described in detail with reference to the aforementioned embodiments for pipe, it will be understood by those of ordinary skill in the art that：Its according to
So can with technical scheme described in the above embodiments is modified, either to which part or all technical features into
Row equivalent replacement；And these modifications or replacements, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme should all cover in the claim of the present invention and the range of specification.Especially, as long as there is no knots
Structure conflict, items technical characteristic mentioned in the various embodiments can be combined in any way.The present invention does not limit to
In specific embodiment disclosed herein, but include all technical solutions fallen within the scope of the appended claims.
Claims (5)
1. a kind of optimization method for seismic exploration data, which is characterized in that its step includes：
Step 1, fine processing is carried out to original seismic exploration data, CDP numeric field datas is obtained, wherein the fine processing is at least wrapped
Include near surface static correction and earth's surface-consistent residual static correction；
1. and the domains CRP number step 2, pre-stack time migration processing is carried out to the CDP numeric field datas, obtains pre-stack time migration data
According to；
Step 3, remaining static correction is carried out to CDP numeric field datas sequence and pre-stack time migration is handled, it is inclined to obtain prestack time
Move data 2.；
Step 4, from the pre-stack time migration data, 1. 2. there is hidden anomaly body institute in middle obtain with pre-stack time migration data
Stratum, and choose interface from the stratum where the hidden anomaly body, will be more than the interface it is described folded
1. 2. preceding time migration data match splicing with positioned at the interface pre-stack time migration data below, obtain final
Seismic exploration data；
Wherein, the remaining static correction includes the following steps：
Step 3.1, from the seismic channel data that offset distance is zero in each CRP trace gathers is extracted in the CRP numeric field datas, to set up
Zero-offset trace gather data altogether；
Step 3.2, from least two underlying stratas in the zero-offset trace gather data altogether below identification hidden anomaly body, and
It is residual caused by hidden anomaly body in aforementioned at least two underlying strata therefrom to obtain each CRP in the CRP numeric field datas
The remaining static state time difference；
The remnants of each CRP static time differences are divided into two parts by step 3.3, and respectively as the big gun of corresponding seismic channel
Then point calibration amount and geophone station correcting value carry out the CDP numeric field datas using the shot correction amount and detection point calibration
Static correction.
2. optimization method according to claim 1, which is characterized in that the step 3.2 includes：
Step 3.2.1, at least two underlying stratas in the zero-offset trace gather data altogether below identification hidden anomaly body,
And the layer position data of aforementioned at least two underlying strata are picked up based on seismic interpretation technique middle level position pick-up method；
Step 3.2.2, the intrinsic construction data of each layer are obtained from each layer position data, and by it from each layer digit
The remnants static time differences of each CRP in the CRP numeric field datas are obtained according to middle removal.
3. optimization method according to claim 2, which is characterized in that in step 3.2.2, with exponential smoothing or trend analysis
Method obtains the intrinsic construction data of each layer from each layer position data.
4. optimization method according to any one of claim 1-3, which is characterized in that in step 3.3, use shot point school
The method that positive quantity and geophone station correcting value recycle earth's surface-consistent residual static correction carries out static correction to CRP numeric field datas.
5. optimization method according to any one of claim 1-3, which is characterized in that the interface be with it is described hidden
The bottom of anomalous body is adjacent or the virtual level position interface that is in contact.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810295757.6A CN108535779B (en) | 2018-03-30 | 2018-03-30 | A kind of optimization method for seismic exploration data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810295757.6A CN108535779B (en) | 2018-03-30 | 2018-03-30 | A kind of optimization method for seismic exploration data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108535779A true CN108535779A (en) | 2018-09-14 |
| CN108535779B CN108535779B (en) | 2019-07-23 |
Family
ID=63482682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810295757.6A Active CN108535779B (en) | 2018-03-30 | 2018-03-30 | A kind of optimization method for seismic exploration data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108535779B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5189643A (en) * | 1992-03-05 | 1993-02-23 | Conoco Inc. | Method of accurate fault location using common reflection point gathers |
| US20020042678A1 (en) * | 2000-08-07 | 2002-04-11 | Dimitri Bevc | Typing picks to horizons in migration velocity analysis |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN106094035A (en) * | 2016-08-18 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of residual static correction computational methods quickly processed based on unit |
-
2018
- 2018-03-30 CN CN201810295757.6A patent/CN108535779B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5189643A (en) * | 1992-03-05 | 1993-02-23 | Conoco Inc. | Method of accurate fault location using common reflection point gathers |
| US20020042678A1 (en) * | 2000-08-07 | 2002-04-11 | Dimitri Bevc | Typing picks to horizons in migration velocity analysis |
| CN101315427A (en) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | Method and system for processing seismic exploration data of complex area |
| CN106094035A (en) * | 2016-08-18 | 2016-11-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of residual static correction computational methods quickly processed based on unit |
Non-Patent Citations (5)
| Title |
|---|
| JIANBO FEI，ET AL.: "The solution of complex near-surface problem using integrated velocity model building in Southern Tianshan Foothills", 《SEG》 * |
| WOOK B. LEE AND LIN ZHANG: "Residual shot profile migration", 《GEOPHYSICS》 * |
| 周俊，等: "塔北油田地震资料的处理技术", 《桂林理工大学学报》 * |
| 朱景修，等: "泌阳凹陷陡坡带三维地震资料叠前时间偏移连片处理实践", 《石油地质与工程》 * |
| 费建博，杨子川: "双界面匹配一体化速度建模技术研究与应用-以天山南山前带阳霞区块为例", 《石油物探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108535779B (en) | 2019-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Malehmir et al. | 3D seismic reflection imaging of volcanic-hosted massive sulfide deposits: Insights from reprocessing Halfmile Lake data, New Brunswick, Canada | |
| US7952960B2 (en) | Seismic imaging with natural Green's functions derived from VSP data | |
| CN102176054B (en) | Near-surface comprehensive information processing explanation method | |
| Marsden | Static corrections—A review | |
| Juhojuntti et al. | 3D seismic survey at the Millennium uranium deposit, Saskatchewan, Canada: Mapping depth to basement and imaging post-Athabasca structure near the orebody | |
| Bradford et al. | Depth characterization of shallow aquifers with seismic reflection, Part II—Prestack depth migration and field examples | |
| Yordkayhun et al. | 3D seismic traveltime tomography imaging of the shallow subsurface at the CO2SINK project site, Ketzin, Germany | |
| Bradford et al. | Imaging complex structure in shallow seismic-reflection data using prestack depth migration | |
| Guo et al. | Becoming effective velocity-model builders and depth imagers, Part 2—The basics of velocity-model building, examples and discussions | |
| CN105425300B (en) | A kind of residual static corrections | |
| CN108535779B (en) | A kind of optimization method for seismic exploration data | |
| Smith et al. | Making seismic monitoring work in a complex desert environment—4D processing | |
| Landrø et al. | Using diving waves for detecting shallow overburden gas layers | |
| Wang | Full-azimuth, high-density, 3D single-point seismic survey for shale gas exploration in a loess plateau area, Southeast of Ordos Basin, China | |
| Palmer | Integrating long and short wavelength time and amplitude statics | |
| Ghazali et al. | Velocity model building and imaging in the presence of shallow gas | |
| CN117930359B (en) | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields | |
| Cook | Deep basement seismic reflection profiling of the Purcell anticlinorium using a land air gun source | |
| Osinowo | Reprocessing of regional 2D marine seismic data of part of Taranaki basin, New Zealand using Latest processing techniques | |
| Zhang | Analysis of factors affecting residual moveout picking and solutions | |
| Moreno et al. | Solution of high velocity anomalies imperceptible to seismic resolution, by means of synthetic models, Penobscot Field, Canada | |
| Vigh et al. | An offshore Gulf of Mexico case study applying full-waveform inversion | |
| Ogunsuyi et al. | Integrating seismic-velocity tomograms and seismic imaging: Application to the study of a buried valley | |
| Simarmata et al. | Estimation of anisotropy parameters VTI (vertical transverse isotropy) using velocity variation with offset (VVO) method, case study: BS oil field | |
| Yang et al. | Velocity model building with a time shift objective function: Gulf of Mexico examples |
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 |