CN109188525A - 一种海相页岩层埋深数据的获取方法及系统 - Google Patents

一种海相页岩层埋深数据的获取方法及系统 Download PDF

Info

Publication number
CN109188525A
CN109188525A CN201810858821.7A CN201810858821A CN109188525A CN 109188525 A CN109188525 A CN 109188525A CN 201810858821 A CN201810858821 A CN 201810858821A CN 109188525 A CN109188525 A CN 109188525A
Authority
CN
China
Prior art keywords
rammell
buried depth
data
earth
elevation
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
Application number
CN201810858821.7A
Other languages
English (en)
Other versions
CN109188525B (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.)
China University of Geosciences
Original Assignee
China University of Geosciences
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 China University of Geosciences filed Critical China University of Geosciences
Priority to CN201810858821.7A priority Critical patent/CN109188525B/zh
Publication of CN109188525A publication Critical patent/CN109188525A/zh
Application granted granted Critical
Publication of CN109188525B publication Critical patent/CN109188525B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/32Transforming one recording into another or one representation into another
    • G01V1/325Transforming one representation into another
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/40Transforming data representation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

一种海相页岩层埋深数据的获取方法及系统,本发明先分别以褶皱构造翼为单位建立地表出露地层单位(组)与其正下方目的页岩层埋深之间的拟合关系,然后将数字地表高程图与区域地质图叠置并通过数字化技术获取出露地层界线所在坐标的地表高程,再将地层界线的地表高程减去对应页岩层埋深的拟合值,获得各界线正下方目的页岩层的构造高度;将构造高度按照数字高程图相同的网格密度进行网格化处理得到高度数据网格,最后以数字地表高程网格减去页岩层构造高度网格,即最终获取到页岩层的埋深数据网格。在缺乏地震资料的地区,本发明能够进行页岩层地下埋深的推算,进而可编绘埋深图。

Description

一种海相页岩层埋深数据的获取方法及系统
技术领域
本发明涉及地质数据获取领域,更具体地说,涉及一种海相页岩层埋深数据的获取方法及系统,尤其应用在地震资料和钻井资料缺乏地区。
背景技术
继美国成为世界上实现页岩气大规模商业性开采的国家之后,我国当前正在大力加强页岩气研究与勘探,并已经建成涪陵、长宁-威远等国家级页岩气产业示范区。我国中上扬子区四川盆地及周缘海相地层中上奥陶统五峰组-下志留统龙马溪组、下寒武统筇竹寺组两套富有机质页岩分布广泛,是目前海相页岩气勘探的主要目的层。页岩层埋深是页岩气勘探及评价的重要参数。但在缺乏地震资料的地区,无法通过地震资料解释地下构造,因此难以准确获取目的层的埋深。对于这类地区,如何有效推测页岩层地下埋深,是一项非常有意义的技术。国内一些学者曾在页岩气的研究中编绘过四川盆地及周缘主要页岩层系的埋深图,但由于广大区域缺乏三维、甚至二维地震资料,埋深图的精度并不高,埋深等值线间距大达500m、甚至1000m,显然主要是仅依据地层分布区域估算的。鲜有考虑地表高程起伏因素而进行更精确推算者。
综上所述,对于地震资料缺乏区,如何充分利用区域地质图有关信息,同时充分结合地表高程数据资料,相对准确的获取海相页岩层的埋深数据,是一项十分有意义的技术,但此技术有待进一步开发。
发明内容
本发明要解决的技术问题在于,针对现有技术中在缺乏地震资料的地区上获取的页岩层埋深精度低、未考虑地表高程因素的技术缺陷,提供一种海相页岩层埋深数据的获取方法及系统,来获取海相页岩层埋深数据。
根据本发明的其中一方面,本发明解决其技术问题所采用的技术方案是:提供一种地震资料缺乏区海相页岩层埋深数据的获取方法,包含如下步骤:
S1、分别建立待数据获取区域的图切横剖面的各个构造翼地表出露地层单位与其正下方页岩层埋深之间的拟合关系,并根据多条剖面建立图幅区域总拟合关系;
S2、将所述待数据获取区域的数字地表高程图与该待数据获取区域的区域地质图叠置,通过数字化技术,获取区域地质图上出露地层界线所在坐标的地表高程;
S3、将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,获得该坐标目的页岩层的构造高度;合并研究区内所有构造翼页岩层的构造高度数据,得到研究区总的页岩层构造高度数据;
S4、将研究区页岩层总的构造高度数据按照数字高程图相同的网格密度进行网格化处理,即得到页岩层构造高度网格数据;再将地表高程网格数据减去页岩层构造高度网格数据,即获得页岩层的埋深网格数据,从而最终获取海相页岩层的埋深数据。
进一步地,在本发明的海相页岩层埋深数据的获取方法中,所述其他有效地质信息包括:地层界线、地层产状、断层、经纬度。
进一步地,在本发明的海相页岩层埋深数据的获取方法中,步骤S1中所述的图切横剖面,包括区域地质图上原有的横剖面图、根据需要新制作的横剖面图。
进一步地,在本发明的海相页岩层埋深数据的获取方法中,步骤S2中所述数字地表高程图为DEM(Digital Elevation Model)数据。
进一步地,在本发明的海相页岩层埋深数据的获取方法中,步骤S3中,将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,是指每一个构造翼具有独自一套埋深拟合值;对于无法判属构造翼的局部地区,用待数据获取区域的总拟合关系确定出露地层对应的拟合值。
根据本发明的另一方面,本发明为解决其技术问题还提供了一种地震资料缺乏区海相页岩层埋深数据的获取系统,包含如下模块:
页岩层埋深数据拟合模块,用于分别建立待数据获取区域的图切横剖面的各个构造翼地表出露地层单位与其正下方页岩层埋深之间的拟合关系,并根据多条剖面建立图幅区域总拟合关系;
地层界线地表高程获取模块,用于将所述待数据获取区域的数字地表高程图与该待数据获取区域的区域地质图叠置,通过数字化技术,获取区域地质图上出露地层界线所在坐标的地表高程;
页岩层构造高度数据获取模块,用于将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,获得该坐标页岩层的构造高度;合并研究区内所有构造翼页岩层构造高度数据,得到研究区总的页岩层构造高度数据;
页岩层埋深网格数据获取模块,用于将研究区页岩层总的构造高度数据按照数字高程图相同的网格密度进行网格化处理,即得到页岩层构造高度网格数据;再将地表高程网格数据减去该构造高度网格数据,即获得页岩层的埋深网格数据,从而最终获取海相页岩层的埋深数据。
进一步地,在本发明的海相页岩层埋深数据的获取系统中,所述其他有效地质信息包括:地层界线、地层产状、断层、经纬度。
进一步地,在本发明的海相页岩层埋深数据的获取系统中,页岩层埋深数据拟合模块中所述的图切横剖面,包括区域地质图上原有的横剖面图、根据需要新制作的横剖面图。
进一步地,在本发明的海相页岩层埋深数据的获取系统中,地层界线地表高程获取模块中所述数字地表高程图为DEM(Digital Elevation Model)数据。
进一步地,在本发明的海相页岩层埋深数据的获取系统中,页岩层构造高度数据获取模块中,将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,是指每一个构造翼具有独自一套埋深拟合值;对于无法判属构造翼的局部地区,用待数据获取区域的总拟合关系确定出露地层对应的拟合值。
实施本发明的海相页岩层埋深数据的获取方法及系统,具有以下有益效果:在缺乏地震资料的地区,本发明能够进行页岩层埋深的推算,且推算出的埋深数据精度较高,并可继续相关平面图的编图,从而为页岩气的初步勘探评价提供重要依据。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明的海相页岩层埋深数据的获取方法一实施例的流程图;
图2是地质图的图切横剖面图及其目的层埋深读取示意图;
图3是目的层埋深同地表出露地层(组)之间的拟合关系图;
图4数字高程等值线图与区域地质图重叠示意图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
参考图1,其为本发明的海相页岩层埋深数据的获取方法一实施例的流程图。本实施例依据地表地质信息及地形高程数据,来获取地震资料缺乏区厚度较为稳定的海相地层富有机质页岩段的地下埋深。可以按以下步骤实施:
S1、对于缺乏地震资料的研究区,无法根据地震资料解释地下岩层的构造特征与埋深。但可以根据该区地质图的图切横剖面图及其他有效地质信息,按照海相页岩层上覆海相地层的厚度趋势与产状变化,来恢复页岩层地下构造形态,并可分别度量地表出露地层界线正下方的页岩层埋深。出露地层越新、其下方页岩层埋深越大,出露地层越老其下方页岩层埋深越小。参考图2,根据待数据获取区的地质图横剖面图及其他有效地质信息,按照目的页岩层上覆海相地层的厚度趋势与产状变化,恢复页岩层地下构造形态,并分别度量地表出露地层界线正下方的页岩层埋深(diF,i=1,2,..n,褶皱构造的翼部数;F=a,b,..q,一翼出露的地层组)。参考图3,再分别建立图切剖面的各个构造翼地表出露地层单位(组)与其正下方页岩层埋深之间的拟合关系,并根据多条剖面建立图幅区域总拟合关系,参考图3;在本实施例中,其他有效地质信息包括地质界面、地层产状及经纬度等。
S2、将所述待数据获取区域的将数字地表高程图(DEM,Digital ElevationModel)与该待数据获取区域的区域地质图叠置,通过数字化技术,获取区域地质图上出露组界线所在的地表高程[hF(x,y,z)];(x:经度;y:纬度;z:海拔/m),参考图4,其为数字高程图与区域地质图重叠示意图。
S3、将出露组界线所在的地表高程(海拔)减去对应的岩石地层单位的目的层地下真埋深的拟合值(hF-diF),即获得各界线正下方目的页岩层的构造高度(海拔)[lF(x,y,z)];合并所有出露地层组(F=a,b,..q)构造高度数据,得到目的层地下构造高度数据[L(x,y,z)]。
S4、将研究区页岩层总的构造高度数据按数字高程图的网格密度进行网格化处理,再将高程网格数据减去构造高度网格数据,即[H(x,y,z)-L(x,y,z)],即获得地下沉积岩层真埋深数据[D(x,y,z)]的网格。
本发明在没有地震资料的地区利用了广泛覆盖的区域地质图及其横切地质剖面图上的信息,能够较好地控制埋深小于4000m页岩层埋深;且本发明利用了DEM(DigitalElevation Model)数据,这类数据精度高,对于求取页岩层的真埋深,减少了地表海拔的误差,提高了精度;本发明可以利用比例尺较大的地质图(只要有地质图就可以自己做剖面,就可以实施方案)进行实施,提高埋深数据的精度。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的。本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。

Claims (10)

1.一种地震资料缺乏区海相页岩层埋深数据的获取方法,其特征在于,包含如下步骤:
S1、分别建立待数据获取区域的图切横剖面的各个构造翼地表出露地层单位与其正下方页岩层埋深之间的拟合关系,并根据多条剖面建立图幅区域总拟合关系;
S2、将所述待数据获取区域的数字地表高程图与该待数据获取区域的区域地质图叠置,通过数字化技术,获取区域地质图上出露地层界线所在坐标的地表高程;
S3、将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,获得该坐标目的页岩层的构造高度;合并研究区内所有构造翼页岩层的构造高度数据,得到研究区总的页岩层构造高度数据;
S4、将研究区页岩层总的构造高度数据按照数字高程图相同的网格密度进行网格化处理,即得到页岩层构造高度网格数据;再将地表高程网格数据减去页岩层构造高度网格数据,即获得页岩层的埋深网格数据,从而最终获取海相页岩层的埋深数据。
2.根据权利要求1所述的海相页岩层埋深数据的获取方法,其特征在于,所述其他有效地质信息包括:地层界线、地层产状、断层、经纬度。
3.根据权利要求1所述的海相页岩层埋深数据的获取方法,其特征在于,步骤S1中所述的图切横剖面,包括区域地质图上原有的横剖面图、根据需要新制作的横剖面图。
4.根据权利要求1所述的海相页岩层埋深数据的获取方法,其特征在于,步骤S2中所述数字地表高程图为DEM数据,即Digital Elevation Model数据。
5.根据权利要求1所述的海相页岩层埋深数据的获取方法,其特征在于,步骤S3中,将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,是指每一个构造翼具有独自一套埋深拟合值;对于无法判属构造翼的局部地区,用待数据获取区域的总拟合关系确定出露地层对应的拟合值。
6.一种地震资料缺乏区海相页岩层埋深数据的获取系统,其特征在于,包含如下模块:
页岩层埋深数据拟合模块,用于分别建立待数据获取区域的图切横剖面的各个构造翼地表出露地层单位与其正下方页岩层埋深之间的拟合关系,并根据多条剖面建立图幅区域总拟合关系;
地层界线地表高程获取模块,用于将所述待数据获取区域的数字地表高程图与该待数据获取区域的区域地质图叠置,通过数字化技术,获取区域地质图上出露地层界线所在坐标的地表高程;
页岩层构造高度数据获取模块,用于将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,获得该坐标页岩层的构造高度;合并研究区内所有构造翼页岩层构造高度数据,得到研究区总的页岩层构造高度数据;
页岩层埋深网格数据获取模块,用于将研究区页岩层总的构造高度数据按照数字高程图相同的网格密度进行网格化处理,即得到页岩层构造高度网格数据;再将地表高程网格数据减去该构造高度网格数据,即获得页岩层的埋深网格数据,从而最终获取海相页岩层的埋深数据。
7.根据权利要求6所述的海相页岩层埋深数据的获取方法,其特征在于,所述其他有效地质信息包括:地层界线、地层产状、断层、经纬度。
8.根据权利要求6所述的海相页岩层埋深数据的获取方法,其特征在于,页岩层埋深数据拟合模块中所述的图切横剖面,包括区域地质图上原有的横剖面图、根据需要新制作的横剖面图。
9.根据权利要求6所述的海相页岩层埋深数据的获取方法,其特征在于,地层界线地表高程获取模块中所述数字地表高程图为DEM数据,即Digital Elevation Model数据。
10.根据权利要求6所述的海相页岩层埋深数据的获取方法,其特征在于,页岩层构造高度数据获取模块中,将出露地层界线所在坐标的地表高程,减去该坐标所在构造翼地层对应页岩层埋深的拟合值,是指每一个构造翼具有独自一套埋深拟合值;对于无法判属构造翼的局部地区,用待数据获取区域的总拟合关系确定出露地层对应的拟合值。
CN201810858821.7A 2018-07-31 2018-07-31 一种海相页岩层埋深数据的获取方法及系统 Expired - Fee Related CN109188525B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810858821.7A CN109188525B (zh) 2018-07-31 2018-07-31 一种海相页岩层埋深数据的获取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810858821.7A CN109188525B (zh) 2018-07-31 2018-07-31 一种海相页岩层埋深数据的获取方法及系统

