CN110716235B - 一种砂泥岩测井孔隙结构反演方法 - Google Patents
一种砂泥岩测井孔隙结构反演方法 Download PDFInfo
- Publication number
- CN110716235B CN110716235B CN201810774378.5A CN201810774378A CN110716235B CN 110716235 B CN110716235 B CN 110716235B CN 201810774378 A CN201810774378 A CN 201810774378A CN 110716235 B CN110716235 B CN 110716235B
- Authority
- CN
- China
- Prior art keywords
- rock
- sand shale
- shale
- sand
- pore structure
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种砂泥岩测井孔隙结构反演方法,其特征在于,包括确定砂泥岩的岩石物理模型;根据砂泥岩孔隙结构特征对孔隙进行分类,结合测井资料,确定砂泥岩的孔隙结构。本发明能够充分结合砂泥岩实际地质特征,采用适用于砂泥岩岩性的Mori‑Tanaka等效介质模型,充分考虑各孔隙之间的相互影响,能够简捷、有效地反演出砂泥岩的孔隙结构,结果精度高,而且能够进行横波预测,满足常规的测井资料分析,提高应用范围,为测井岩石物理分析提供帮助。
Description
技术领域
本发明涉及一种砂泥岩测井孔隙结构反演方法,属于岩石物理研究领域。
背景技术
地震岩石物理方法在非常规油气藏勘探开发中起着重要的技术支持作用,能够有效的对储层物性、复杂性孔隙结构和含流体性进行描述。物性、孔隙结构、流体这三方面是储层预测和描述过程中重要的三个因素,孔隙度和孔隙结构对岩石弹性性质的影响是岩石勘探中重要方面之一。由于岩石物理分析方法是建立孔隙度、孔隙结构与弹性性质关联的有效途径,因此基于岩石物理方法进行孔隙反演是可靠的,也越来越受到重视。
地震岩石物理研究介质通常为双相介质,即岩石骨架和流体。在储层预测中孔隙度和流体是两个重要参数,孔隙结构影响着流体特性且制约着储层预测和开发的可行性。为了准确模拟和分析岩石中的孔隙结构,现今岩石物理学家通过将孔隙结构假想为理想椭球体,将孔隙(椭球体)的纵横比设定为一个固定的参数应用到岩石物理模型中。进一步地,Xu-white模型将岩石中的孔隙假设为两种不同纵横比来描述不同的孔隙类型,达到模拟复杂孔隙度结构的目的。Yan(2002)等利用神经网络算法分别反演了Xu-White模型的两种孔隙结构。Xu等基于微分等效理论对Xu-white模型进行修改,使其适用于碳酸盐岩储层并反演了储层的孔隙类型。Guo等进一步对Xu-Payne提出的碳酸盐岩模型进行改进将其应用于页岩储层,反演得到页岩储层的孔隙纵横比。但是,目前已有的孔隙结构反演方法考虑的孔隙结构较为简单,而且方法的局限性较大,不同岩性和地区采用的模型方法来反演孔隙结构也不同。由于砂泥岩微观结构较为复杂,矿物组分复杂,因此,已有的孔隙结构反演方法很难适用,需要构建合适的孔隙结构反演方法,为砂泥岩储层预测和开发提供技术支撑。
发明内容
针对上述问题,本发明的目的是提供一种砂泥岩测井孔隙结构反演方法,其能够简捷有效地确定砂泥岩的孔隙结构,结果精度高,而且能够进行横波预测,提高应用范围。
为实现上述目的,本发明采取以下技术方案:一种砂泥岩测井孔隙结构反演方法,包括确定砂泥岩的岩石物理模型;根据砂泥岩孔隙结构特征对孔隙进行分类,结合测井资料,确定砂泥岩的孔隙结构。
在一个具体的实施例中,根据砂泥岩的孔隙结构,结合砂泥岩的岩石物理模型确定砂泥岩岩石的横波速度。
在一个具体的实施例中,砂泥岩的岩石物理模型包括利用Voigt-Reuss-Hill平均模型,获取砂泥岩岩石矿物组分的等效弹性属性,利用Mori-Tanaka等效介质模型,获取砂泥岩包含各孔隙类型干岩石骨架的等效弹性属性,利用Gassmann方程,获取砂泥岩饱和流体时岩石的弹性属性。
在一个具体的实施例中,根据砂泥岩孔隙结构特征将孔隙分为三种类型,孔隙纵横比为0.01的微裂缝,孔隙纵横比为0.1的粒间孔,孔隙纵横比为0.9的硬孔隙。
在一个具体的实施例中,利用Voigt-Reuss-Hill平均模型确定砂泥岩岩石矿物组分的等效弹性属性所依据的公式:
式中,Mm为砂泥岩岩石矿物组分的等效弹性模量,Mi为砂泥岩岩石矿物组分第i个组成成分的模量;fi为砂泥岩岩石矿物组分第i个组成成分的体积分量;Mv为采用Voigt上限方法(公式(2))计算获得的砂泥岩的岩石模量;MR为采用Reuss下限方法(公式(3))计算获得的砂泥岩的岩石模量。
在一个具体的实施例中,利用Mori-Tanaka等效介质模型确定包含有各孔隙类型的砂泥岩干岩石骨架等效弹性模量所依据的公式:
VS+VR+VC=Φt (6)
式中,K0为砂泥岩岩石矿物组分的有效体积模量,μ0为砂泥岩岩石矿物组分的有效剪切模量,φt为孔隙度,Kd为砂泥岩干岩石骨架的弹性体积模量、μd为砂泥岩干岩石骨架的剪切模量,VS为硬孔隙所占孔隙度,VR为粒间孔所占孔隙度,VC为微孔隙所占孔隙度,P,Q为具有不同纵横比孔隙的几何因子。
在一个具体的实施例中,利用Gassmann方程确定砂泥岩饱和流体时岩石的弹性模量所依据的公式:
μsat=μdry (8)
式中,Kdry为砂泥岩干岩石骨架的有效体积模量,Ksat为砂泥岩饱和岩石的有效体积模量,Kf为砂泥岩饱和孔隙流体的有效体积模量。
在一个具体的实施例中,根据砂泥岩的岩石物理模型,结合各类孔隙结构所占比例,计算砂泥岩岩石的纵波速度,并进行砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差值判断,确定砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差范围是否在误差范围内,否则调整各孔隙结构所占比例,继续计算砂泥岩岩石的纵波速度,直至满足误差范围,确定砂泥岩岩石的孔隙结构。
在一个具体的实施例中,确定砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差小于30%。
在一个具体的实施例中,砂泥岩岩石矿物组分包括石英和泥质。
本发明由于采取以上技术方案,其具有以下优点:本发明能够充分结合砂泥岩实际地质特征,采用适用于砂泥岩岩性的Mori-Tanaka等效介质模型,充分考虑各孔隙之间的相互影响,能够简捷、有效地反演出砂泥岩的孔隙结构,结果精度高,而且能够进行横波预测,满足常规的测井资料分析,提高应用范围,为测井岩石物理分析提供帮助。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分的从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚的说明本发明实施例中的技术方案,下面将对实施例描述中所需要的附图做简单的介绍:
图1是本发明的一个具体实施例的结构示意图;
图2是本发明的一个具体实施例孔隙结构反演确定的纵波速度和测井纵波速度的结构示意图;
图3是本发明的一个具体实施例孔隙结构反演确定的纵波速度与测井纵波速度误差的结构示意图;
图4是本发明的一个具体实施例确定的孔隙结构硬孔隙孔隙度的结构示意图;
图5是本发明的一个具体实施例确定的孔隙结构中粒间孔隙度的结构示意图;
图6是本发明的一个具体实施例确定的孔隙结构微裂缝孔隙度的结构示意图;
图7是本发明的一个具体实施例粘土含量估算值与实际测井粘土含量值的结构示意图;
图8是本发明的一个具体实施例横波速度预测值的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
如图1所示,本发明提出的砂泥岩测井孔隙结构反演方法,包括
1)确定砂泥岩的岩石物理模型。
砂泥岩的岩石物理模型包括利用Voigt-Reuss-Hill平均模型,获取砂泥岩岩石矿物组分的等效弹性属性,利用Mori-Tanaka等效介质模型,获取砂泥岩包含各孔隙类型干岩石骨架的等效弹性属性,利用Gassmann方程,获取砂泥岩饱和流体时岩石的弹性属性。
在此基础上,确定的砂泥岩的岩石物理模型能够符合砂泥岩的实际岩性。
2)根据砂泥岩孔隙结构特征对孔隙进行分类,结合测井资料,进行孔隙结构反演,确定砂泥岩目标孔隙结构。
根据砂泥岩孔隙结构特征将孔隙分为三种类型,孔隙纵横比为0.01的微裂缝,孔隙纵横比为0.1的粒间孔,孔隙纵横比为0.9的硬孔隙。
测井资料包括砂泥岩的孔隙度测井数据,密度测井数据和岩石矿物组分中泥质含量测井数据,以及砂泥岩的测井纵波速度数据。
在一个优选的实施例中,砂泥岩岩石矿物组分包括石英和泥质。
首先,利用Voigt-Reuss-Hill平均模型确定砂泥岩岩石矿物组分的等效弹性属性所依据的公式:
式中,Mm为砂泥岩岩石矿物组分的等效弹性模量,Mi为砂泥岩岩石矿物组分第i个组成成分的模量;fi为砂泥岩岩石矿物组分第i个组成成分的体积分量;Mv为采用Voigt上限方法(公式(2))计算获得的砂泥岩的岩石模量;MR为采用Reuss下限方法(公式(3))计算获得的砂泥岩的岩石模量。
接着,利用Mori-Tanaka等效介质模型确定包含有各孔隙类型的砂泥岩干岩石骨架等效弹性模量所依据的公式:
VS+VR+VC=Φt (6)
式中,K0为砂泥岩岩石矿物组分的有效体积模量,μ0为砂泥岩岩石矿物组分的有效剪切模量,φt为孔隙度,Kd为砂泥岩干岩石骨架的弹性体积模量、μd为砂泥岩干岩石骨架的剪切模量,VS为硬孔隙所占孔隙度,VR为粒间孔所占孔隙度,VC为微孔隙所占孔隙度,P,Q为具有不同纵横比孔隙的几何因子。
然后,利用Gassmann方程确定砂泥岩饱和流体时岩石的弹性模量所依据的公式:
μsat=μdry (8)
式中,Kdry为砂泥岩干岩石骨架的有效体积模量,Ksat为砂泥岩饱和岩石的有效体积模量,Kf为砂泥岩饱和孔隙流体的有效体积模量。
在公式4)和公式5)中,计算P、Q的公式:
公式,Tiijj、Tijij分别为椭球体包含物的四阶弹性张量。
其中,确定P、Q所需要相关标量的公式:
Tiijj=3F1/F2 (11)
F6=1+A[1+f-R(f+θ)]+B(1-θ)×(3-4R) (18)
F9=A[(R-1)f-Rθ]+Bθ(3-4R) (21)
R=(1-2υsc)/2(1-υsc) (24)
式中,μj为砂泥岩组成岩石孔隙成分的剪切模量,Kj为砂泥岩组成岩石孔隙成分的体积模量,α为孔隙的纵横比,Ksc为砂泥岩岩石矿物组分的体积模量,μsc为砂泥岩岩石矿物组分的剪切模量,υsc为砂泥岩岩石矿物组分的泊松比,F1~F9、f、θ、R、A、B分别为中间变量。
根据确定的砂泥岩的岩石物理模型,结合各类孔隙结构所占比例,计算砂泥岩岩石的纵波速度。
假设砂泥岩的初始孔隙纵横比为0.1的孔隙含量为100%,纵横比为0.01孔隙含量为0%,纵横比为0.9孔隙含量为0%,调整三种孔隙所占百分含量,通过砂泥岩的岩石物理模型计算砂泥岩岩石的纵波速度值。
进行砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差判断。确定砂泥岩岩石的纵波速度与测井纵波速度的相对误差是否在误差范围内(相对误差<30%),如果满足误差范围,则确定当前的孔隙结构配比。当误差范围不满足时,重新调整砂泥岩的孔隙结构比例,继续计算砂泥岩岩石的纵波速度,直至满足误差范围。最后,确定砂泥岩岩石的孔隙结构。
3)根据砂泥岩确定的目标孔隙结构,结合砂泥岩的岩石物理模型确定砂泥岩岩石的横波速度。
根据砂泥岩反演确定的孔隙结构,结合砂泥岩的岩石物理模型计算出砂泥岩岩石的横波速度,能够达到预测横波的目的,从而满足常规的测井资料分析。
下面列举一具体实施例
本实施结合实际的砂泥岩测井资料,进行孔隙结构的反演。图2、图3所示,可以看出估算的纵波速度与原始测井纵波速度值吻合较好,误差整体为0.5,预测精度高。如图4~6所示,从反演结果的分布上看,微裂隙含量较少,硬孔隙含量最多,这符合砂泥岩特征,反演结果合理。如图7、图8所示,得到的横波速度可用于常规的测井资料分析,结合孔隙结构反演结果弥补了测井横波测试的难题,为测井岩石物理分析提供了帮助。
虽然本发明所公开的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所公开的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (7)
1.一种砂泥岩测井孔隙结构反演方法,其特征在于,包括
确定砂泥岩的岩石物理模型;砂泥岩的岩石物理模型包括利用Voigt-Reuss-Hill平均模型,获取砂泥岩岩石矿物组分的等效弹性属性,利用Mori-Tanaka等效介质模型,获取砂泥岩包含各孔隙类型干岩石骨架的等效弹性属性,利用Gassmann方程,获取砂泥岩饱和流体时岩石的弹性属性;
根据砂泥岩孔隙结构特征对孔隙进行分类,结合测井资料,确定砂泥岩的孔隙结构;其中,根据砂泥岩孔隙结构特征对孔隙进行分类包括根据砂泥岩孔隙结构特征将孔隙分为三种类型,孔隙纵横比为0.01的微裂缝,孔隙纵横比为0.1的粒间孔,孔隙纵横比为0.9的硬孔隙;
其中,利用Mori-Tanaka等效介质模型确定包含有各孔隙类型的砂泥岩干岩石骨架等效弹性模量所依据的公式:
VS+VR+VC=Φt (6)
式中,K0为砂泥岩岩石矿物组分的有效体积模量,μ0为砂泥岩岩石矿物组分的有效剪切模量,φt为孔隙度,Kd为砂泥岩干岩石骨架的弹性体积模量、μd为砂泥岩干岩石骨架的剪切模量,VS为硬孔隙所占孔隙度,VR为粒间孔所占孔隙度,VC为微孔隙所占孔隙度,P,Q为具有不同纵横比孔隙的几何因子。
2.根据权利要求1所述的一种砂泥岩测井孔隙结构反演方法,其特征在于,根据砂泥岩的孔隙结构,结合砂泥岩的岩石物理模型确定砂泥岩岩石的横波速度。
5.根据权利要求4所述的一种砂泥岩测井孔隙结构反演方法,其特征在于:根据砂泥岩的岩石物理模型,结合各类孔隙结构所占比例,计算砂泥岩岩石的纵波速度,并进行砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差值判断,确定砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差范围是否在误差范围内,否则调整各孔隙结构所占比例,继续计算砂泥岩岩石的纵波速度,直至满足误差范围,确定砂泥岩岩石的孔隙结构。
6.根据权利要求5所述的一种砂泥岩测井孔隙结构反演方法,其特征在于:确定砂泥岩岩石的纵波速度与砂泥岩测井纵波速度的相对误差小于30%。
7.根据权利要求1所述的一种砂泥岩测井孔隙结构反演方法,其特征在于,砂泥岩岩石矿物组分包括石英和泥质。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810774378.5A CN110716235B (zh) | 2018-07-13 | 2018-07-13 | 一种砂泥岩测井孔隙结构反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810774378.5A CN110716235B (zh) | 2018-07-13 | 2018-07-13 | 一种砂泥岩测井孔隙结构反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110716235A CN110716235A (zh) | 2020-01-21 |
CN110716235B true CN110716235B (zh) | 2022-03-08 |
Family
ID=69208699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810774378.5A Active CN110716235B (zh) | 2018-07-13 | 2018-07-13 | 一种砂泥岩测井孔隙结构反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110716235B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111595752B (zh) * | 2020-06-22 | 2021-07-06 | 中国科学技术大学 | 一种岩石有效孔隙度确定方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2490140A1 (de) * | 2010-08-19 | 2012-08-22 | PART Engineering GmbH | Computerimplementierters Verfahren zur Bestimmung des elasto-plastischen, mechanischen Verhaltens von inhomogenen, anisotropen Werkstoffen |
CN104007482A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 一种基于各向异性有效场的泥页岩岩石物理模型方法 |
CN104570129A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气股份有限公司 | 碳酸盐岩岩石储层横波速度预测方法 |
CN106443780A (zh) * | 2016-08-31 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司 | 一种页岩气地层的横波速度估算方法 |
CN108072902A (zh) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | 一种碳酸盐岩岩石物理模型的建模方法 |
-
2018
- 2018-07-13 CN CN201810774378.5A patent/CN110716235B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2490140A1 (de) * | 2010-08-19 | 2012-08-22 | PART Engineering GmbH | Computerimplementierters Verfahren zur Bestimmung des elasto-plastischen, mechanischen Verhaltens von inhomogenen, anisotropen Werkstoffen |
CN104007482A (zh) * | 2014-04-16 | 2014-08-27 | 孙赞东 | 一种基于各向异性有效场的泥页岩岩石物理模型方法 |
CN104570129A (zh) * | 2014-12-19 | 2015-04-29 | 中国石油天然气股份有限公司 | 碳酸盐岩岩石储层横波速度预测方法 |
CN106443780A (zh) * | 2016-08-31 | 2017-02-22 | 中国石油集团川庆钻探工程有限公司 | 一种页岩气地层的横波速度估算方法 |
CN108072902A (zh) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | 一种碳酸盐岩岩石物理模型的建模方法 |
Non-Patent Citations (3)
Title |
---|
基于储层砂岩微观孔隙结构特征的弹性波频散响应分析;邓继新 等;《地球物理学报》;20150930;第58卷(第9期);第3389-3400页 * |
基于岩石物理模型的页岩孔隙结构反演及横波速度预测;逄硕 等;《吉林大学学报(地球科学版)》;20170331;第47卷(第2期);第606-615页 * |
部分饱和砂岩储层岩石物理和地震响应研究;司文朋;《中国博士学位论文全文数据库•基础科学辑》;20180215(第2期);第12-13,21页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110716235A (zh) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
Yu et al. | Porosity estimation in kerogen-bearing shale gas reservoirs | |
CN103235338B (zh) | 一种反演岩石裂隙参数的方法 | |
CN105182424A (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN105426612A (zh) | 一种地层组分最优化确定方法及装置 | |
CN112255688B (zh) | 一种基于岩石物理理论的三维地震反演地层压力的方法 | |
CN105653815A (zh) | 一种基于岩石物理模型理论的定量解释油藏流体分布方法 | |
CN111399044B (zh) | 一种储层渗透率预测方法、装置及存储介质 | |
WO2006062612A2 (en) | Integrated anisotropic rock physics model | |
CN104570129A (zh) | 碳酸盐岩岩石储层横波速度预测方法 | |
CN103643949A (zh) | 一种储层含油气的定量预测方法及装置 | |
CN112505772B (zh) | 一种利用孔隙、裂隙介质弹性波理论反演岩石孔隙分布特征的方法 | |
CN102913240B (zh) | 一种储层流体识别方法 | |
CN110703330B (zh) | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 | |
CN112965103B (zh) | 一种多重孔隙储层叠前地震概率化多道反演方法 | |
CN106154343A (zh) | 计算致密油储层的含油饱和度的方法 | |
CN109839676B (zh) | 一种基于记忆模拟退火的基质模量估算方法及电子设备 | |
CN110716235B (zh) | 一种砂泥岩测井孔隙结构反演方法 | |
CN107132574B (zh) | 一种海相地层碎屑岩岩性组合的预测方法 | |
CN106568918B (zh) | 页岩有机碳含量toc预测方法 | |
CN114063162B (zh) | 一种基于小样本机器学习SVM Vp/Vs预测方法 | |
Khoshdel et al. | Permeability estimation using rock physics modeling and seismic inversion in a carbonate reservoir | |
CN110007348A (zh) | 一种灰质背景浊积岩储层的岩石物理建模方法 | |
CN115685335A (zh) | 一种纵横波速度预测方法、装置、计算机设备及存储介质 | |
CN114675328A (zh) | 基于变化孔隙纵横比Xu-white模型的富有机质页岩横波预测方法 |
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 |