CN103487830B - 一种地震振幅类属性勘探中补偿差异压实的方法 - Google Patents
一种地震振幅类属性勘探中补偿差异压实的方法 Download PDFInfo
- Publication number
- CN103487830B CN103487830B CN201310421916.XA CN201310421916A CN103487830B CN 103487830 B CN103487830 B CN 103487830B CN 201310421916 A CN201310421916 A CN 201310421916A CN 103487830 B CN103487830 B CN 103487830B
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- compaction
- sample spot
- numerical value
- 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.)
- Expired - Fee Related
Links
Abstract
本发明公开了一种地震振幅类属性勘探中补偿差异压实的方法,该方法的步骤包括:对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释,沿地震解释层提取地震属性平面图;统计样品点坐标与样品点之间的深度差异数据;按x,y大地坐标与所提取坐标寻求对应关系;得到在一个地震属性数值随深度增加而增加的关系式;运用得到的关系式对导出的地震属性数据体进行二次运算,得到补偿压实作用之后的平面地震属性矢量图。本发明通过这种补偿压实作用对地震属性工作影响的方法,最大限度的体现了地层由于岩性变化所引起阻抗值的变化,排除由于构造压实作用影响岩层密度从而引起的反演结果不准确,从而能真实反映地层发育特征与分布范围。
Description
技术领域
本发明属于地球物理技术领域,尤其涉及一种地震振幅类属性勘探中补偿差异压实的方法。
背景技术
根据目前国内外已知的补偿压实的方法,整体上说只是一个通过缩小研究工区面积降低构造压实对岩层密度的影响,而目前这类方法并不能真正起到真实反映波阻抗值由岩性变化而变化的趋势,存在相对较大的误差。如果只是传统的勘探方法,不仅投资大,而且延误工期,甚至有时因场地无法布置钻探而影响勘探精度。因此单一地质勘察方法已不能满足需求,需要应用新的勘察方法,以求的更高的经济效益。
地震属性是指那些由叠前或叠后地震数据经过数学变换而导出的有关地震波的几何形态、运动学特征和统计学特征的特殊测量值。根据不同的算法、不同的属性具有不同的地质意义,有些属性值可能擅长于揭示不易探测到的岩性变化,而有些属性则可以直接用于烃类检测。
地震属性的研究和使用始于20世纪70年代,最早主要是以振幅为基础的瞬时属性,用来直接指示油气与地下不同岩层的展布情况。70年代早期,NigelAnstey发现了含气砂岩波阻抗的异常变化,使用了反射波振幅变化特征-“亮点”与“暗点”对含气砂岩储集体进行预测。80年代,地震属性的数量迅速增加,多元属性分析技术首次被使用。80年代晚期,多维属性如倾角和方位角的初步发展,导致了三维连续属性在90年代的广泛应用。究其原因,它的成功得益于清楚和明确的地质含义。90年代以来,基于精细储层描述以及三维地震数据体解释的需要,地震属性分析技术更是急剧发展。现今,地震属性技术已广泛应用于储层预测、油气藏动态监测、油气藏特征描述等领域,并取得了很好的效果。
综合地震属性的数学、物理及地质意义,也可以将地震属性分类两类:第一类是基本属性,包括反射时间、振幅及相位,他们各自具有独立的物理意义和地质意义;第二类是延伸属性,如对多种属性进行处理(叠合、求差,相关、求主成分等)得到的复合属性以及根据单一地震属性再提取出的新的级联属性。其中振幅类属性主要是用于体现地下不同岩性地层的分布特征与展布规律。然而传统的勘探技术不易探测到岩性的变化,且不能直接用于烃类检测,对于寻找剩余油、提高采收率成为油田开发工作的重点及难点,也是稳产的关键。
发明内容
本发明的目的在于提供一种地震振幅类属性勘探中补偿差异压实的方法,旨在解决传统的勘探方法,不仅投资大,而且延误工期,甚至有时因场地无法布置钻探而影响勘探精度的问题。
本发明是这样实现的,一种地震振幅类属性勘探中补偿差异压实的方法,该方法的步骤包括:
步骤一、对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释,沿地震解释层提取地震属性平面图;
根据测井、录井数据完成井上层位标定;
根据井上层位标定结合地震反射特征完成地震解释;
根据地震数据体进行属性体提取的二次运算;
根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片,得到沿地震解释层位地震属性平面图。
步骤二、在对研究区进行岩心观察、录井解释的基础上进行小层对比,挑选相同岩石类型,相同地层沉积厚度但埋深不同的样品点,统计样品点坐标与样品点之间的深度差异数据;
根据对取心井观察与数据统计,挑选相同岩石类型与沉积厚度但不同埋深的样品点;
根据所述样品点选择,统计其大地坐标和样品点埋深数据;
根据所述埋深数值二次运算得到埋深数据之间的差值。
步骤三、将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的地震属性数值,按x,y大地坐标与步骤二中所提取样品点坐标寻求对应关系;
将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的包括大地坐标的地震属性数值;
根据所述地震属性数值中的大地坐标和井位坐标寻求相同坐标的样品点,将这部分样品点作为有效样品点挑选出来;
根据所挑选的有效样品点,计算其地震属性数值差异。
步骤四、根据工区研究对象精度要求,设定补偿压实作用差值最小计算单元,然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析,得到在一个地震属性数值随深度增加而增加的关系式;
根据研究要求,设定补偿差异压实作用最小计算单元;
根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析,得到一个关于两组差值数据之间的关系式:Y=aX+b,其中Y为地震属性差异数值,X为深度差异数值,a、b为根据交汇分析所得常量;
步骤五、运用得到的关系式对导出的地震属性数据体进行二次运算,对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值,然后重新平面成图,从而得到补偿压实作用之后的沿地震解释层位的平面地震属性矢量图。
进一步,沿地震解释层提取地震属性平面图中的反演工作要求较为精确的地震层位以及断层解释,然后合理的合成记录制作以及时深转换;
进一步,统计样品点坐标与样品点之间的深度差异数据中数据提取工作要求合理的小层对比方案,尽可能多的提取数据点,保证样品点的代表性与非偶然性;
进一步,沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值中寻求数据对应的工作要求阻抗数值的大地坐标与井点坐标两者尽可能的接近甚至重合,要求平面上距离最短,需考虑方位角和偏移距等可能引起误差的各种因素。方位角和偏移距有误差使研究区沿地震解释层的波阻抗数值不精确,且导致阻抗数值的大地坐标与井点坐标两者出现偏差,不能达到平面上距离最短的要求;
进一步,为得到在一个地震属性数值随深度增加而增加的关系式需根据研究目标精度要求设定数据单元,且数据单元应该不小于研究目标要求精度;
进一步,为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰,尽可能的反映由于岩性变化而引起的属性变化空间展布和变化规律。
本发明提供的地震振幅类属性勘探中补偿差异压实的方法,该实验方法真实性高,通过这种补偿压实作用对地震属性工作影响的方法,最大限度的体现了地层由于岩性变化所引起阻抗值的变化,排除由于构造压实作用影响岩层密度从而引起的反演结果不准确,从而能真实反映地层发育特征与分布范围。
本发明提供的地震振幅类属性勘探中补偿差异压实的方法,该实验工作效率高,与目前国内外主流的研究区分块,分段地震解释后提取属性的方法相比(目前国内一般是20平方公里,美国一般是20平方英里),现有方法在针对较大工区整体工作时具有无与伦比的效率优势和技术优势。首先,整个工区集中在一起研究能够从总体上把握整个研究区的构造与沉积背景,在地质上有一个整体的认识,从而避免了从局部无法把握全局的场面;其次,在地震属性提取的过程中整个工区一起进行能够最大限度的避免对地震体分块而导致的数据丢失,保证数据的完整性与精度;再次,集中属性提取能够在曲线归一化的基础上尽可能多的运用井资料和地震资料结合,相互影响,控制因素较多,最大限度的保证了研究工作的质量和准确性;最后,整体属性分析能够避免最后全区综合评价时期的一个拼接工作,提高工作效率。
附图说明
图1是本发明提供的地震振幅类属性勘探中补偿差异压实的方法。
具体实施方式
本发明是这样实现的,一种地震振幅类属性勘探中补偿差异压实的方法,该方法的步骤包括:
在步骤S101中:对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释,沿地震解释层提取地震属性平面图;
根据测井、录井数据完成井上层位标定;
根据井上层位标定结合地震反射特征完成地震解释;
根据地震数据体进行属性体提取的二次运算;
根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片,得到沿地震解释层位地震属性平面图。
在步骤S102中:在对研究区进行岩心观察、录井解释的基础上进行小层对比,挑选相同岩石类型,相同地层沉积厚度但埋深不同的样品点,统计样品点坐标与样品点之间的深度差异数据;
根据对取心井观察与数据统计,挑选相同岩石类型与沉积厚度但不同埋深的样品点;
根据所述样品点选择,统计其大地坐标和样品点埋深数据;
根据所述埋深数值二次运算得到埋深数据之间的差值。
在步骤S103中:将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的地震属性数值,按x,y大地坐标与步骤二中所提取样品点坐标寻求对应关系;
将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的包括大地坐标的地震属性数值;
根据所述地震属性数值中的大地坐标和井位坐标寻求相同坐标的样品点,将这部分样品点作为有效样品点挑选出来;
根据所挑选的有效样品点,计算其地震属性数值差异。
在步骤S104中:根据工区研究对象精度要求,设定补偿压实作用差值最小计算单元,然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析,得到在一个地震属性数值随深度增加而增加的关系式;
根据研究要求,设定补偿差异压实作用最小计算单元;
根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析,得到一个关于两组差值数据之间的关系式:Y=aX+b,其中Y为地震属性差异数值,X为深度差异数值,a、b为根据交汇分析所得常量;
在步骤S105中:运用得到的关系式对导出的地震属性数据体进行二次运算,对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值,然后重新平面成图,从而得到补偿压实作用之后的沿地震解释层位的平面地震属性矢量图。
进一步,沿地震解释层提取地震属性平面图中的反演工作要求较为精确的地震层位以及断层解释,然后合理的合成记录制作以及时深转换;
进一步,统计样品点坐标与样品点之间的深度差异数据中数据提取工作要求合理的小层对比方案,尽可能多的提取数据点,保证样品点的代表性与非偶然性;
进一步,沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值中寻求数据对应的工作要求阻抗数值的大地坐标与井点坐标两者尽可能的接近甚至重合,要求平面上距离最短,需考虑方位角和偏移距等可能引起误差的各种因素。方位角和偏移距有误差使研究区沿地震解释层的波阻抗数值不精确,且导致阻抗数值的大地坐标与井点坐标两者出现偏差,不能达到平面上距离最短的要求。
进一步,为得到在一个地震属性数值随深度增加而增加的关系式需根据研究目标精度要求设定数据单元,且数据单元应该不小于研究目标要求精度。
进一步,为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰,尽可能的反映由于岩性变化而引起的属性变化空间展布和变化规律。
本发明提供的地震振幅类属性勘探中补偿差异压实的方法,该方法包括:对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释,沿地震解释层提取地震属性平面图;在对研究区进行岩心观察、录井解释的基础上进行小层对比,挑选相同岩石类型,相同地层沉积厚度但埋深不同的样品点,统计样品点坐标与样品点之间的深度差异数据;将沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值,按x,y大地坐标与所提取坐标寻求对应关系;根据工区研究对象精度要求,设定补偿压实作用差值单元,然后将所对应好的数据进行交汇分析,得到在一个地震属性数值随深度增加而增加的关系式;运用得到的关系式对导出的地震属性数据体进行二次运算,对全区不同深度点的地震属性数据乘以不同系数之后得出新的属性值,然后重新平面成图,从而得到补偿压实作用之后的平面地震属性矢量图。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种地震振幅类属性勘探中补偿差异压实的方法,其特征在于,该方法的步骤包括:
步骤一、对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释,沿地震解释层提取地震属性平面图;
根据测井、录井数据完成井上层位标定;
根据井上层位标定结合地震反射特征完成地震解释;
根据地震数据体进行属性体提取的二次运算;
根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片,得到沿地震解释层位地震属性平面图;
步骤二、在对研究区进行岩心观察、录井解释的基础上进行小层对比,挑选相同岩石类型,相同地层沉积厚度但埋深不同的样品点,统计样品点坐标与样品点之间的深度差异数据;
根据对取心井观察与数据统计,挑选相同岩石类型与沉积厚度但不同埋深的样品点;
根据所述样品点选择,统计其大地坐标和样品点埋深数据;
根据所述埋深数值二次运算得到埋深数据之间的差值;
步骤三、将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的地震属性数值,按x,y大地坐标与步骤二中所提取样品点坐标寻求对应关系;
将沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的包括大地坐标的地震属性数值;
根据所述地震属性数值中的大地坐标和井位坐标寻求相同坐标的样品点,将这部分样品点作为有效样品点挑选出来;
根据所挑选的有效样品点,计算其地震属性数值差异;
步骤四、根据工区研究对象精度要求,设定补偿压实作用差值最小计算单元,然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析,得到在一个地震属性数值随深度增加而增加的关系式;
根据研究要求,设定补偿差异压实作用最小计算单元;
根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析,得到一个关于两组差值数据之间的关系式:Y=aX+b,其中Y为地震属性差异数值,X为深度差异数值,a、b为根据交汇分析所得常量;
步骤五、运用得到的关系式对导出的地震属性数据体进行二次运算,对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值,然后重新平面成图,从而得到补偿压实作用之后的沿地震解释层位的平面地震属性矢量图。
2.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法,其特征在于,沿地震解释层提取地震属性平面图,该步工作要求地震层位以及断层解释,然后合成记录制作以及时深转换公式。
3.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法,其特征在于,统计样品点坐标与样品点之间的深度差异数据,该步数据提取工作要求合理的小层对比方案;提取样品的时候需要在不同岩性区域,不同构造单元中提取样品点。
4.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法,其特征在于,沿地震解释层位提取的属性图按数值格式导出,得到一组研究区沿地震解释层的地震属性数值,然后将该组数值中的大地坐标寻求与权利要求1步骤2中的样品大地坐标进行对应,该步工作要求地震属性数值的大地坐标与样品坐标两者尽可能的接近甚至重合,要求平面上距离最短,需考虑样品坐标中方位角和偏移距引起的误差。
5.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法,其特征在于,为得到在一个地震振幅类属性数值随深度增加而增加的关系式,需根据研究目标精度要求设定数据单元,且数据单元应该不小于研究目标要求精度。
6.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法,其特征在于,为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰,尽可能的反映由于岩性变化而引起的振幅类属性变化的空间展布和变化规律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310421916.XA CN103487830B (zh) | 2013-09-17 | 2013-09-17 | 一种地震振幅类属性勘探中补偿差异压实的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310421916.XA CN103487830B (zh) | 2013-09-17 | 2013-09-17 | 一种地震振幅类属性勘探中补偿差异压实的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487830A CN103487830A (zh) | 2014-01-01 |
CN103487830B true CN103487830B (zh) | 2016-07-06 |
Family
ID=49828199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310421916.XA Expired - Fee Related CN103487830B (zh) | 2013-09-17 | 2013-09-17 | 一种地震振幅类属性勘探中补偿差异压实的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487830B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103852787B (zh) * | 2014-02-24 | 2016-08-24 | 长江大学 | 一种砂岩储层成岩地震相表征方法 |
CN104834008B (zh) * | 2015-05-04 | 2017-10-17 | 中国石油天然气股份有限公司 | 计算碳酸盐岩缝洞型储层充填程度的方法 |
CN106569277B (zh) * | 2015-10-12 | 2018-08-07 | 中国石油化工股份有限公司 | 地震属性叠后振幅补偿方法 |
CN112860820B (zh) * | 2019-11-27 | 2023-12-26 | 中国石油天然气集团有限公司 | 一种通用地球物理数据存储方法及装置 |
CN112965112B (zh) * | 2021-02-25 | 2022-06-10 | 中国海洋石油集团有限公司 | 地下储层微构造三维散点的生成方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
EP1914568A2 (en) * | 2006-10-20 | 2008-04-23 | BHP Billiton Innovation Pty Ltd | Method for determining impedence coefficients of a seismic trace |
CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
CN102841378A (zh) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | 一种地震反演数据预测储层的方法 |
-
2013
- 2013-09-17 CN CN201310421916.XA patent/CN103487830B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6058073A (en) * | 1999-03-30 | 2000-05-02 | Atlantic Richfield Company | Elastic impedance estimation for inversion of far offset seismic sections |
EP1914568A2 (en) * | 2006-10-20 | 2008-04-23 | BHP Billiton Innovation Pty Ltd | Method for determining impedence coefficients of a seismic trace |
CN102116871A (zh) * | 2011-01-24 | 2011-07-06 | 中国石油天然气股份有限公司 | 一种提高地震反演数据预测储层能力的方法 |
CN102841378A (zh) * | 2011-06-23 | 2012-12-26 | 中国石油天然气股份有限公司 | 一种地震反演数据预测储层的方法 |
Non-Patent Citations (3)
Title |
---|
AVO约束稀疏脉冲反演技术在濮城油田储层预测中的应用;张铭,等;《石油天然气学报》;20120430;第34卷(第4期);第78-82页 * |
稀疏脉冲和基于模型反演在王家岗沙四段砂岩油藏精细勘探中的应用;刘成斋,等;《石油地球物理勘探》;20020831;第37卷(第4期);第401-406页 * |
霸县凹陷文安斜坡沙二段-东营组地震反演研究;李进;《中国优秀硕士学位论文全文数据库 基础科学辑》;20130115(第01期);第20-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103487830A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
McCarthy et al. | Anatomy of a metamorphic core complex: Seismic refraction/wide‐angle reflection profiling in southeastern California and western Arizona | |
CN103487830B (zh) | 一种地震振幅类属性勘探中补偿差异压实的方法 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103454686A (zh) | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 | |
CN104678434A (zh) | 一种预测储层裂缝发育参数的方法 | |
CN103217714A (zh) | 基于频率-尺度匹配的地震沉积学解释方法 | |
CN104793247A (zh) | 近地表速度建模方法及建模装置 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN104834008B (zh) | 计算碳酸盐岩缝洞型储层充填程度的方法 | |
Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
CN105093311A (zh) | 叠合盆地多期差异隆升区地层剥蚀厚度测量方法 | |
CN108121009B (zh) | 开发后期复杂断块变速构造成图方法 | |
CN103257363A (zh) | 一种探测地下裂缝性储层中裂缝倾角的方法 | |
CN103399345B (zh) | 一种潜山裂缝分布的勘测方法与装置 | |
Lin et al. | Predicting the azimuth of natural fractures and in situ horizontal stress: A case study from the Sichuan Basin, China | |
CN105445787A (zh) | 一种最优方位子体相干的裂缝预测方法 | |
CN103529473B (zh) | 一种地震反演过程中差异压实的补偿方法 | |
YU et al. | Constrained joint inversion of gravity and seismic data using the simulated annealing algorithm | |
Wilson et al. | Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site | |
Jiskani et al. | Fault orientation modeling of Sonda-Jherruck coalfield, Pakistan | |
Chen et al. | Description technology of fractured-vuggy carbonate reservoir in halahatang oil field, tarim basin—take the ha 7 test area as an example | |
Deng et al. | Hydrocarbon accumulation conditions and key exploration and development technologies for PL 19–3 oilfield | |
Zhang et al. | Architecture characteristics and characterization methods of fault-controlled karst reservoirs: A case study of the Shunbei 5 fault zone in the Tarim Basin, China | |
Nikitenko et al. | High-frequency LWD electromagnetic measurements applied to studying macroanisotropic properties of thin-bedded sedimentary deposits | |
Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20180917 |