Publications (2)

Publication Number Publication Date
CN109188525A true CN109188525A (zh) 2019-01-11
CN109188525B CN109188525B (zh) 2019-09-17

Family

ID=64937555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810858821.7A Expired - Fee Related CN109188525B (zh) 2018-07-31 2018-07-31 一种海相页岩层埋深数据的获取方法及系统

Country Status (1)

Country Link
CN (1) CN109188525B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112554876A (zh) * 2019-09-26 2021-03-26 中国石油天然气集团有限公司 地层凹陷区域选取方法及装置
CN114384585A (zh) * 2021-12-30 2022-04-22 西北核技术研究所 基于相对位置及最小埋深的山地地下爆炸绝对定位方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017003517A1 (en) * 2015-06-30 2017-01-05 Halliburton Energy Services, Inc. Correcting the effects of deviation and dispersion on sonic log measurements of deviated wells in laminated formations
CN107102354A (zh) * 2016-12-21 2017-08-29 中国石油化工股份有限公司江汉油田分公司物探研究院 一种页岩甜点地震综合评价方法
CN107590513A (zh) * 2017-09-11 2018-01-16 西南石油大学 一种页岩气富集成藏分类方法
US20180016896A1 (en) * 2016-07-15 2018-01-18 Geoisochem Corporation Assessing Permeability
CN107817520A (zh) * 2017-09-20 2018-03-20 中国石油化工股份有限公司 海相泥页岩地层的压力系数预测方法及系统
CN107966546A (zh) * 2017-11-21 2018-04-27 西南石油大学 一种页岩岩相平面分布编制方法及页岩勘探体系

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017003517A1 (en) * 2015-06-30 2017-01-05 Halliburton Energy Services, Inc. Correcting the effects of deviation and dispersion on sonic log measurements of deviated wells in laminated formations
US20180016896A1 (en) * 2016-07-15 2018-01-18 Geoisochem Corporation Assessing Permeability
CN107102354A (zh) * 2016-12-21 2017-08-29 中国石油化工股份有限公司江汉油田分公司物探研究院 一种页岩甜点地震综合评价方法
CN107590513A (zh) * 2017-09-11 2018-01-16 西南石油大学 一种页岩气富集成藏分类方法
CN107817520A (zh) * 2017-09-20 2018-03-20 中国石油化工股份有限公司 海相泥页岩地层的压力系数预测方法及系统
CN107966546A (zh) * 2017-11-21 2018-04-27 西南石油大学 一种页岩岩相平面分布编制方法及页岩勘探体系

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
朱定伟 等: "鄂尔多斯盆地东南部古构造恢复及地质意义", 《特种油气藏》 *
武瑾 等: "渝东北地区下志留统龙马溪组页岩气勘探前景", 《特种油气藏》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112554876A (zh) * 2019-09-26 2021-03-26 中国石油天然气集团有限公司 地层凹陷区域选取方法及装置
CN114384585A (zh) * 2021-12-30 2022-04-22 西北核技术研究所 基于相对位置及最小埋深的山地地下爆炸绝对定位方法

