CN101950034B - 一种用于地震反演的自然电位曲线基线漂移校正方法 - Google Patents
一种用于地震反演的自然电位曲线基线漂移校正方法 Download PDFInfo
- Publication number
- CN101950034B CN101950034B CN201010254765XA CN201010254765A CN101950034B CN 101950034 B CN101950034 B CN 101950034B CN 201010254765X A CN201010254765X A CN 201010254765XA CN 201010254765 A CN201010254765 A CN 201010254765A CN 101950034 B CN101950034 B CN 101950034B
- Authority
- CN
- China
- Prior art keywords
- baseline
- data
- spontaneous potential
- curve
- mud stone
- 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
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于地震反演的自然电位曲线基线漂移校正方法，包括：获取砂岩段的顶深和厚度数据；生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成泥岩段对应的采样点的深度值；对自然电位测井曲线数据进行检索，确定泥岩段的各采样点深度值对应的自然电位幅度值，生成泥岩段的采样点深度值与自然电位幅度值对应的数据点；将各个泥岩段对应的数据点相连，生成自然电位曲线基线；将自然电位曲线与基线做校正处理，生成校正后的自然电位曲线。本发明实施例能够自动计算并校正自然电位曲线基线，具有较强的准确性，而且能够进行大量井数据的自动批量化处理，有利于提高地震反演测井曲线校正的精度和效率。
Description
技术领域
本发明涉及油藏地球物理开发地震反演所用测井资料的处理技术，尤其涉及一种用于地震反演的自然电位曲线基线漂移校正方法及装置。
背景技术
现有的声波测井曲线不能刻画岩性是砂泥岩储层类型油田普遍存在的现象，为此，在开展地震反演的过程中有必要重构储层物理特征曲线。但是，由于老油田密度曲线和伽马曲线测量较少，基础井网基本没有这两种测井曲线，而且电阻率类测井曲线，如微电极曲线、微电位曲线等受流体变化的影响较大，鉴于自然电位曲线能够反映储层泥质含量的变化，与孔隙度有一定关系，为此采用自然电位曲线开展用于地震反演的储层物理特征曲线重构是很好的方法。但是自然电位曲线存在严重的基线漂移现象，因此需要开展自然电位曲线基线漂移校正。
目前现有技术主要采取人工绘制基线方法进行校正，即先在软件窗口中显示自然电位曲线，然后用鼠标点击测井曲线上泥岩段上的点，点击多个点并进行连接即为自然电位的基线，然后再根据相关公式进行计算，校正掉自然电位的基线漂移量。
现有技术的这种采取在软件窗口中人工绘制自然电位基线的校正方法，能够满足油田勘探阶段已知井较少情况下的地震反演工作，但是对于高含水油田开展的开发地震工作而言，由于一个开发区块至少有上百口井，有的几百口井，甚至上千口井，这种每口井都要手工绘制基线的方法，工作效率非常低下，严重影响工作进度，从而需要研发适合高含水油田多井条件下的开发地震反演所需的快速自动测井曲线基线漂移校正方法。
发明内容
本发明实施例提供一种用于地震反演的自然电位曲线基线漂移校正方法，包括：获取砂岩段的顶深和厚度数据；根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值；根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点；将各个泥岩段对应的数据点相连，生成自然电位曲线基线；将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线；其中，将所述自然电位曲线与所述自然电位曲线基线做校正处理，包括：利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。。
本发明实施例还提供一种用于地震反演的自然电位曲线基线漂移校正装置，包括：数据获取单元，用于获取砂岩段的顶深和厚度数据；深度值生成单元，用于根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值；数据点生成单元，用于根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点；基线生成单元，用于将各个泥岩段对应的数据点相连，生成自然电位曲线基线；校正曲线生成单元，用于将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线，其包括：
利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。
本发明实施例的用于地震反演的自然电位曲线基线漂移校正方法及装置，能够自动计算自然电位曲线基线，具有较强的准确性，而且实现了多口井的自动批量化处理，大大提高了工作效率，节省了科技人员宝贵的人力。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1为本发明实施例的用于地震反演的自然电位曲线基线漂移校正方法的流程图；
图2为本发明实施例的用于地震反演的自然电位曲线基线漂移校正装置的结构示意图；
图3为采用本发明实施例的用于地震反演的自然电位基线漂移校正方法来确定的曲线基线的示意图；
图4为本发明实施例的用于地震反演的自然电位曲线基线校正后的对比示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例的用于地震反演的自然电位曲线基线漂移校正方法的流程图，如图所示，所述方法包括：
S101，获取砂岩段的顶深和厚度数据。本实施例中，可从油田开发单井小层数据表中来获取砂岩段的顶深和厚度数据。单井小层数据表为一个数据库文件，其结构主要包括：jh表示井号、yczmc表示油层组名称、xch表示小层号、xfch表示细分层号、syds表示砂岩顶深、syhd表示砂岩厚度等，如下表一所示。
表一：单井小层数据库
| jh | yczmc | xch | xfch | syds | syhd |
| D50 | S2 | 1 | 1 | 882.20 | 1.60 |
| D50 | S2 | 1 | 2 | 884.60 | 0.60 |
| D50 | S2 | 1 | 3 | 886.20 | 0.30 |
| D50 | S2 | 2 | 1 | 887.10 | 1.20 |
| D50 | S2 | 2 | 2 | 888.90 | 0.70 |
S102，根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值。根据步骤S101中的单井小层数据库所示，这个数据表仅记录砂岩的信息，而不是砂岩的层段就是泥岩，具体地讲就是在D50这口井上地下882.2m到883.8m(882.2+1.6)这一1.6m的深度范围的岩性是砂岩，而从883.8m(882.2+1.6)到884.6m这一深度范围就是泥岩。因此，泥岩段的深度值等于所述泥岩段对应的上层砂岩层的顶深数据加上厚度数据，即：
泥岩顶深＝syds1(其上砂岩的顶深)+syhd(其上砂岩的厚度)；
泥岩厚度＝syds2(其下砂岩的顶深)-泥岩顶深；
泥岩底深＝nyds(泥岩顶深)+nyhd(泥岩厚度)。
在泥岩顶深和泥岩底深范围内按着测井曲线采样间隔的深度值进行计算，既可以得到泥岩段的各个深度点。如测井曲线采样间隔通常为0.125m，那么下一个点的深度值就是泥岩顶深+0.125m。
S103，从根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点。当在步骤S102中获取到泥岩段的深度值后，即确定了泥岩段的位置后，到自然电位测井曲线上检索到这一段泥岩段对应的测井曲线，自然电位测井曲线由多个自然电位测井曲线数据组成，自然电位测井曲线数据是一个深度值对应一个自然电位曲线的自然电位幅度值，因此根据步骤S102中确定的深度值，可以得到对应的自然电位幅度值，即可以生成所述泥岩段深度值与所述自然电位幅度值的对应数据点。
S104，将各个泥岩段对应的数据点相连，生成自然电位曲线基线。本实施例中，此步骤具体为根据所述的数据点，在所述的自然电位曲线中相连，生成自然电位曲线基线。
S105，将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线并输出显示。在步骤S104确定了自然电位曲线基线后，利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。
对应于上述实现方法，图2为本发明实施例的用于地震反演的自然电位曲线基线漂移校正装置的结构示意图。如图所示，自然电位曲线基线漂移校正装置包括：
数据获取单元101，用于获取砂岩段的顶深和厚度数据。本实施例中，数据获取单元101可从油田开发单井小层数据表中来获取砂岩段的顶深和厚度数据。如上述表一所示，单井小层数据表为一个数据库文件，其结构主要包括：井号、油层组名称、小层号、细分层号、砂岩顶深、砂岩厚度等。
深度值生成单元102，用于根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值。根据单井小层数据表(表一)所示，这个数据表仅记录砂岩的信息，而不是砂岩的层段就是泥岩，具体地讲就是在表一中的D50这口井上地下882.2m到883.8m(882.2+1.6)这一1.6m的深度范围的岩性是砂岩，而从883.8m(882.2+1.6)到884.6m这一深度范围就是泥岩。因此，泥岩段的深度值等于所述泥岩段对应的上层砂岩层的顶深数据加上厚度数据，即：
泥岩顶深＝syds1(其上砂岩的顶深)+syhd(其上砂岩的厚度)；
泥岩厚度＝syds2(其下砂岩的顶深)-泥岩顶深；
泥岩底深＝nyds(泥岩顶深)+nyhd(泥岩厚度)。
在泥岩顶深和泥岩底深范围内按着测井曲线采样间隔的深度值进行计算，既可以得到泥岩段的各个深度点。如测井曲线采样间隔通常为0.125m，那么下一个点的深度值就是泥岩顶深+0.125m。
数据点生成单元103，用于根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点。当在深度值生成单元102获取到泥岩段的深度值后，即确定了泥岩段的位置后，到自然电位测井曲线上检索到这一段泥岩段对应的测井曲线，自然电位测井曲线由多个自然电位测井曲线数据组成，自然电位测井曲线数据是一个深度值对应一个自然电位曲线的自然电位幅度值，因此根据深度值生成单元102确定的深度值，可以得到对应的自然电位幅度值，即可以生成所述泥岩段深度值与所述自然电位幅度值的相对应的数据点。
基线生成单元104，用于将各个泥岩段对应的数据点相连，生成自然电位曲线基线。本实施例中，基线生成单元104在所述的自然电位曲线中将生成的数据点进行相连，生成的曲线即为自然电位曲线基线。
校正曲线生成单元105，用于将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线。
在基线生成单元104确定了自然电位曲线基线后，利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。
图3为采用本发明实施例的用于地震反演的自然电位基线漂移校正方法来确定的曲线基线的示意图。曲线1为根据本发明实施例确定的自然电位曲线基线，在曲线回返较小的薄泥岩段做忽略处理。
图4为本发明实施例的用于地震反演的自然电位曲线基线校正后的对比示意图。曲线2为未校正的自然电位曲线，曲线3为校正后的自然电位曲线，可以看到，根据本发明的方法及装置，自然电位曲线的基线漂移现象得到了很好的校正。
本发明实施例的自然电位曲线基线漂移校正方法及装置，能够自动计算自然电位曲线基线，具有较强的准确性，而且实现了多口井的自动批量化处理，大大提高了工作效率，节省了科技人员宝贵的人力。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种用于地震反演的自然电位曲线基线漂移校正方法，其特征在于，所述方法包括：
获取砂岩段的顶深和厚度数据；
根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值；
根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点；
将各个泥岩段对应的数据点相连，生成自然电位曲线基线；
将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线；
其中，将所述自然电位曲线与所述自然电位曲线基线做校正处理，包括：
利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。
2.如权利要求1所述的方法，其特征在于，所述获取砂岩段的顶深和厚度数据，包括：
从油田开发单井小层数据表中获取砂岩段的顶深和厚度数据。
3.如权利要求1所述的方法，其特征在于，所述根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，包括：
泥岩段的顶深数据等于所述泥岩段对应的上层砂岩层的顶深数据加上厚度数据。
4.一种用于地震反演的自然电位曲线基线漂移校正装置，其特征在于，所述装置包括：
数据获取单元，用于获取砂岩段的顶深和厚度数据；
深度值生成单元，用于根据所述砂岩段的顶深和厚度数据生成泥岩段的顶深和厚度数据，并按自然电位测井曲线的采样间隔生成所述泥岩段对应的采样点的深度值；
数据点生成单元，用于根据所述的泥岩段的采样点的深度值对所述自然电位测井曲线数据进行检索，确定所述泥岩段的各采样点深度值对应的自然电位幅度值，生成所述泥岩段的采样点深度值与自然电位幅度值对应的数据点；
基线生成单元，用于将各个泥岩段对应的数据点相连，生成自然电位曲线基线；
校正曲线生成单元，用于将所述自然电位曲线与所述自然电位曲线基线做校正处理，生成校正后的自然电位曲线，其包括：
利用公式SP_NEWi＝SPi-((Di-D1)/(Dn-D1))×(Vn-V1)-V1进行校正计算，其中，
SP_NEWi：第i采样点处基线漂移校正后的自然电位幅度值；
SPi：第i采样点处校正前的自然电位幅度值；
Di：第i采样点处的深度值；
D1：基线起点深度值；
V1：基线起点自然电位幅度值；
Dn：基线终点深度值；
Vn：基线终点自然电位幅度值。
5.如权利要求4所述的装置，其特征在于，所述的数据获取单元从油田开发单井小层数据表中获取砂岩段的顶深和厚度数据。
6.如权利要求4所述的装置，其特征在于，所述的深度值生成单元用于生成泥岩段的深度值，包括：
泥岩段的顶深数据等于所述泥岩段对应的上层砂岩层的顶深数据加上厚度数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010254765XA CN101950034B (zh) | 2010-08-16 | 2010-08-16 | 一种用于地震反演的自然电位曲线基线漂移校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010254765XA CN101950034B (zh) | 2010-08-16 | 2010-08-16 | 一种用于地震反演的自然电位曲线基线漂移校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101950034A CN101950034A (zh) | 2011-01-19 |
| CN101950034B true CN101950034B (zh) | 2011-11-09 |
Family
ID=43453597
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201010254765XA Active CN101950034B (zh) | 2010-08-16 | 2010-08-16 | 一种用于地震反演的自然电位曲线基线漂移校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101950034B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967881A (zh) * | 2012-11-12 | 2013-03-13 | 中国石油天然气股份有限公司 | 用于地震资料解释的地质小层数据深时转换方法及装置 |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562051B (zh) * | 2011-12-31 | 2015-07-15 | 中国石油集团川庆钻探工程有限公司 | 处理气测图版原始数据的方法 |
| CN104295293B (zh) * | 2014-10-23 | 2017-04-12 | 中国石油天然气股份有限公司 | 一种获取测井密度曲线的方法 |
| CN109143370B (zh) * | 2018-07-25 | 2020-03-31 | 中国地震局地球物理研究所 | 地震动加速度记录基线漂移的校正方法 |
| CN110488349B (zh) * | 2019-08-20 | 2021-11-02 | 福建省建筑设计研究院有限公司 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6473696B1 (en) * | 2001-03-13 | 2002-10-29 | Conoco Inc. | Method and process for prediction of subsurface fluid and rock pressures in the earth |
| US7529626B1 (en) * | 2004-05-14 | 2009-05-05 | Leroy Ellis | Method of integration and displaying of information derived from a mud gas isotope logging interpretative process in association with geophysical and other logs from oil and gas drilling operations |
| WO2006113502A2 (en) * | 2005-04-14 | 2006-10-26 | Swanson Consulting, Inc. | Modeling clastic reservoirs |
| CN100562764C (zh) * | 2007-12-28 | 2009-11-25 | 西安西凌地球物理技术开发有限责任公司 | 自然电位连续剖面勘探方法 |
| CN101290357B (zh) * | 2008-06-13 | 2010-06-02 | 杨辉 | 基于小循环平面多极同步基点的地面自然电位数据采集处理方法 |
-
2010
- 2010-08-16 CN CN201010254765XA patent/CN101950034B/zh active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102967881A (zh) * | 2012-11-12 | 2013-03-13 | 中国石油天然气股份有限公司 | 用于地震资料解释的地质小层数据深时转换方法及装置 |
| CN102967881B (zh) * | 2012-11-12 | 2015-07-08 | 中国石油天然气股份有限公司 | 用于地震资料解释的地质小层数据深时转换方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101950034A (zh) | 2011-01-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101950034B (zh) | 一种用于地震反演的自然电位曲线基线漂移校正方法 | |
| CN109061765B (zh) | 非均质薄砂岩互层油藏的圈闭评价方法 | |
| CN102646141B (zh) | 一种非等深钻孔的地质剖面图自动成图方法 | |
| CN102879821B (zh) | 一种针对地震叠前道集的同相轴精细拉平处理方法 | |
| CN104677734A (zh) | 一种多尺度信息融合的致密地层脆性指数测量方法 | |
| CN102967881B (zh) | 用于地震资料解释的地质小层数据深时转换方法及装置 | |
| CN107515290B (zh) | 岩石矿物组分含量定量计算方法 | |
| CN109799540B (zh) | 基于地质信息约束下的火山岩型铀矿床磁化率反演方法 | |
| CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
| CN106646613A (zh) | 深度域多尺度井控建模与成像联合处理方法 | |
| CN110579802A (zh) | 一种天然气水合物储层物性参数的高精度反演方法 | |
| CN105588883A (zh) | 三维岩石力学参数获取方法和系统 | |
| CN104142516A (zh) | 一种薄单砂层厚度预测方法 | |
| CN115327616B (zh) | 一种海量数据驱动的矿山微震震源自动定位方法 | |
| CN105649616A (zh) | 一种低渗气藏井下节流条件下气井动储量评价方法 | |
| CN103852789B (zh) | 用于地震数据的非线性层析方法及其装置 | |
| CN108957554B (zh) | 一种地球物理勘探中的地震反演方法 | |
| CN112444856B (zh) | 一种基于深度学习的地震数据分辨率提高方法 | |
| CN102841382B (zh) | 模型约束静校正法的校正量计算方法 | |
| CN201747364U (zh) | 一种声波时差曲线重构设备 | |
| CN106842289B (zh) | 一种适用于测井约束反演的波阻抗曲线去压实处理方法 | |
| CN106285653A (zh) | 确定页岩气储层的有机孔隙度的方法 | |
| CN106761666A (zh) | 一种四探头散射伽马测井及非线性数据反演的方法及装置 | |
| CN107632137B (zh) | 一种铀矿勘查中土壤氡气瞬时测量数据校正方法 | |
| CN104570100B (zh) | 多子波克希霍夫地震数据偏移方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |