CN110954463B - Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 - Google Patents

Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 Download PDF

Info

Publication number
CN110954463B
CN110954463B CN201911277298.XA CN201911277298A CN110954463B CN 110954463 B CN110954463 B CN 110954463B CN 201911277298 A CN201911277298 A CN 201911277298A CN 110954463 B CN110954463 B CN 110954463B
Authority
CN
China
Prior art keywords
corrosion
cement
well cementation
area
permeability
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
Application number
CN201911277298.XA
Other languages
English (en)
Other versions
CN110954463A (zh
Inventor
冯福平
丛子渊
韩旭
严茂森
刘圣源
胡小康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeast Petroleum University
Original Assignee
Northeast Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northeast Petroleum University filed Critical Northeast Petroleum University
Priority to CN201911277298.XA priority Critical patent/CN110954463B/zh
Publication of CN110954463A publication Critical patent/CN110954463A/zh
Application granted granted Critical
Publication of CN110954463B publication Critical patent/CN110954463B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Abstract

本发明涉及的是CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,它包括:一、根据油井或水井被腐蚀处固井水泥石所对应的实际地层温压条件、CO2含量、离子含量及固井水泥石参数,通过数值模拟得出不同时间固井水泥石内各腐蚀区域厚度的模拟结果以及未腐蚀区域厚度;二、计算固井水泥石各腐蚀区域内每种矿物不同时间的体积分数;三、计算固井水泥石各腐蚀区域的孔隙度;四、计算固井水泥石各腐蚀区域的渗透率;五、结合步骤一模拟结果,计算不同腐蚀阶段垂直腐蚀方向水泥石的等效渗透率。本发明考虑了水泥石试件尺寸对垂直腐蚀方向等效渗透率的影响,能够真实评价CO2腐蚀对井下不同尺寸环空水泥石垂直腐蚀方向等效渗透率的影响程度。

Description

CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法
技术领域:
本发明涉及的是水泥环完整性评价技术,具体涉及的是CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法。
背景技术:
含CO2气藏及CO2驱注井均面临固井水泥石被腐蚀的问题,含CO2的酸性流体与固井水泥石发生化学反应,改变水泥的微观物质组成及宏观力学性质,降低环空水泥环的密封完整性,极易造成CO2逃逸现象的发生。一方面,逃逸的CO2会破坏注采井网的完整性及污染地下水水质;另一方面,还存在CO2气体沿水泥环空泄漏至地面的可能性,造成井口环空带压,增加了泄压成本和安全生产的风险性。
对于CO2腐蚀对水泥石物性的影响,国内外学者均进行了大量的研究。 Ashok通过实验发现由于Ca(OH)2的溶解和反应速率要高于CSH,因此在反应过程中主要是Ca(OH)2与H+结合生成致密的膨胀性腐蚀产物CaCO3,导致腐蚀后水泥石的渗透率减小,孔隙度降低;然而,在Hunt、Kim等人的研究中认为,在富含CO2的地层条件下,随着碳酸溶液的不断侵蚀,初期反应生成的致密性腐蚀产物CaCO3在H2O及CO2的作用下会转变为易溶性的Ca(HCO3)2,并不断消耗水泥石内部的Ca(OH)2形成淋滤作用,使得水泥石的孔隙度和渗透率增大,抗压强度降低。CO2腐蚀对水泥渗透率的影响主要通过实验测得,但目前的实验研究存在如下问题:
1、CO2腐蚀对水泥石渗透率变化规律的影响认识不清,不同学者得到的结果甚至相反。由于CO2对水泥石的腐蚀受到实验温度、离子浓度、压力、水泥石原始渗透率、流体性质、水泥类型等多因素的影响,致使CO2与水泥石的反应速率是不同的,具体表现为其孔隙度、渗透率以及抗压强度的变化随腐蚀时间的不同呈现较大的差异,导致不同学者在实验中得出的CO2对水泥石腐蚀不同时间后孔隙度、渗透率以及抗压强度变化规律相反的结论,使得CO2腐蚀对水泥石渗透率变化规律的影响认识不清。
2、现有CO2腐蚀后水泥石渗透率的实验评价方法忽略了各腐蚀区域及水泥石的试件尺寸对整体渗透率的影响。CO2与固井水泥石发生长期化学反应后,由水泥石表面向内侧将依次形成无定型SiO2富集区、CaCO3沉积区、Ca(OH)2消耗区和未反应区4个区域。由于不同腐蚀区内物质组成和微观结构不同,导致各腐蚀区的渗透特性存在较大差异。目前CO2腐蚀对水泥石渗透率的影响主要通过腐蚀反应后实验测得,然而腐蚀区域深度相对于实验水泥石试件较小,所测出的渗透率结果受实验水泥石试件的尺寸影响较大。在相同的实验条件下,水泥石尺寸越大,所测得腐蚀对其整体渗透率的影响就越小。这种评价方法由于实验条件和水泥石试件尺寸的不统一,无法真实的评价CO2腐蚀对水泥石渗透率的影响程度。
3、CO2对水泥石的腐蚀速率相对较慢而且影响因素较大,受实验条件及经费的限制,只能得出部分条件下短时间内的实验结果,不能全面的反映不同井下工况条件以及长期腐蚀后水泥石等效渗透率的变化规律。
4、现有的研究无法描述CO2腐蚀后垂直于腐蚀方向水泥石等效渗透率的动态变化。水泥石受CO2腐蚀后垂直腐蚀方向的渗透率发生改变,而垂直于腐蚀方向上的渗透率决定了水泥环沿井筒方向上的封隔能力,垂直腐蚀方向水泥石渗透率的大小是评价CO2腐蚀后水泥环完整性的重要参数。由于腐蚀区域的厚度较小,目前的室内实验无法进行垂直于腐蚀方向的渗透率测定;理论研究方面,也尚未提出垂直腐蚀方向水泥石等效渗透率的评价方法,因此现有的研究无法描述CO2腐蚀后垂直于腐蚀方向水泥石等效渗透率的动态变化。
含CO2气藏及CO2驱注井中,CO2对水泥石的腐蚀不可避免。真实准确的评价不同CO2腐蚀阶段垂直腐蚀方向水泥石的等效渗透率是井筒完整性评价和补救措施确定的重要依据,然而目前国内外尚未提出垂直于CO2腐蚀方向水泥石等效渗透率的评价方法。因此,亟需一种科学的评价方法解决上述问题,揭示不同CO2腐蚀阶段垂直腐蚀方向水泥石等效渗透率的动态变化规律。
发明内容:
本发明的目的是提供CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,这种CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法用于解决由于腐蚀区域的厚度较小,目前的室内实验无法进行垂直于腐蚀方向的渗透率测定的问题。
本发明解决其技术问题所采用的技术方案是:这种CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法:
步骤一、根据油井或水井被腐蚀处固井水泥石所对应的实际地层温压条件、 CO2含量、离子含量及固井水泥石参数,通过数值模拟得出不同时间固井水泥石内各腐蚀区域厚度的模拟结果hi以及未腐蚀区域厚度h4,i为腐蚀区域编号,i=1, 2,3,h1为无定型SiO2富集区、h2为CaCO3沉积区,h3为Ca(OH)2消耗区;
步骤二、计算固井水泥石各腐蚀区域内每种矿物不同时间的体积分数:
φj,i,t+Δt=φj,i,t+Vj×rj,i,t+Δt×Δt
式中:j为水泥中矿物种类编号;φj,i,t为第j种矿物在i区域t时间的体积分数; Vj为矿物j的摩尔体积;rj,i,t+Δt为j物质在i区域t时间的反应速率;Δt为时间步长;
步骤三、计算固井水泥石各腐蚀区域的孔隙度φi
Figure BDA0002315904120000031
式中:Nj为腐蚀区域矿物的数量;φj,i,t+Δt为各种矿物的体积分数;
步骤四、计算固井水泥石各腐蚀区域的渗透率ki
Figure BDA0002315904120000032
式中:φ为固井水泥石腐蚀前各腐蚀区域的孔隙度、φi为固井水泥石腐蚀后各腐蚀区域的孔隙度;k为固井水泥石腐蚀前各腐蚀区域的渗透率,ki为固井水泥石腐蚀后各腐蚀区域的渗透率;
步骤五、结合步骤一模拟结果,计算不同腐蚀阶段垂直腐蚀方向水泥石的等效渗透率:
Figure BDA0002315904120000033
式中:h为固井水泥石试件长度,k1为固井水泥石腐蚀后无定型SiO2富集区的渗透率、k2为固井水泥石腐蚀后CaCO3沉积区的渗透率,k3固井水泥石腐蚀后为Ca(OH)2消耗区的渗透率。
本发明具有以下有益效果:
(1)本发明提供的CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法主要基于数值模拟的结果,只需进行少量的室内实验进行验证即可,节省了大量的实验时间和实验成本。
(2)本发明提供的CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,可以得出不同腐蚀时间对垂直腐蚀方向水泥石等效渗透率的影响,能够进行较长腐蚀时间下的垂直腐蚀方向水泥石等效渗透率动态评价。
(3)本发明提供的CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,可以进行不同腐蚀条件下的垂直腐蚀方向水泥石等效渗透率评价,能够通过对比分析真实准确的得出影响垂直腐蚀方向水泥石等效渗透率的主要因素。
(4)本发明提供的CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,考虑了水泥石试件尺寸对垂直腐蚀方向等效渗透率的影响,能够更为真实的评价CO2腐蚀对井下不同尺寸环空水泥石垂直腐蚀方向等效渗透率的影响程度。
(5)本发明提供的CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,填补了CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价的空白,为定量评价CO2腐蚀后水泥石的密封完整性提供了理论依据。
具体实施方式
下面对本发明做进一步的说明:
这种CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法:
步骤一、根据油井或水井被腐蚀处固井水泥石所对应的实际地层温压条件、 CO2含量、离子含量及固井水泥石参数,通过数值模拟得出不同时间固井水泥石内各腐蚀区域厚度的模拟结果hi以及未腐蚀区域厚度h4,i为腐蚀区域编号,i=1, 2,3,h1为无定型SiO2富集区、h2为CaCO3沉积区,h3为Ca(OH)2消耗区;
步骤二、计算固井水泥石各腐蚀区域内每种矿物不同时间的体积分数:
φj,i,t+Δt=φj,i,t+Vj×rj,i,t+Δt×Δt
式中:j为水泥中矿物种类编号;φj,i,t为第j种矿物在i区域t时间的体积分数; Vj为矿物j的摩尔体积,单位为m3·mol-1;rj,i,t+Δt为j物质在i区域t时间的反应速率,单位为mol·m-3·s-1;Δt为时间步长,单位为s;
步骤三、计算固井水泥石各腐蚀区域的孔隙度φi
Figure BDA0002315904120000051
式中:Nj为腐蚀区域矿物的数量;φj,i,t+Δt为各种矿物的体积分数;
步骤四、计算固井水泥石各腐蚀区域的渗透率ki
Figure BDA0002315904120000052
式中:φ为固井水泥石腐蚀前各腐蚀区域的孔隙度、φi为固井水泥石腐蚀后各腐蚀区域的孔隙度;k为固井水泥石腐蚀前各腐蚀区域的渗透率,单位为10-3um2, ki为固井水泥石腐蚀后各腐蚀区域的渗透率,单位为10-3um2,其中i=1,2,3, k1为固井水泥石腐蚀后无定型SiO2富集区的渗透率、k2为固井水泥石腐蚀后 CaCO3沉积区的渗透率,k3固井水泥石腐蚀后为Ca(OH)2消耗区的渗透率。
步骤五、结合步骤一模拟结果,计算不同腐蚀阶段垂直腐蚀方向水泥石的等效渗透率:
Figure BDA0002315904120000053
式中:h为固井水泥石试件长度,单位为mm,k1为固井水泥石腐蚀后无定型SiO2富集区的渗透率、k2为固井水泥石腐蚀后CaCO3沉积区的渗透率,k3固井水泥石腐蚀后为Ca(OH)2消耗区的渗透率,h1为无定型SiO2富集区、h2为 CaCO3沉积区,h3为Ca(OH)2消耗区,h4为未腐蚀区域厚度。
这种等效渗透率评价方法考虑了不同CO2腐蚀阶段水泥石内部微观腐蚀区域的特性及试样尺寸对垂直腐蚀方向等效渗透率的影响,更为真实准确的反映了CO2腐蚀后垂直腐蚀方向水泥石等效渗透率的变化规律,同时能够显著的降低室内实验的时间和成本。

Claims (1)

1.一种CO2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法,其特征在于:
步骤一、根据油井或水井被腐蚀处固井水泥石所对应的实际地层温压条件、CO2含量、离子含量及固井水泥石参数,通过数值模拟得出不同时间固井水泥石内各腐蚀区域厚度的模拟结果hi以及未腐蚀区域厚度h4,i为腐蚀区域编号,i=1,2,3,h1为无定型SiO2富集区、h2为CaCO3沉积区,h3为Ca(OH)2消耗区;
步骤二、计算固井水泥石各腐蚀区域内每种矿物不同时间的体积分数:
φj,i,t+Δt=φj,i,t+Vj×rj,i,t+Δt×Δt
式中:j为水泥中矿物种类编号;φj,i,t为第j种矿物在i区域t时间的体积分数;Vj为矿物j的摩尔体积;rj,i,t+Δt为j物质在i区域t时间的反应速率;Δt为时间步长;
步骤三、计算固井水泥石各腐蚀区域的孔隙度φi
Figure FDA0003554841400000011
式中:Nj为腐蚀区域矿物的数量;φj,i,t+Δt为各种矿物的体积分数;
步骤四、计算固井水泥石各腐蚀区域的渗透率ki
Figure FDA0003554841400000012
式中:φ为固井水泥石腐蚀前各腐蚀区域的孔隙度、φi为固井水泥石腐蚀后各腐蚀区域的孔隙度;k为固井水泥石腐蚀前各腐蚀区域的渗透率,ki为固井水泥石腐蚀后各腐蚀区域的渗透率;
步骤五、结合步骤一模拟结果,计算不同腐蚀阶段垂直腐蚀方向水泥石的等效渗透率:
Figure FDA0003554841400000013
式中:h为固井水泥石试件长度,k1为固井水泥石腐蚀后无定型SiO2富集区的渗透率、k2为固井水泥石腐蚀后CaCO3沉积区的渗透率,k3为固井水泥石腐蚀后Ca(OH)2消耗区的渗透率,k4为固井水泥石腐蚀后未腐蚀区域的渗透率。
CN201911277298.XA 2019-12-12 2019-12-12 Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法 Active CN110954463B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911277298.XA CN110954463B (zh) 2019-12-12 2019-12-12 Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911277298.XA CN110954463B (zh) 2019-12-12 2019-12-12 Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法

Publications (2)

Publication Number Publication Date
CN110954463A CN110954463A (zh) 2020-04-03
CN110954463B true CN110954463B (zh) 2022-04-22

Family

ID=69981292

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911277298.XA Active CN110954463B (zh) 2019-12-12 2019-12-12 Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法

Country Status (1)

Country Link
CN (1) CN110954463B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5156205A (en) * 1991-07-08 1992-10-20 Prasad Raj K Method of determining vertical permeability of a subsurface earth formation
AU1722992A (en) * 1991-06-27 1993-03-04 Schlumberger Technology B.V. Apparatus for determining horizontal and/or vertical permeability of an earth formation
CN102830057A (zh) * 2012-09-10 2012-12-19 西南石油大学 一种固井水泥石腐蚀评价方法
CN103091234A (zh) * 2011-11-03 2013-05-08 中国石油天然气股份有限公司 一种评价酸性气藏气井固井水泥环腐蚀状况的方法
CN104391342A (zh) * 2014-12-09 2015-03-04 中国石油天然气股份有限公司 基于裂缝孔隙度反演定量预测储层渗透性的方法及装置
CN106354983A (zh) * 2016-10-20 2017-01-25 中国石油化工股份有限公司 一种确定co2埋存泄漏风险监测点的方法
CN106644921A (zh) * 2016-10-26 2017-05-10 东北石油大学 用于二氧化碳腐蚀试验的岩芯夹持器
CN107366534A (zh) * 2017-08-10 2017-11-21 中国石油天然气股份有限公司 粗化渗透率的确定方法和装置

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU1722992A (en) * 1991-06-27 1993-03-04 Schlumberger Technology B.V. Apparatus for determining horizontal and/or vertical permeability of an earth formation
US5156205A (en) * 1991-07-08 1992-10-20 Prasad Raj K Method of determining vertical permeability of a subsurface earth formation
CN103091234A (zh) * 2011-11-03 2013-05-08 中国石油天然气股份有限公司 一种评价酸性气藏气井固井水泥环腐蚀状况的方法
CN102830057A (zh) * 2012-09-10 2012-12-19 西南石油大学 一种固井水泥石腐蚀评价方法
CN104391342A (zh) * 2014-12-09 2015-03-04 中国石油天然气股份有限公司 基于裂缝孔隙度反演定量预测储层渗透性的方法及装置
CN106354983A (zh) * 2016-10-20 2017-01-25 中国石油化工股份有限公司 一种确定co2埋存泄漏风险监测点的方法
CN106644921A (zh) * 2016-10-26 2017-05-10 东北石油大学 用于二氧化碳腐蚀试验的岩芯夹持器
CN107366534A (zh) * 2017-08-10 2017-11-21 中国石油天然气股份有限公司 粗化渗透率的确定方法和装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CO2 对水泥石腐蚀机理及密封性的影响研究进展;冯福平 等;《硅酸盐学报》;20180228;第46卷(第2期);第247-255页 *
CO2埋存过程中储层岩石孔渗变化规律研究;韩旭;《中国优秀硕博士学位论文全文数据库(硕士) 工程科技Ⅰ辑》;20190115(第1期);第13-45页 *
Development of a Novel Method To Evaluate Well Integrity During CO2 Underground Storage;Mingxing Bai 等;《SPE Journal》;20150630;第628-641页 *
Fully Coupled Modeling of Long Term Cement Well Seal Stability in the Presence of CO2;James C Wilson 等;《Energy Procedia》;20111231;第4卷;第5162-5169页 *
Wellbore Cement Porosity Evolution in Response to Mineral Alteration during CO2 Flooding;Michael C. Cheshire 等;《Environmental Science & Technology》;20161213;第51卷;第692-698页 *

Also Published As

Publication number Publication date
CN110954463A (zh) 2020-04-03

Similar Documents

Publication Publication Date Title
Mangat et al. Prediction of long term chloride concentration in concrete
Ducasse-Lapeyrusse et al. Effect of calcium gluconate, calcium lactate, and urea on the kinetics of self-healing in mortars
CN110954462B (zh) Co2腐蚀后沿腐蚀方向水泥石等效渗透率评价方法
CN105388249B (zh) 一种测试酸蚀裂缝导流能力的装置及方法
Reiterman et al. Crystalline coating and its influence on the water transport in concrete
Costa et al. Evaluation of density influence on resistance to carbonation process in oil well cement slurries
CN101701899B (zh) 流水和压力水作用下水泥石溶出性侵蚀测试装置及测试方法
Duraisamy Strength and stiffness improvement of bio-cemented Sydney sand
Huang et al. Micromechanical properties of biocemented shale soils analyzed using nanoindentation test
Adeoye et al. Physical and chemical alterations in engineered cementitious composite under geologic CO2 storage conditions
Tan et al. Experimental investigation on the effects of different fracturing fluids on shale surface morphology
Gao et al. Experimental study on repair of fractured rock mass by microbial induction technology
CN110954463B (zh) Co2腐蚀后垂直腐蚀方向水泥石等效渗透率评价方法
Condor et al. Experimental results of diffusion coefficient of sulfate ions in cement type 10 and class G
Joshi et al. Performance of bacterial mediated mineralization in concrete under carbonation and chloride induced corrosion
Turnbull et al. Effectiveness of corrosion inhibitors in retarding the rate of propagation of localized corrosion
Yujie et al. Experimental study on the reaction process and engineering characteristics of marine calcareous sand reinforced by eco-friendly methods
Lv et al. Experimental study on the durability and microstructure of marine concrete covered with barnacles
Roux et al. Influence of a biopolymer admixture on corrosion behaviour of steel rebars in concrete
Dong et al. Effect of oxygen on CO2 corrosion and erosion-corrosion behavior of N80 steel under high temperature and high pressure
CN109932498B (zh) 一种分析酸蚀裂缝壁面稳定性的实验方法
Ma et al. Study of the Initial Corrosion of X80 Steel in a Saturated Saline Soil Co-Contaminated with Cl-, SO42-and HCO3
Turner et al. Treatment of fractured concrete via microbially induced carbonate precipitation: From micro-scale characteristics to macro-scale behaviour
Niu et al. Corrosion rules for ordinary concrete exposed to sulfur dioxide-containing environments
Inazumi et al. Assessment of potassium silicate based surface penetration materials with low viscosity in the repair of concrete structures

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