CN106019400B - 一种获取塑性指数的方法 - Google Patents
一种获取塑性指数的方法 Download PDFInfo
- Publication number
- CN106019400B CN106019400B CN201510117142.0A CN201510117142A CN106019400B CN 106019400 B CN106019400 B CN 106019400B CN 201510117142 A CN201510117142 A CN 201510117142A CN 106019400 B CN106019400 B CN 106019400B
- Authority
- CN
- China
- Prior art keywords
- cap rock
- analyzed
- confining pressure
- rock
- plasticity
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种获取塑性指数的方法,该方法包括:获取待分析盖层所处地层区域的地层参数,根据地层参数确定待分析盖层的盖层深度;根据待分析盖层的盖层深度,确定待分析盖层的围压;根据待分析盖层的围压,基于预设塑性围压模型,确定待分析盖层的塑性指数。本发明所提供的方法能够获得不同地质历史时期、不同埋深条件下盖层的塑性指数。
Description
技术领域
本发明涉及油气勘探技术领域,具体地说,涉及一种获取塑性指数的方法。
背景技术
石油天然气的盖层封闭性主要取决于岩石物理性质,以往的盖层封闭性定量评价指标主要以物性评价为主,如孔隙度、渗透率和突破压力。但是,越来越多的证据显示,盖层的裂缝对油气的保存起到了至关重要的作用,而盖层裂缝发育的类型、规模与密度与盖层的脆塑性密切关系。
影响盖层脆塑性的因素主要包括:盖层的自身矿物组成、排列方式及其成岩演化程度、盖层的地质埋深、应力环境等。目前可以用三轴岩石力学实验来测试得到单个样品的塑性指数。但是这种方法仅能够用来评价现今盖层的塑性,不能确定在某一地质历史时期盖层塑性指数的变化。
基于上述情况,亟需一种能够获取地质历史时期的盖层塑性指数的方法。
发明内容
本发明所要解决的技术问题是获取特定地质历史时期的盖层塑性指数。为解决上述问题,本发明的实施例首先提供了一种获取泥岩盖层塑性指数的方法,这种方法既可以定量评价现今埋藏条件下泥岩盖层的脆塑性,也可以定量评价不同地质历史时期泥岩盖层的脆塑性。所述方法包括:
获取待分析盖层所处地层区域的地层参数,根据所述地层参数确定所述待分析盖层的盖层深度;
根据所述待分析盖层的盖层深度,确定所述待分析盖层的围压;
根据所述待分析盖层的围压,基于预设塑性围压模型,确定所述待分析盖层的塑性指数。
根据本发明的一个实施例,确定所述待分析盖层的盖层深度的步骤包括:
根据所述地层参数,恢复得到所述待分析盖层所处地层区域的地层埋藏史;
根据所述地层埋藏史,确定所述待分析盖层的盖层深度。
根据本发明的一个实施例,确定所述待分析盖层的围压的步骤包括:
获取所述待分析盖层所处区域以及相邻区域的测井曲线,根据所述测井曲线确定所述待分析盖层的上覆地层的岩石密度;
获取待分析盖层的上覆地层的水密度,结合所述上覆地层的岩石密度和待分析盖层的盖层深度,确定出待分析盖层的围压。
根据本发明的一个实施例,根据如下表达式确定所述待分析盖层的围压:
σc=0.010133×(ρd-ρw)×H
其中,σc和H分别表示待分析盖层的围压和盖层深度,ρd和ρw分别表示待分析盖层的上覆地层的岩石密度和水密度。
根据本发明的一个实施例,构建所述预设塑性围压模型的步骤包括:
根据所述待分析盖层所述地层区域的盖层分布特性,获取实验盖层岩心,并确定出所述实验盖层岩心的峰值抗压强度和残余抗压强度;
根据所述峰值抗压强度和残余抗压强度,分别确定围压与峰值抗压强度和残余抗压强度的函数关系;
根据围压与峰值抗压强度和残余抗压强度的函数关系,构建得到所述预设塑性围压模型。
根据本发明的一个实施例,对所述实验盖层岩心施加不同的围压来进行盖层破裂测试,根据测试结果分析得到所述实验盖层岩心的峰值抗压强度和残余抗压强度;
根据本发明的一个实施例,对所述实验盖层岩心施加至少5组不同的围压,来进行盖层破裂测试。
根据本发明的一个实施例,利用多项式拟合分别对围压与峰值抗压强度和残余抗压强度进行非线性回归,分别得到围压与峰值抗压强度和残余抗压强度的函数关系。
根据本发明的一个实施例,所述预设塑性围压模型为:
其中,fd表示塑性指数,σp(σc)表示峰值抗压强度σp关于围压σc的函数,σr(σc)表示残余抗压强度σr关于围压σc的函数。
根据本发明的一个实施例,所述地层参数包括以下所列项中的至少一项:
地层厚度、不整合剥蚀量和古地温。
现有的盖层塑性指数获取方法都是基于现今盖层样品的静态分析,仅仅能够获取到现今盖层的塑性指数。而本发明所提供的盖层塑性指数获取方法则是在盖层岩石力学测试的基础上构建出了盖层深度与围压的关系以及围压与塑性指数的关系,这样也就可以获得不同地质历史时期、不同埋深条件下盖层的塑性指数。利用获取到的盖层塑性指数,可以对盖层产生破裂的难以程度进行分析,这为分析油气成藏与破坏历史提供了参考与支撑,也为寻找有利勘探区带提供了依据。
同时,本发明所提供的盖层塑性指数获取方法不但能够获取到某个特定地质历史时期的盖层塑性指数,而且还能够与盆地的构造演化和地层的升降历史相结合来动态分析某个地质历史阶段内盖层塑性指数的变化特征。该方法具有操作方便和实用价值高等特点,其可以应用到含油气盆地的盖层分析中,以分析盖层封闭性动态演化特点。
本发明所提供的盖层塑性指数获取方法弥补了传统方法仅从物性封闭分析盖层的缺陷,可以定量评价盖层的韧脆性,预测盖层产生裂缝的风险性,从而更加全面地分析盖层封闭性。利用该方法确定出的泥岩盖层塑性形变深度,对页岩气储层的压裂也有指导意义。
本发明的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要的附图做简单的介绍:
图1是根据本发明一个实施例的获取盖层塑性指数的流程图;
图2是根据本发明一个实施例的确定待分析盖层的围压的流程图;
图3是根据本发明一个实施例的构建预设塑性围压模型的流程图;
图4是根据本发明一个实施例的三轴应力应变曲线图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式,借此对本发明如何应用技术手段来解决技术问题,并达成技术效果的实现过程能充分理解并据以实施。需要说明的是,只要不构成冲突,本发明中的各个实施例以及各实施例中的各个特征可以相互结合,所形成的技术方案均在本发明的保护范围之内。
同时,在以下说明中,出于解释的目的而阐述了许多具体细节,以提供对本发明实施例的彻底理解。然而,对本领域的技术人员来说显而易见的是,本发明可以不用这里的具体细节或者所描述的特定方式来实施。
目前关于油气盖层的塑性分析都是基于现今样本的静态评价,本发明的目的在于解决上述现有技术中存在的缺陷,以提供一种能够获取不同地质历史时期盖层的塑性指数的方法。
图1示出了本实施例所提供的盖层塑性指数获取方法的流程图。
如图1所示,本实施例所所提供的方法在步骤S101中获取待分析盖层所处地层区域的地层参数,并根据获取到的地层参数恢复待分析盖层所处地层区域的地层埋藏史。
本实施例中,在步骤S101中所获取到的地层参数包括地层厚度、不整合剥蚀量和古地温等参数。基于这些参数,利用盆地模拟器恢复得到待分析盖层所处地层区域的地层埋藏史。
当然,在本发明的其他实施例中,在步骤S101中所获取到的地层参数可以包括一项所列项中的任一项或几项,也可以包括未列出的合理参数,本发明不限于此。
随后在步骤S102中,根据步骤S101确定出的地层埋藏史图,读取出待分析盖层所处地质历史时期的盖层深度。需要说明的是,在本发明的其他实施例中,还可以利用所获取的地层参数,采用其他合理的方式来确定待分析盖层的盖层深度,本发明不限于此。
得到待分析盖层的盖层深度后,在步骤S103中根据待分析盖层的盖层深度,确定出待分析盖层的围压。图2示出了本实施例中确定待分析盖层的围压的流程图。
如图2所示,首先在步骤S201中从待分析盖层所处区域以及相邻区域中选取能够表征这些区域地质特性的典型钻井,并利用这些钻井得到声波密度测井曲线。随后在步骤S202中,利用步骤S201中所得到的声波密度测井曲线,确定出待分析盖层的上覆地层的岩石密度。在步骤S203中,获取待分析盖层的上覆地层的水密度,并在步骤S204中根据水密度、岩石密度和盖层深度,确定出待分析盖层的围压。
本实施例中,根据如下表达式计算待分析盖层的围压:
σc=0.010133×(ρd-ρw)×H (1)
其中,σc和H分别表示待分析盖层的围压和盖层深度,ρd和ρw分别表示待分析盖层的上覆地层的岩石密度和水密度。
再次如图1所示,得到待分析盖层的围压后,在步骤S104中根据待分析盖层的围压,基于预设塑性围压模型,确定出待分析盖层的塑性指数。预设塑性围压模型能够表征待分析盖层所述地层区域中围压与塑性指数的函数关系,它是通过岩石力学实验构建得到的。
图3示出了本实施例中构建预设塑性围压模型的流程图。
如图3所示,本实施例在构建预设塑性围压模型时,首先在步骤S301中根据待分析盖层所述地层区域的盖层分布特性,获取实验盖层岩心,并确定出实验盖层岩心的峰值抗压强度和残余抗压强度。具体地,本实施例利用三轴应力实验,对实验盖层岩心施加5组不同的围压来进行盖层破裂测试,分析得到实验盖层岩心的峰值抗压强度和残余抗压强度。
图4示出了本实施例中一组典型的三轴应力应变曲线图,该图示出了在以确定围压条件下,应力与应变的关系曲线。从图4中可以看出,在一组三轴应力实验中,能够分别确定出一组峰值抗压强度数据和残余抗压强度数据。这样,通过改变围压来进行多组实验,便可以获得多个围压条件下的峰值抗压强度和残余抗压强度。
当然,在本发明的其他实施例中,对实验盖层岩心所进行的盖层破裂测试的次数还可以为其他合理值,例如大于5组等,本发明不限于此。同时,还可以利用实验盖层岩心,根据其他合理的方式来确定实验盖层岩心的峰值抗压强度和残余抗压强度,本发明同样不限于此。
得到各个围压所对应的峰值抗压强度和残余抗压强度后,在步骤S302中,利用多项式拟合分别对围压与峰值抗压强度以及围压与残余抗压强度进行非线性回归,分别确定出围压为峰值抗压强度的函数关系以及围压与残余抗压强度的函数关系。
本实施例中,围压与峰值抗压强度的函数关系可以表示为:
围压与残余抗压强度的函数关系可以表示为:
其中,σp和σr分别表示峰值抗压强度和残余抗压强度,σc表示围压,a1、a2、b1、b2、c1和c2表示常数。
需要说明的是,在本发明的其他实施例中,还可以利用所得围压数据以及对应的峰值抗压强度数据和残余抗压强度数据,采用其他合理的方法来确定围压与峰值抗压强度的函数关系以及围压与残余抗压强度的函数关系。
在步骤S303中,根据步骤S302中所得到的围压与峰值抗压强度的函数关系以及围压与残余抗压强度的函数关系,结合塑性指数的物理含义,构建得到塑性指数与围压的函数关系,即预设塑性围压模型。
本实施例中,塑性指数与围压的函数关系可以采用如下表达式表示:
其中,fd表示塑性指数,σp(σc)表示峰值抗压强度σp关于围压σc的函数,σr(σc)表示残余抗压强度σr关于围压σc的函数。
再次如图1所示,本实施例所提供的方法在步骤S103中确定出了待分析盖层的围压,而预先构建得到的塑性围压模型能够表示出盖层的塑性指数与围压的函数关系,因此将待分析盖层的围压代入预设塑性围压模型即可确定出待分析盖层的塑性指数。为了进一步说明本实施例所提供的盖层塑性指数获取方法的可靠性和实用性,本实施例还利用该方法分别构建出了四川盆地东部志留系盖层的塑性指数与围压的函数关系以及围压与盖层深度的函数关系,并最终得到了盖层的塑性指数与盖层深度的关系,其结果如下所示:
本实施例分别获取到了四川盆地东部现今不同埋深条件下和X井200Ma以来志留系盖层的塑性指数,其结果如表1所示。
表1
序号 | 地质年代(Ma) | 深度(m) | 围压(MPa) | fd |
1 | 现今 | 0 | 0 | 0 |
2 | 现今 | 1000 | 17 | 0.24 |
3 | 现今 | 2000 | 34 | 0.37 |
4 | 现今 | 3000 | 51 | 0.53 |
5 | 现今 | 4000 | 68 | 0.79 |
6 | 现今 | 5000 | 85 | 1 |
7 | 200 | 3950 | 67.15 | 0.77 |
8 | 150 | 4490 | 76.33 | 1 |
9 | 100 | 5760 | 97.92 | 1 |
10 | 50 | 4400 | 74.8 | 0.96 |
从表1所示的数据可以看出,利用本实施例所提供的方法获取到的盖层塑性指数与油气勘探实践得到的数据相吻合。与传统的方法相比,本实施例所提供的方法由于采用具体的岩石力学参数做支撑,因此具有更高的准确性和可靠性。
从上述描述中可以看出,现有的盖层塑性指数获取方法都是基于现今盖层样品的静态分析,仅仅能够获取到现今盖层的塑性指数。而本实施例所提供的盖层塑性指数获取方法则是在盖层岩石力学测试的基础上构建出了盖层深度与围压的关系以及围压与塑性指数的关系,这样也就可以获得不同地质历史时期、不同埋深条件下盖层的塑性指数。利用获取到的盖层塑性指数,可以对盖层产生破裂的难以程度进行分析,这为分析油气成藏与破坏历史提供了参考与支撑,也为寻找有利勘探区带提供了依据。
同时,本实施例所提供的盖层塑性指数获取方法不但能够获取到某个特定地质历史时期的盖层塑性指数,而且还能够与盆地的构造演化和地层的升降历史相结合来动态分析某个地质历史阶段内盖层塑性指数的变化特征。本实施例所提供的方法具有操作方便和实用价值高等特点,其可以应用到含油气盆地的盖层分析中,以分析盖层封闭性动态演化特点。
应该理解的是,本发明所公开的实施例不限于这里所公开的特定结构、处理步骤或材料,而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是,在此使用的术语仅用于描述特定实施例的目的,而并不意味着限制。
说明书中提到的“一个实施例”或“实施例”意指结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中。因此,说明书通篇各个地方出现的短语“一个实施例”或“实施例”并不一定均指同一个实施例。
为了方便,在此使用的多个项目、结构单元、组成单元和/或材料可出现在共同列表中。然而,这些列表应解释为该列表中的每个元素分别识别为单独唯一的成员。因此,在没有反面说明的情况下,该列表中没有一个成员可仅基于它们出现在共同列表中便被解释为相同列表的任何其它成员的实际等同物。另外,在此还可以连同针对各元件的替代一起来参照本发明的各种实施例和示例。应当理解的是,这些实施例、示例和替代并不解释为彼此的等同物,而被认为是本发明的单独自主的代表。
虽然上述示例用于说明本发明在一个或多个应用中的原理,但对于本领域的技术人员来说,在不背离本发明的原理和思想的情况下,明显可以在形式上、用法及实施的细节上作各种修改而不用付出创造性劳动。因此,本发明由所附的权利要求书来限定。
Claims (8)
1.一种获取塑性指数的方法,其特征在于,所述方法包括:
获取待分析盖层所处地层区域的地层参数,根据所述地层参数确定所述待分析盖层的盖层深度;
根据所述待分析盖层的盖层深度,确定所述待分析盖层的围压;
根据所述待分析盖层的围压,基于预设塑性围压模型,确定所述待分析盖层的塑性指数;
其中,构建所述预设塑性围压模型的步骤包括:
根据所述待分析盖层所处地层区域的盖层分布特性,获取实验盖层岩心,并确定出所述实验盖层岩心的峰值抗压强度和残余抗压强度;
根据所述峰值抗压强度和残余抗压强度,分别确定围压与峰值抗压强度和残余抗压强度的函数关系;
根据围压与峰值抗压强度和残余抗压强度的函数关系,构建得到所述预设塑性围压模型;
其中,所述预设塑性围压模型为:
其中,fd表示塑性指数,σp(σc)表示峰值抗压强度σp关于围压σc的函数,σr(σc)表示残余抗压强度σr关于围压σc的函数。
2.如权利要求1所述的方法,其特征在于,确定所述待分析盖层的盖层深度的步骤包括:
根据所述地层参数,恢复得到所述待分析盖层所处地层区域的地层埋藏史;
根据所述地层埋藏史,确定所述待分析盖层的盖层深度。
3.如权利要求1或2所述的方法,其特征在于,确定所述待分析盖层的围压的步骤包括:
获取所述待分析盖层所处地层区域以及相邻区域的测井曲线,根据所述测井曲线确定所述待分析盖层的上覆地层的岩石密度;
获取待分析盖层的上覆地层的水密度,结合所述上覆地层的岩石密度和待分析盖层的盖层深度,确定出待分析盖层的围压。
4.如权利要求3所述的方法,其特征在于,根据如下表达式确定所述待分析盖层的围压:
σc=0.010133×(ρd-ρw)×H
其中,σc和H分别表示待分析盖层的围压和盖层深度,ρd和ρw分别表示待分析盖层的上覆地层的岩石密度和水密度。
5.如权利要求1所述的方法,其特征在于,对所述实验盖层岩心施加不同的围压来进行盖层破裂测试,根据测试结果分析得到所述实验盖层岩心的峰值抗压强度和残余抗压强度。
6.如权利要求5所述的方法,其特征在于,对所述实验盖层岩心施加至少5组不同的围压,来进行盖层破裂测试。
7.如权利要求5或6所述的方法,其特征在于,利用多项式拟合分别对围压与峰值抗压强度和残余抗压强度进行非线性回归,分别得到围压与峰值抗压强度和残余抗压强度的函数关系。
8.如权利要求1所述的方法,其特征在于,所述地层参数包括以下所列项中的至少一项:
地层厚度、不整合剥蚀量和古地温。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510117142.0A CN106019400B (zh) | 2015-03-17 | 2015-03-17 | 一种获取塑性指数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510117142.0A CN106019400B (zh) | 2015-03-17 | 2015-03-17 | 一种获取塑性指数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106019400A CN106019400A (zh) | 2016-10-12 |
CN106019400B true CN106019400B (zh) | 2018-11-20 |
Family
ID=57082344
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510117142.0A Active CN106019400B (zh) | 2015-03-17 | 2015-03-17 | 一种获取塑性指数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106019400B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108572391A (zh) * | 2017-03-08 | 2018-09-25 | 中国石油化工股份有限公司 | 一种反演盖层塑性的方法 |
CN110231407B (zh) * | 2018-03-06 | 2022-02-15 | 中国石油化工股份有限公司 | 一种判断碳酸盐岩盖层有效性的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8515721B2 (en) * | 2009-10-01 | 2013-08-20 | Schlumberger Technology Corporation | Method for integrated inversion determination of rock and fluid properties of earth formations |
CN103206206B (zh) * | 2012-01-16 | 2016-08-03 | 中国石油化工股份有限公司 | 用于计算超深井地层破裂压力的方法 |
CN102749660B (zh) * | 2012-06-26 | 2015-01-07 | 中国人民解放军第二炮兵工程设计研究所 | 高地应力地区近水平岩层岩爆的综合预报方法 |
FR3007533B1 (fr) * | 2013-06-20 | 2015-07-24 | Ifp Energies Now | Procede pour predire la quantite et la composition des fluides produits par des reactions minerales operant dans un bassin sedimentaire |
CN104252009B (zh) * | 2013-06-27 | 2017-01-25 | 中国石油化工股份有限公司 | 一种恢复石油与天然气盖层脆塑性演化史的方法 |
-
2015
- 2015-03-17 CN CN201510117142.0A patent/CN106019400B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106019400A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xu et al. | Characterization of rock mechanical properties using lab tests and numerical interpretation model of well logs | |
Luffel et al. | Travis Peak core permeability and porosity relationships at reservoir stress | |
Yielding | Using probabilistic shale smear modelling to relate SGR predictions of column height to fault-zone heterogeneity | |
An et al. | Modeling anisotropic permeability of coal and its effects on CO2 sequestration and enhanced coalbed methane recovery | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
Germay et al. | The continuous-scratch profile: a high-resolution strength log for geomechanical and petrophysical characterization of rocks | |
CN107290799B (zh) | 一种岩石可压性的确定方法 | |
da Silva Jr et al. | Petrophysical studies of north American carbonate rock samples and evaluation of pore-volume compressibility models | |
CN106019400B (zh) | 一种获取塑性指数的方法 | |
CN109577969B (zh) | 一种基于岩石压缩系数计算碳酸盐岩地层孔隙压力的方法 | |
Qian et al. | Extraction and identification of spectrum characteristics of coal and rock hydraulic fracturing and uniaxial compression signals | |
Ma et al. | Experimental and numerical study on fracture propagation near open-hole horizontal well under hydraulic pressure | |
Salman et al. | Permeability prediction in one of Iraqi carbonate reservoir using statistical, hydraulic flow units, and ANN methods | |
Shin et al. | On computation of strain‐dependent permeability of rocks and rock‐like porous media | |
US9989667B2 (en) | Pore size classification in subterranean formations based on nuclear magnetic resonance (NMR) relaxation distributions | |
Guo et al. | Numerical Simulation and Experimental Studies of Karst Caves Collapse Mechanism in Fractured‐Vuggy Reservoirs | |
CN109184676A (zh) | 一种页岩气藏有效改造体积评价方法 | |
Lopez et al. | Estimates of stress dependent properties in tight reservoirs: Their use with drill cuttings data | |
Liu et al. | A new method for estimating the free-to-adsorbed ratio in shale gas reservoirs using acoustic amplitude attenuation and porosity | |
CN111257134A (zh) | 剪切应力的连续深度处理方法 | |
Hasanov et al. | Graphical 3D Detailing of the Bottomhole Part of the Well with the Identification of Drainage Zones within of Productive Interval | |
CN108132342B (zh) | 岩石静态杨氏模量的确定方法 | |
Lobo et al. | Methodology for Petrophysical and Geomechanical Analysis of Shale Plays. Case Study: La Luna and Capacho Formations, Maracaibo Basin. | |
Kitamura et al. | Potential evaluation of CO2 reservoir using the measured petrophysical parameter of rock samples in the Gundih CCS Project, Indonesia | |
Xu et al. | Damage evolution characteristics of heterogeneous fractured sandstone reservoir under different fracturing fluids |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |