CN110794038B - 基于声波波速各向异性的岩心室内定向方法 - Google Patents
基于声波波速各向异性的岩心室内定向方法 Download PDFInfo
- Publication number
- CN110794038B CN110794038B CN201911026888.5A CN201911026888A CN110794038B CN 110794038 B CN110794038 B CN 110794038B CN 201911026888 A CN201911026888 A CN 201911026888A CN 110794038 B CN110794038 B CN 110794038B
- Authority
- CN
- China
- Prior art keywords
- core
- wave
- rock
- orientation
- anisotropy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000011161 development Methods 0.000 claims abstract description 14
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims description 47
- 238000012216 screening Methods 0.000 claims description 7
- 239000011229 interlayer Substances 0.000 claims description 6
- 238000007914 intraventricular administration Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 2
- 238000013508 migration Methods 0.000 abstract description 2
- 230000005012 migration Effects 0.000 abstract description 2
- 239000011707 mineral Substances 0.000 abstract description 2
- 239000000203 mixture Substances 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/02—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by mechanically taking samples of the soil
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/011—Velocity or travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
- G01N2291/0232—Glass, ceramics, concrete or stone
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Mining & Mineral Resources (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Geology (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Acoustics & Sound (AREA)
- Soil Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Signal Processing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于声波波速各向异性的岩心室内定向方法。本发明基于钻井岩心与取心段地层的声波波速各向异性的相一致，建立了矿场声波测井与室内多方位岩样声波测试相结合实现岩心地面定向的方法，本发明不仅操作简单、且不受岩样岩性及矿物组成的影响，同时还可对取自同一层段的岩心实现批量定位。本发明具有结果可靠、操作便捷、经济高效、适用性广的特点。本发明对地层构造(裂缝)发育解释、油气运移方向分析、开发方案合理优化、老油田剩余油开发井部署优化等具有重要的基础指导价值。
Description
技术领域
本发明涉及石油天然气钻探技术领域，具体涉及一种基于声波波速各向异性的岩心室内定向方法。
背景技术
钻井取心是研究深部地层地质特征信息的直接、重要载体。岩心定向是基于岩心研究地层的构造特征、沉积相、产状、裂缝发育特征以及地应力岩心测试评价的基础；在注水、注气开发的油田，岩心定向有助于科学直观地认识、分析油、气、水的流动方向，为井网部署、优化提供依据。因此，钻井岩心的科学定向对油气勘探开发具有重要意义。
然而，大量钻井岩心为非定向钻取,取心过程中经过复杂旋转,无法直接获的其原始方位,致使岩心的众多地质信息无法有效利用。从20世纪60年代逐渐发展起来的定向取心(机械法)技术、古地磁岩心定向以及基于成像测井图像岩心定向等方法、技术，虽然可以确定岩心方位，但应用实践表明，但这些方法都不同程度存在操作工艺复杂、经济成本高、分析结果可靠性差等问题。
发明内容
针对现有技术中的上述不足，本发明提供的一种基于声波波速各向异性的岩心室内定向方法解决了岩心定向存在操作工艺复杂、经济成本高、分析结果可靠性差的问题。
为了达到上述发明目的，本发明采用的技术方案为：一种基于声波波速各向异性的岩心室内定向方法，包括以下步骤：
S1、获取新鲜完整的圆柱形钻井岩心；
S2、对岩心进行表观结构观察，根据岩心的表观结构观察结果筛选无裂缝发育的岩心；
S3、沿径向标记岩心的相对0°方位，并设置测试方位间隔，从相对0°方位根据测试方位间隔沿顺时针方向依次确定测试方位；
S4、从相对0°方位开始，通过投射法按照测试方位对岩心进行沿直径方向的声波波速进行测试，得到不同测试方位下的声波波速；
S5、由不同测试方位下的声波波速分析波速与测试方位的关系，确定岩心的最大波速相对方位；
S6、基于声波测试井信息分析获取井周地层声波波速的各项异性特征，确定波速最大方位；
S7、将岩心的最大波速相对方位与波速最大方位进行对比，确定岩心的地理方位。
进一步地：所述步骤S1中圆柱形钻井岩心的直径为4in或2in，高度为2in～4in。
进一步地：所述步骤S2中的具体步骤为：观察岩心的层理及其他弱结构面的发育情况，筛选结构完整、层理均匀且无夹层、裂缝及其它结构面发育的岩心。
进一步地：所述步骤S3中测试方位间隔为5°、10°或15°。
进一步地：所述步骤S6中的声波测试井信息包括XMAC、SonicScan或DSI声波测试井信息。
进一步地：所述步骤S6中的井周地层为岩心对应的无断层、无夹层和无其他非均质结构地层，且该井周地层的厚度不超过25m。
进一步地：所述步骤S6中的波速最大方位的计算方法为：从声波测试井信息中提取快波数据和慢波数据，通过快波数据和慢波数据计算各向异性系数，根据各向异性系数统计层段的快波方位和慢波方位，根据快波方位和慢波方位分析目标层段的各向异性，并获取井周地层的波速最大方位。
进一步地：所述步骤S6中各向异性系数的计算公式为：
上式中，S为各向异性系数，Vfast为快波数据，VSlow为慢波数据。
进一步地：所述步骤S7中岩心的地理方位的确定方法为：令岩心的最大波速相对方位减去波速最大方位，得到岩心的地理方位。
本发明的有益效果为：本发明基于钻井岩心与取心段地层的声波波速各向异性的相一致，建立了矿场声波测井与室内多方位岩样声波测试相结合实现岩心地面定向的方法，本发明不仅操作简单、且不受岩样岩性及矿物组成的影响，同时还可对取自同一层段的岩心实现批量定位。本发明具有结果可靠、操作便捷、经济高效、适用性广的特点。本发明对地层构造(裂缝)发育解释、油气运移方向分析、开发方案合理优化、老油田剩余油开发井部署优化等具有重要的基础指导价值。
附图说明
图1为本发明流程图；
图2为本发明实施例中岩样波速随测试方位的变化及其最大波速方位、最小波速方位；
图3为本发明实施例中地层快波方位统计图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，一种基于声波波速各向异性的岩心室内定向方法，包括以下步骤：
S1、获取新鲜完整的圆柱形钻井岩心；圆柱形钻井岩心的直径为4in或2in，高度为2in～4in。
S2、对岩心进行表观结构观察，根据岩心的表观结构观察结果筛选无裂缝发育的岩心；具体步骤为：观察岩心的层理及其他弱结构面的发育情况，筛选结构完整、层理均匀且无夹层、裂缝及其它结构面发育的岩心。
S3、沿径向标记岩心的相对0°方位，并设置测试方位间隔，测试方位间隔为5°、10°或15°，从相对0°方位根据测试方位间隔沿顺时针方向依次确定测试方位；
S4、从相对0°方位开始，通过投射法按照测试方位对岩心进行沿直径方向的声波波速进行测试，得到不同测试方位下的声波波速；
S5、由不同测试方位下的声波波速分析波速与测试方位的关系，确定岩心的最大波速相对方位；
S6、基于声波测试井信息分析获取井周地层声波波速的各项异性特征，确定波速最大方位；声波测试井信息包括XMAC、SonicScan或DSI声波测试井信息。井周地层为岩心对应的无断层、无夹层和无其他非均质结构地层，且该井周地层的厚度不超过25m。从声波测试井信息中提取快波数据和慢波数据，通过快波数据和慢波数据计算各向异性系数，根据各向异性系数统计层段的快波方位和慢波方位，根据快波方位和慢波方位分析目标层段的各向异性，并获取井周地层的波速最大方位。
S7、将岩心的最大波速相对方位与波速最大方位进行对比，令岩心的最大波速相对方位减去波速最大方位，得到岩心的地理方位。
本发明的实施例为：
S1、对需定向的钻井取心进行尺寸测试，包括不同方位直径(至少两个正交方位)；
S2、获取矿场泥页岩岩样并对其进行观察描述，筛选结构完整、均匀的岩心段；
S3、对已筛选好的岩心，选定某一径向标定为相对0°方位，并沿顺时针方向隔一定角度标定一个方向直至标完一周(例如间隔30°，分别标定0°，30°，60°，90°，120°，150°，180°)。
S4、采用透射法，从相对0°方位开始，对结构完整的标记段岩心沿标定方位进行声波波速测试。
S5、由不同方位的声波波速测试结果，分析波速随测试方位的变化，确定岩心的最大波速相对方位与最小波速相对方位分别为90°、0°，如图2所示。
S6、利用DSI声波测井数据，分别提取快波、横波，计算分析波速各向异性，并统计快波方位，图所示快波方位为90°～105°，中心方位97.5°，如图3所示。
S7、将室内测试得到的岩心最大波速相对方位与测井数据分析得到的井周地层最大波速方位进行对比，确定岩心的地理方位，实现岩心的定向。
因此，所分析岩心最大波速相对方位90°对应着井下实际地理方位97.5°，即岩心相对方位-7.5°对应正北方向，相对方位82.5°对应正东方向，从而实现了该岩心的室内定向。
Claims (7)
1.一种基于声波波速各向异性的岩心室内定向方法，其特征在于，包括以下步骤：
S1、获取新鲜完整的圆柱形钻井岩心；
S2、对岩心进行表观结构观察，根据岩心的表观结构观察结果筛选无裂缝发育的岩心；
S3、沿径向标记岩心的相对0°方位，并设置测试方位间隔，从相对0°方位根据测试方位间隔沿顺时针方向依次确定测试方位；
S4、从相对0°方位开始，通过透射法按照测试方位对岩心进行沿直径方向的声波波速进行测试，得到不同测试方位下的声波波速；
S5、由不同测试方位下的声波波速分析波速与测试方位的关系，确定岩心的最大波速相对方位；
S6、基于声波测试井信息分析获取井周地层声波波速的各项异性特征，确定波速最大方位；
所述步骤S6中的波速最大方位的计算方法为：从声波测试井信息中提取快波数据和慢波数据，通过快波数据和慢波数据计算各向异性系数，根据各向异性系数统计层段的快波方位和慢波方位，根据快波方位和慢波方位分析目标层段的各向异性，并获取井周地层的波速最大方位；
各向异性系数的计算公式为：
上式中，S为各向异性系数，Vfast为快波数据，VSlow为慢波数据；
S7、将岩心的最大波速相对方位与波速最大方位进行对比，确定岩心的地理方位。
2.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S1中圆柱形钻井岩心的直径为4in或2in，高度为2in～4in。
3.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S2中的具体步骤为：观察岩心的层理及其他弱结构面的发育情况，筛选结构完整、层理均匀且无夹层、裂缝及其它结构面发育的岩心。
4.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S3中测试方位间隔为5°、10°或15°。
5.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S6中的声波测试井信息包括XMAC、SonicScan或DSI声波测试井信息。
6.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S6中的井周地层为岩心对应的无断层、无夹层和无其他非均质结构地层，且该井周地层的厚度不超过25m。
7.根据权利要求1所述的基于声波波速各向异性的岩心室内定向方法，其特征在于，所述步骤S7中岩心的地理方位的确定方法为：令岩心的最大波速相对方位减去波速最大方位，得到岩心的地理方位。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911026888.5A CN110794038B (zh) | 2019-10-26 | 2019-10-26 | 基于声波波速各向异性的岩心室内定向方法 |
| ZA2020/06605A ZA202006605B (en) | 2019-10-26 | 2020-10-23 | Laboratory core orientation method based on acoustic velocity anisotropy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911026888.5A CN110794038B (zh) | 2019-10-26 | 2019-10-26 | 基于声波波速各向异性的岩心室内定向方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110794038A CN110794038A (zh) | 2020-02-14 |
| CN110794038B true CN110794038B (zh) | 2021-03-19 |
Family
ID=69441511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911026888.5A Active CN110794038B (zh) | 2019-10-26 | 2019-10-26 | 基于声波波速各向异性的岩心室内定向方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN110794038B (zh) |
| ZA (1) | ZA202006605B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113533519B (zh) * | 2021-07-15 | 2022-09-20 | 厦门大学 | 一种非接触无损评估材料各向异性的方法和装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1321815A2 (ru) * | 1985-07-14 | 1987-07-07 | Кольский Филиал Им.С.М.Кирова Ан Ссср | Способ определени напр женного состо ни массива горных пород |
| CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
| CN103852785A (zh) * | 2012-11-28 | 2014-06-11 | 中国石油集团长城钻探工程有限公司 | 地层各向异性的评价方法 |
| CN104823048A (zh) * | 2013-05-06 | 2015-08-05 | 雪佛龙美国公司 | 用于从无定向的传统岩心确定储藏地质体的定向的系统和方法 |
| CN105136362A (zh) * | 2015-09-25 | 2015-12-09 | 中国石油大学(华东) | 一种基于岩石波速各向异性确定地应力方向的测量装置及方法 |
| CN109406019A (zh) * | 2018-11-16 | 2019-03-01 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种重现波速的地应力测试方法及装置 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105095631B (zh) * | 2014-05-21 | 2018-08-07 | 中国石油化工股份有限公司 | 一种页岩各向异性岩石物理建模方法 |
| CN104515738A (zh) * | 2014-12-09 | 2015-04-15 | 中国石油大学(北京) | 岩石各向异性检测系统和检测方法 |
-
2019
- 2019-10-26 CN CN201911026888.5A patent/CN110794038B/zh active Active
-
2020
- 2020-10-23 ZA ZA2020/06605A patent/ZA202006605B/en unknown
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1321815A2 (ru) * | 1985-07-14 | 1987-07-07 | Кольский Филиал Им.С.М.Кирова Ан Ссср | Способ определени напр женного состо ни массива горных пород |
| CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
| CN103852785A (zh) * | 2012-11-28 | 2014-06-11 | 中国石油集团长城钻探工程有限公司 | 地层各向异性的评价方法 |
| CN104823048A (zh) * | 2013-05-06 | 2015-08-05 | 雪佛龙美国公司 | 用于从无定向的传统岩心确定储藏地质体的定向的系统和方法 |
| CN105136362A (zh) * | 2015-09-25 | 2015-12-09 | 中国石油大学(华东) | 一种基于岩石波速各向异性确定地应力方向的测量装置及方法 |
| CN109406019A (zh) * | 2018-11-16 | 2019-03-01 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种重现波速的地应力测试方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 川西致密储层现今地应力场特征及石油工程地质应用研究;周文;《中国优秀博士学位论文全文数据库 基础科学辑》;20061215(第12期);第27-28，31-33页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110794038A (zh) | 2020-02-14 |
| ZA202006605B (en) | 2021-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9075158B2 (en) | Using a drill bit as a seismic source for SET velocity analysis | |
| US8154951B2 (en) | Model-based relative bearing estimation of three-component receivers | |
| NO339700B1 (no) | Fremgangsmåte og anordning for nedihulls formasjonsavbildning ved hjelp av retningsbestemt akustisk bølge-måling | |
| CN112068206B (zh) | 一种基于井孔与弹性波相互作用理论的邻井探测方法 | |
| CN107784159B (zh) | 一种储层电阻率各向异性系数的确定方法 | |
| CN108957548B (zh) | 一种多波多分量联合观测地震页岩气富集区预测方法 | |
| CN102681013A (zh) | 碳酸盐岩储层空间模型的建立方法及装置 | |
| CN107065014A (zh) | 一种基于微动探测结合地质钻探的地铁区间孤石群探测方法 | |
| CN108680965B (zh) | 一种适用于戈壁荒漠浅覆盖区快速找矿方法 | |
| CN108693562B (zh) | 一种用电成像测井测定古流向的方法 | |
| CN105804731B (zh) | 一种岩石地应力检测方法和系统 | |
| CN110794038B (zh) | 基于声波波速各向异性的岩心室内定向方法 | |
| Paillet et al. | Effects of lithology on televiewer-log quality and fracture interpretation | |
| US9081112B1 (en) | Borehole seismic system | |
| CN114442191B (zh) | 基于岩芯和井眼结构恢复确定井下现今地应力方向的方法 | |
| CN109657258B (zh) | 地层岩石力学特征剖面的处理方法和装置 | |
| CN111768104B (zh) | 单斜砂泥岩互层型坝基岩性变化程度定量评价方法 | |
| Perdana et al. | Characterizing subsurfac e structures and determining in-situ stress orientation using fullbore formatio n MicroImager (FMI) and sonic scanner: a case study from FORGE Well 21-31, fallon, Nevada, USA | |
| Hill et al. | Techniques for determining subsurface stress direction and assessing hydraulic fracture azimuth | |
| Croon et al. | Formation evaluation case study: Glyde unconventional Middle Proterozoic play in the McArthur Basin, northern Australia | |
| CN112147698B (zh) | 裂缝发育带识别与特征确定的方法及系统 | |
| CN117250658B (zh) | 建立研究区的地震数据集的方法 | |
| Das et al. | Integrated analysis of acoustic images and shear wave anisotropy used for geomechanical modelling in a tight sand oil reservoir | |
| CN115079268A (zh) | 一种基于椭圆短轴拟合的宽方位地震资料储层预测方法 | |
| Wilczynski et al. | Drill‐bit position monitoring using seismic‐while‐drilling data; numerical and field examples from Sweden |
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 |