CN106872594A - 一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 - Google Patents
一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 Download PDFInfo
- Publication number
- CN106872594A CN106872594A CN201710082826.0A CN201710082826A CN106872594A CN 106872594 A CN106872594 A CN 106872594A CN 201710082826 A CN201710082826 A CN 201710082826A CN 106872594 A CN106872594 A CN 106872594A
- Authority
- CN
- China
- Prior art keywords
- overbar
- delta
- concentration distribution
- tau
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明涉及一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,本发明利用考虑油相膨胀的多孔介质一维轴向扩散模型,通过某一时间点多孔介质中的CO2实测浓度分布与数学模型计算得到的浓度分布进行拟合,从有量纲时间量和无量纲时间量的关系中求解出CO2在饱和油多孔介质中的扩散系数D。其中需要使用的实测浓度分布是在某一时间下,进行流体取样、分离并做气相色谱分析测得的。本发明能够直接测量得到某一时间点下,含油多孔介质中CO2的扩散浓度分布,避免了数模预测浓度的误差,同时求解方便,消除了压降曲线法在高温高压下压力急剧波动的影响以及复杂的求解、修正过程。
Description
技术领域
本发明涉及一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,属于油气田开发工程的技术领域。
背景技术
CO2溶于原油后,能够使原油体积膨胀,并有效减小原油粘度和油水间界面张力,是理想的注入流体,能大幅提高原油采收率。注CO2提高原油采收率技术的发展已有50多年的历史,目前已经得到了广泛应用。对于低渗、致密油藏,注水开发难度较大,注CO2开采更是有效的提高采收率方法。根据《石油与天然气杂志》的2012年“全球EOR调查”结果显示,在全球EOR项目中,气驱项目数量占54%,而CO2相关项目数量占到气驱项目的77%,已成为最为重要的提高原油采收率技术之一。
利用CO2提高采收率的关键是CO2在原油中大量溶解以改善原油物性,CO2在油藏中的扩散距离及扩散速度对CO2驱油、CO2吞吐作业的效果起决定性作用。而CO2在原油中的传质扩散是浓差作用下的自发过程,受到储层温度、压力、渗透率及含油饱和度等多种因素的影响,难以准确测量。
目前的研究中大多采取数值计算方法,对油藏中不同时间点、空间点的CO2扩散浓度进行预测,以《Industrial&Engineering Chemistry Research》杂志中2009年40期的《Experimental Study of Carbon Dioxide Diffusion in Oil-Saturated Porous Mediaunder Reservoir Conditions》一文采用的方法较为成熟。其利用压降曲线法通过数学模型描述了CO2在饱和原油的低渗多孔介质中的扩散,计算出多孔介质内部CO2无量纲浓度变化规律;并通过实测扩散压降曲线将无量纲浓度有量纲化,计算出了考虑原油膨胀现象的CO2扩散系数以及不同时间点处多孔介质中CO2浓度分布。这种方法也能模拟储层的高温高压环境,所得扩散系数更接近实际值。但是此方法的CO2浓度分布由纯计算获得,数学模型中涉及状态方程、压缩因子、膨胀系数等参数的计算,可能与油藏真实状态有较大误差,从而影响原油中CO2浓度的最终结果。
发明内容
针对现有技术存在的问题,本发明提供一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法。本发明能够对多孔介质介质内某处CO2浓度进行取样,直接测量获得某一时间点CO2浓度分布,并计算出CO2扩散系数。
本发明的技术方案如下:
一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,包括步骤如下:
1)对饱和好油样的岩心进行CO2轴向扩散;
2)扩散完毕后,对岩心轴向不同部位的流体进行取样,并记录实际扩散的有量纲时间t;
3)读出取得油样的体积Vo(m3),对取样得到的油样充分加热并收集所得气体,通过气相色谱分析获得CO2溶解量N(mol),计算得到CO2在岩心中的浓度分布mol/m3;
4)根据将测得的CO2浓度分布无量纲化,式②中c0为实验条件下CO2在油样中的最大溶解量,mol/m3;为CO2无量纲浓度;
5)化简一维轴向扩散的连续性方程得到式③,并进行无量纲化得到式④;
式③中,x为到岩心中点的距离,m;t为实际扩散时间,s;u为原油膨胀引起的轴向流动速度,m/s;D为CO2扩散系数,m2/s;无量纲化方法如式④所示:
式④中为无量纲距离;为无量纲流动速度;x0为岩心半长,m;τ为无量纲时间;
6)将无量纲化连续性方程结合边界条件形成数学模型⑤式,并通过对空间的中心差分和对时间的向前差分⑥式,形成差分方程⑦式;
差分格式如下:
式⑥中i为空间点序号,n为时间点序号,下同;
得到差分方程如下:
式⑦中,a、b、e分别为差分方程组中的系数,无实际物理含义;f为差分方程组中已知常数项,无实际物理含义;
7)通过高斯-赛德尔迭代求解微分方程组⑧,获得不同无量纲时间点下,岩心中的CO2无量纲浓度分布;
8)从所求解中确认达到实测无量纲浓度分布时的时间点τ,通过无量纲时间量与有量纲时间量的关系式求解得到实验条件下的CO2扩散系数D,式⑨中τ为无量纲时间;t为实际扩散的有量纲时间,s;x0为岩心半长,m。
根据本发明优选的,所述方法还包括前期对岩心饱和油的预处理:
(1-1)对岩心夹持器进行清洗并烘干;
(1-2)对长岩心进行抽真空,饱和油。
根据本发明优选的,所述方法还包括前期对岩心轴向扩散CO2的预处理:
(1-3)将岩心夹持器放入恒温箱,并调整温度至实验所需温度,稳定8-10h;
(1-4)将CO2通入岩心夹持器两端进行岩心轴向扩散,扩散过程中始终保持气体压力稳定在实验要求压力上;扩散过程中通过岩心夹持器对岩心侧面施加围压,围压保持高于扩散压力2~3MPa;
(1-5)扩散进行不少于10h,然后对岩心轴向不同部位的流体进行取样,并记录实际扩散时间t。
本发明的优点在于:
本发明能够直接测量得到某一扩散时间下,含油多孔介质中CO2的扩散浓度分布,避免了压降曲线法中数模预测浓度的误差,同时求解方便,消除了压降曲线法在高温高压下压力急剧波动的影响以及求解过程中的复杂修正过程,测试结果更精确。
本发明提出的方法可直接获得CO2浓度分布,进而通过无量纲时间量与有量纲时间量的关系,直接求得使求得CO2扩散系数,避免了压降曲线法中通过数学模型进行反复修正所带来的误差,精度更高且运算方便;同时本发明也可以满足高温高压实验条件,符合实际油藏环境,从而对CO2在多孔介质中的实际扩散情况进行研究,对现场应用给予理论指导。
附图说明
图1为本发明实施例中实测多孔介质在某一时间点的CO2浓度分布曲线;
图2为图1中数据点无量纲化后的无量纲浓度分布曲线;
图3为数学模型计算得到的无量纲浓度分布曲线。
具体实施方式
下面根据实施例和附图对本发明用法进行详述,但不限于此。
实施例1、
一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,包括步骤如下:
1)对饱和好油样的岩心进行CO2轴向扩散,首先对岩心饱和油的预处理:
(1-1)对岩心夹持器进行清洗并烘干;
(1-2)对长岩心进行抽真空,保持10MPa的压力饱和油。
其中,所述岩心为尺寸45×45×300mm的矩形岩心。
其次对岩心轴向扩散CO2的预处理:
(1-3)将岩心夹持器放入恒温箱,并调整温度至实验所需温度70℃,稳定8-10h;
(1-4)将CO2通入岩心夹持器两端进行岩心轴向扩散,扩散过程中始终保持气体压力稳定在实验要求压力上,此处优选15MPa;扩散过程中通过岩心夹持器对岩心侧面施加围压,围压保持高于扩散压力2~3MPa,阻止径向扩散;
(1-5)扩散进行不少于10h,然后对岩心轴向不同部位的流体进行取样,并记录实际扩散时间t,本实施例中,扩散进行13小时后,开启针阀取样器,对距离岩心中点距离分别为X0,X1,……,X5等6个部位的流体进行取样,并记录实际扩散时间t;扩散是从岩心两端同时进行的,因此两侧浓度分布关于岩心中点对称,只测一侧浓度即可。
2)扩散完毕后,对岩心轴向不同部位的流体进行取样,并记录实际扩散的有量纲时间t;
3)读出取得油样的体积Vo(m3),对取样得到的油样充分加热并收集所得气体,通过气相色谱分析获得CO2溶解量N(mol),计算得到CO2在岩心中的浓度分布mol/m3;
4)根据将测得的CO2浓度分布无量纲化,式②中c0为实验条件下CO2在油样中的最大溶解量,mol/m3;为CO2无量纲浓度;式②中c0为实验条件70℃,15MPa下CO2在油样中的最大溶解量,4788.24mol/m3;实测的浓度数据形成表1,并根据表1中的数据绘制成图1;
表1实验进行13小时后实测浓度分布数据
5)化简一维轴向扩散的连续性方程得到式③,并进行无量纲化得到式④;
式③中,x为到岩心中点的距离,m;t为实际扩散时间,s;u为原油膨胀引起的轴向流动速度,m/s;D为CO2扩散系数,m2/s;无量纲化方法如式④所示:
式④中为无量纲距离;为无量纲流动速度;x0为岩心半长,m;τ为无量纲时间;并离散求解,得到不同无量纲时间点下的浓度分布;根据式④中方法将图1中数据无量纲化,绘制成图2;
6)将无量纲化连续性方程结合边界条件形成数学模型⑤式,并通过对空间的中心差分和对时间的向前差分⑥式,形成差分方程⑦式;
差分格式如下:
式⑥中i为空间点序号,n为时间点序号,下同;
得到差分方程如下:
式⑦中,a、b、e分别为差分方程组中的系数,无实际物理含义;f为差分方程组中已知常数项,无实际物理含义;
7)通过高斯-赛德尔迭代求解微分方程组⑧,获得不同无量纲时间点下,岩心中的CO2无量纲浓度分布;
8)从所求解中确认达到实测无量纲浓度分布时的时间点τ,通过无量纲时间量与有量纲时间量的关系式求解得到实验条件下的CO2扩散系数D,式⑨中τ为无量纲时间;t为实际扩散的有量纲时间,s;x0为岩心半长,m。本实施例中,求解得到的不同时间点下无量纲浓度数据形成表2,并根据表2中数据绘制图3。由于CO2浓度分布关于岩心中点对称,只给出一侧的浓度分布;将图2与图3进行对比,找到图3中与图2重合度最好的一条曲线,得到该曲线所对应的无量纲时间。
从所求解中确认无量纲时间为τ=0.00075时,岩心中CO2浓度分布达到与实测值基本相同的状态。通过无量纲时间量与有量纲时间量的关系式求解得到实验条件下的CO2扩散系数:
表2通过数学模型计算得到的无量纲浓度分布数据
Claims (3)
1.一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,其特征在于,该方法包括步骤如下:
1)对饱和好油样的岩心进行CO2轴向扩散;
2)扩散完毕后,对岩心轴向不同部位的流体进行取样,并记录实际扩散的有量纲时间t;
3)读出取得油样的体积Vo(m3),对取样得到的油样充分加热并收集所得气体,通过气相色谱分析获得CO2溶解量N(mol),计算得到CO2在岩心中的浓度分布mol/m3;
4)根据将测得的CO2浓度分布无量纲化,式②中c0为实验条件下CO2在油样中的最大溶解量,mol/m3;为CO2无量纲浓度;
5)化简一维轴向扩散的连续性方程得到式③,并进行无量纲化得到式④;
式③中,x为到岩心中点的距离,m;t为实际扩散时间,s;u为原油膨胀引起的轴向流动速度,m/s;D为CO2扩散系数,m2/s;无量纲化方法如式④所示:
式④中为无量纲距离;为无量纲流动速度;x0为岩心半长,m;τ为无量纲时间;
6)将无量纲化连续性方程结合边界条件形成数学模型⑤式,并通过对空间的中心差分和对时间的向前差分⑥式,形成差分方程⑦式;
差分格式如下:
式⑥中i为空间点序号,n为时间点序号,下同;
得到差分方程如下:
式⑦中,a、b、e分别为差分方程组中的系数,无实际物理含义;f为差分方程组中已知常数项,无实际物理含义;
7)通过高斯-赛德尔迭代求解微分方程组⑧,获得不同无量纲时间点下,岩心中的CO2无量纲浓度分布;
8)从所求解中确认达到实测无量纲浓度分布时的时间点τ,通过无量纲时间量与有量纲时间量的关系式求解得到实验条件下的CO2扩散系数D,式⑨中τ为无量纲时间;t为实际扩散的有量纲时间,s;x0为岩心半长,m。
2.根据权利要求1所述的一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,其特征在于,所述方法还包括前期对岩心饱和油的预处理:
(1-1)对岩心夹持器进行清洗并烘干;
(1-2)对长岩心进行抽真空,饱和油。
3.根据权利要求1或2所述的一种测试饱和油多孔介质中CO2浓度分布和扩散系数的方法,其特征在于,所述方法还包括前期对岩心轴向扩散CO2的预处理:
(1-3)将岩心夹持器放入恒温箱,并调整温度至实验所需温度,稳定8-10h;
(1-4)将CO2通入岩心夹持器两端进行岩心轴向扩散,扩散过程中始终保持气体压力稳定在实验要求压力上;扩散过程中通过岩心夹持器对岩心侧面施加围压,围压保持高于扩散压力2~3MPa;
(1-5)扩散进行不少于10h,然后对岩心轴向不同部位的流体进行取样,并记录实际扩散时间t。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710082826.0A CN106872594B (zh) | 2017-02-16 | 2017-02-16 | 一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710082826.0A CN106872594B (zh) | 2017-02-16 | 2017-02-16 | 一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106872594A true CN106872594A (zh) | 2017-06-20 |
CN106872594B CN106872594B (zh) | 2018-11-16 |
Family
ID=59165997
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710082826.0A Active CN106872594B (zh) | 2017-02-16 | 2017-02-16 | 一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106872594B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108362614A (zh) * | 2018-01-19 | 2018-08-03 | 中国石油大学(华东) | 测量页岩油二氧化碳吞吐过程中扩散系数的装置及方法 |
CN108548835A (zh) * | 2018-02-09 | 2018-09-18 | 上海大学 | 一种评价泡沫在岩心驱油过程中动态稳定性的方法 |
CN113008739A (zh) * | 2021-03-23 | 2021-06-22 | 中国石油大学(华东) | 一种稠油热采条件下气体扩散系数的确定方法 |
CN113552024A (zh) * | 2021-06-24 | 2021-10-26 | 中国石油化工股份有限公司 | 针对co2在油水同时存在情况下的扩散传质的计算方法 |
CN113984587A (zh) * | 2021-10-29 | 2022-01-28 | 大连理工大学 | 一种原位测量多孔介质内co2-水扩散系数的方法 |
CN114859010A (zh) * | 2022-03-24 | 2022-08-05 | 中国矿业大学 | 监测气藏储层岩石注co2过程中co2波及动态的装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102121891A (zh) * | 2010-12-16 | 2011-07-13 | 西南石油大学 | 高温高压多组分油气体系分子扩散系数的测试装置及方法 |
CN102809526A (zh) * | 2012-08-27 | 2012-12-05 | 中国石油大学(华东) | 一种测量二氧化碳在饱和油岩心中扩散系数的方法 |
CN104502236A (zh) * | 2014-12-15 | 2015-04-08 | 中国石油大学(华东) | 一种测量co2从水相向油相扩散过程中扩散系数和平衡浓度的方法 |
CN106021778A (zh) * | 2016-05-31 | 2016-10-12 | 中国海洋石油总公司 | 一种模拟co2驱替动态混相压力的确定方法 |
-
2017
- 2017-02-16 CN CN201710082826.0A patent/CN106872594B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102121891A (zh) * | 2010-12-16 | 2011-07-13 | 西南石油大学 | 高温高压多组分油气体系分子扩散系数的测试装置及方法 |
CN102809526A (zh) * | 2012-08-27 | 2012-12-05 | 中国石油大学(华东) | 一种测量二氧化碳在饱和油岩心中扩散系数的方法 |
CN104502236A (zh) * | 2014-12-15 | 2015-04-08 | 中国石油大学(华东) | 一种测量co2从水相向油相扩散过程中扩散系数和平衡浓度的方法 |
CN106021778A (zh) * | 2016-05-31 | 2016-10-12 | 中国海洋石油总公司 | 一种模拟co2驱替动态混相压力的确定方法 |
Non-Patent Citations (4)
Title |
---|
SONGYAN LI 等,: "Diffusion coefficients of supercritical CO2 in oil-saturated cores under low permeability reservoir conditions", 《JOURNAL OF CO2 UTILIZATION》 * |
ZHAOWEN LI 等: "Experimental Study of Carbon Dioxide Diffusion in Oil-Saturated Porous Media under Reservoir Conditions", 《IND. ENG. CHEM. RES.》 * |
赵仁保 等: "CO2在原油中的扩散规律及变扩散系数计算方法", 《中国石油大学学报(自然科学版)》 * |
郭彪 等: "CO2在多孔介质中扩散系数的测定", 《石油化工高等学校学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108362614A (zh) * | 2018-01-19 | 2018-08-03 | 中国石油大学(华东) | 测量页岩油二氧化碳吞吐过程中扩散系数的装置及方法 |
CN108362614B (zh) * | 2018-01-19 | 2020-03-27 | 中国石油大学(华东) | 测量页岩油二氧化碳吞吐过程中扩散系数的装置及方法 |
CN108548835A (zh) * | 2018-02-09 | 2018-09-18 | 上海大学 | 一种评价泡沫在岩心驱油过程中动态稳定性的方法 |
CN108548835B (zh) * | 2018-02-09 | 2020-08-04 | 上海大学 | 一种评价泡沫在岩心驱油过程中动态稳定性的方法 |
CN113008739A (zh) * | 2021-03-23 | 2021-06-22 | 中国石油大学(华东) | 一种稠油热采条件下气体扩散系数的确定方法 |
CN113552024A (zh) * | 2021-06-24 | 2021-10-26 | 中国石油化工股份有限公司 | 针对co2在油水同时存在情况下的扩散传质的计算方法 |
CN113984587A (zh) * | 2021-10-29 | 2022-01-28 | 大连理工大学 | 一种原位测量多孔介质内co2-水扩散系数的方法 |
CN114859010A (zh) * | 2022-03-24 | 2022-08-05 | 中国矿业大学 | 监测气藏储层岩石注co2过程中co2波及动态的装置及方法 |
CN114859010B (zh) * | 2022-03-24 | 2023-04-21 | 中国矿业大学 | 监测气藏储层岩石注co2过程中co2波及动态的装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106872594B (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106872594B (zh) | 一种测试饱和油多孔介质中co2浓度分布和扩散系数的方法 | |
CN102644459B (zh) | 多组分气液体系在岩心中分子扩散系数的测定装置及方法 | |
CN108490156B (zh) | 高温高压条件下混合气体驱油埋存量实验测试方法 | |
Xiao et al. | Experimental investigation on CO2 injection in the Daqing extra/ultra-low permeability reservoir | |
CN103926267B (zh) | 一种定量评价应力敏感过程中孔喉变化程度的方法 | |
CN104359819A (zh) | 低渗致密岩心气水相对渗透率测定装置及测定方法 | |
CN109959672B (zh) | 一种二氧化碳置换页岩甲烷的定量评价方法及其应用 | |
CN105300849B (zh) | 一种多孔介质中气体扩散系数的测试装置及方法 | |
CN103674799A (zh) | 一种测定气体在多孔隙介质中轴向扩散系数的装置及方法 | |
CN102507414A (zh) | 地层压力条件下岩心渗透率实验测试方法及其装置 | |
CN202502089U (zh) | 成岩作用模拟实验装置 | |
CN105759015B (zh) | 一种在线测量注气原油体积系数和溶解度的装置及方法 | |
CN106066291A (zh) | Co2置换页岩气及页岩对页岩气或co2的吸附解析模拟测试系统与方法 | |
CN112198093B (zh) | 测试气体在饱和活油岩心中扩散系数的装置与方法 | |
CN114544463B (zh) | 一种裂缝性油藏地质体co2埋存评估方法及装置 | |
CN110879196B (zh) | 富含油凝析气藏油水相渗测试方法 | |
CN104389592B (zh) | 底水带油环凝析气藏水淹层油损失评价实验测试方法 | |
CN116150559A (zh) | 一种计算co2 在三相体系中的扩散系数的方法 | |
CN106840973A (zh) | 一种测试多孔介质中co2扩散浓度和扩散系数的装置及其测试方法 | |
CN105474011A (zh) | 确定油组合物与流体的最小混溶压力 | |
CN105548180B (zh) | 一种驱油用泡沫性能的在线检测方法 | |
CN105548465B (zh) | 一种驱油用泡沫性能的在线检测装置 | |
CN104316448B (zh) | 一种高阶煤岩气相渗透率动态变化的预测方法 | |
CN105929135B (zh) | 一种低渗透气藏储量动用程度的测试方法 | |
CN205620412U (zh) | 一种在线测量注气原油体积系数和溶解度的装置 |
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 |