Also Published As

Publication number Publication date
CN109188525B (zh) 2019-09-17

Similar Documents

Publication Publication Date Title
Nowroozi et al. Saltwater intrusion into the freshwater aquifer in the eastern shore of Virginia: a reconnaissance electrical resistivity survey
Blaikie et al. Three-dimensional potential field modelling of a multi-vent maar-diatreme—The Lake Coragulac maar, Newer Volcanics Province, south‐eastern Australia
Sweet Interaction between aeolian, fluvial and playa environments in the Permian Upper Rotliegend Group, UK southern North Sea
Mancini et al. A physical stratigraphy model for seismic microzonation of the Central Archaeological Area of Rome (Italy)
Miller et al. Back to basics of sequence stratigraphy: early Miocene and mid-Cretaceous examples from the New Jersey paleoshelf
Lindenmaier et al. Structure and genesis of the Cubango Megafan in northern Namibia: implications for its hydrogeology
Li et al. Seismic analysis of polygonal fault systems in the Great South Basin, New Zealand
McFarlan Jr et al. Subsurface geology of the late Tertiary and Quaternary deposits, coastal Louisiana and the adjacent continental shelf
Gallerini et al. 3D modeling using geognostic data: The case of the low valley of Foglia river (Italy)
CN109188525B (zh) 一种海相页岩层埋深数据的获取方法及系统
Souei et al. Hydrogeological framework and geometry modeling via joint gravity and borehole parameters, the Nadhour-Sisseb-El Alem basin (central-eastern Tunisia)
Livani et al. Subsurface geological and geophysical data from the Po Plain and the northern Adriatic Sea (north Italy)
Bosch et al. Airborne electromagnetic measurements as basis for a 3D geological model of an Elsterian incision
Aizebeokhai et al. Field evaluation of 3D geo-electrical resistivity imaging for environmental and engineering studies using parallel 2D profiles
Sweetkind Three-dimensional hydrogeologic framework model of the Rio Grande transboundary region of New Mexico and Texas, USA, and northern Chihuahua, Mexico
Avalos et al. Two-dimensional seismic refraction tomography of a buried bedrock valley at Hallsands beach, Devon, United Kingdom
Caran et al. Quaternary stratigraphy and paleoenvironments of the Texas Rolling Plains
Lachaal et al. Seismic, gravity, and wire line logging characterization of the Zéramdine fault corridor and its influence in the deep Miocene aquifers distribution (east-central Tunisia)
Wylie Jr et al. Log-curve amplitude slicing: Visualization of log data and depositional trends in the Middle Devonian Traverse Group, Michigan basin, United States
Fadillah et al. Hydrogeological interpretation using electrical resistivity tomography: methodology and conceptual model in andesitic volcanic context
Sweetkind et al. Construction of 3-D geologic framework and textural models for Cuyama Valley groundwater basin, California
Khaled et al. Geoelectrical and hydrogeological study to delineate the geological structures affecting the groundwater occurrence in Wadi El Khariq Basin, Northwest El Maghara, north Sinai, Egypt
Stamos-Pfeiffer et al. Hydrogeology and simulation of groundwater flow in the Lucerne Valley groundwater basin, California
Lau et al. Modeling the sediment fill of the upper Troy pre-glacial Bedrock Valley, McHenry County, Illinois, USA
Sweetkind et al. Geologic Characterization of Young Alluvial Basin-Fill Deposits from Drill Hole Data in Yucca Flat, Nye County, Nevada.

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
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: 20190917

Termination date: 20210731