CN104007485A - 一种建立复杂孔隙泥页岩岩石物理模型方法 - Google Patents

一种建立复杂孔隙泥页岩岩石物理模型方法 Download PDF

Info

Publication number
CN104007485A
CN104007485A CN201410153372.8A CN201410153372A CN104007485A CN 104007485 A CN104007485 A CN 104007485A CN 201410153372 A CN201410153372 A CN 201410153372A CN 104007485 A CN104007485 A CN 104007485A
Authority
CN
China
Prior art keywords
kerogen
hole
aspect ratio
pores
modulus
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
CN201410153372.8A
Other languages
English (en)
Other versions
CN104007485B (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 Petroleum Beijing
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201410153372.8A priority Critical patent/CN104007485B/zh
Publication of CN104007485A publication Critical patent/CN104007485A/zh
Application granted granted Critical
Publication of CN104007485B publication Critical patent/CN104007485B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明一种建立复杂孔隙泥页岩岩石物理模型方法,基于改进的Xu-Payne理论的复杂孔隙泥页岩岩石物理模型,该模型通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布,引入有机质分布,并将孔隙类型分为基质孔隙和干酪根孔隙,基质孔隙分为粒内孔、粒间孔和裂缝,分别用孔隙纵横比0.8、0.15、0.01来表征,干酪根孔隙为微裂隙且充填油气,利用孔隙纵横比0.01来表征,干酪根与油气混合物通过Brown-Korringa方程来实现固体替代,并且假设其纵横比为0.8,从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。

Description

一种建立复杂孔隙泥页岩岩石物理模型方法
技术领域
本发明属于岩石物理技术领域,特别涉及一种建立复杂孔隙泥页岩岩石物理模型的方法,是一种基于改进Xu-Payne理论的考虑了复杂基质孔隙结构及干酪根微裂隙影响建立的可以预测含有机质饱和流体岩石弹性张量的岩石物理模型。
背景技术
泥页岩由于其复杂的孔隙特征,以及有机质的分布形态及特征,传统的岩石物理模型既没有考虑复杂的基质孔隙形态,更没有考虑有机质的存在及其分布特征,实验表明有机质的存在会显著影响岩石的弹性性质,且有机质的分布形态及特征也是评价页岩气气源条件的重要因素,而且镜下观测表明有机质可能以充填物的方式存在于孔隙中,且干酪根中也存在充填油气的微裂隙,迫切需要研发一种能够模拟有机质及复杂孔隙响应的岩石物理模型。
发明内容
本发明的目的在于提供一种建立复杂孔隙泥页岩岩石物理模型方法,通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布,引入有机质分布,并将孔隙类型分为基质孔隙和干酪根孔隙,基质孔隙分为粒内孔、粒间孔和裂缝,分别用孔隙纵横比0.8、0.15、0.01来表征,干酪根孔隙为微裂隙且充填油气,利用孔隙纵横比0.01来表征,干酪根与油气混合物通过Brown-Korringa方程来实现固体替代,并且假设其纵横比为0.8,从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。
为了实现上述目的,本发明一种建立复杂孔隙泥页岩岩石物理模型方法,包括以下步骤:
步骤1,结合测井曲线,岩心数据及地质地化报告,得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征;
步骤2,根据步骤1得到的矿物组分的体积分数及弹性参数,通过Voigt-Reuss-Hill平均公式(1952)得到岩石骨架的体积模量和剪切模量; 
步骤3,将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙,通过DEM(Differential Effective Medium)模型得到干酪根与干酪根孔隙的混合物,设定干酪根孔隙的纵横比为0.01,并通过Gassmann(1951)方程在微裂隙中充满油气;
步骤4,通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙,设定其纵横比为0.8,并通过Brown-Korringa(1975)方程实现固体替代;
步骤5,利用DEM模型向前面形成的混合物中加入基质孔隙,分为粒内孔、粒间孔、裂缝,得到干岩石的体积模量与剪切模量;
步骤6,通过Gassmann(1951)方程进行流体替代,计算饱和流体岩石的体积模量与剪切模量;
步骤7,利用计算的体积模量与剪切模量,结合计算的总密度和岩石速度之间的关系,得到泥页岩的纵波速度和横波速度,并与实测数据进行对比,判断是否吻合,如果吻合,结束,如果不吻合,返回步骤5,调整基质孔隙中不同孔隙形态所占的体积分数,直至吻合。
本发明与已有技术相比达到的有益效果是:
基于改进的Xu-Payne(2009)理论的复杂孔隙泥页岩岩石物理模型,该模型通过改进Xu-Payne模型针对碳酸盐储层复杂孔隙特征分布,引入有机质分布,并将孔隙类型分为基质孔隙和干酪根孔隙,基质孔隙分为粒内孔、粒间孔和裂缝,分别用孔隙纵横比0.8、0.15、0.01来表征,干酪根孔隙为微裂隙且充填油气,利用孔隙纵横比0.01来表征,干酪根与油气混合物通过Brown-Korringa方程来实现固体替代,并且假设其纵横比为0.8,从而建立考虑有机质与复杂孔隙类型影响的泥页岩岩石物理模型。
附图说明
图1是基于改进Xu-Payne理论的复杂孔隙泥页岩岩石物理模型建模流程图。
图2是根据测井曲线计算得到的干酪根与地化分析得到的干酪根数据准确度分析图。
图3是根据本模型预测得到的纵波速度与实测纵波速度准确度分析图。
图4(a)是含有机质岩石纵波阻抗随干酪根体积分数的变化示意图。
图4(b)是含有机质岩石纵波阻抗与横波阻抗随干酪根体积分数的变化示意图。
图5(a)是含有机质岩石纵波阻抗与纵横波速度比随干酪根体积分数的变化示意图。
图5(b)是含有机质岩石杨氏模量与脆性因子随干酪根体积分数的变化示意图。
图6基于改进Xu-Payne理论的复杂孔隙形态泥页岩岩石物理模型得到的页岩气井综合分析评价图。
具体实施方式
以下结合实例与附图说明本发明具体实施方式。
一种建立复杂孔隙泥页岩岩石物理模型方法,包括以下步骤:
步骤1,结合测井曲线,岩心数据及地质地化报告,得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征;
步骤2,根据步骤1得到的矿物组分的体积分数及弹性参数,通过Voigt-Reuss-Hill平均公式(1952)得到岩石骨架的体积模量和剪切模量; 
步骤3,将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙,通过DEM(Differential Effective Medium)模型得到干酪根与干酪根孔隙的混合物,设定干酪根孔隙的纵横比为0.01,并通过Gassmann(1951)方程在微裂隙中充满油气;
步骤4,通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙,设定其纵横比为0.8,并通过Brown-Korringa(1975)方程实现固体替代;
步骤5,利用DEM模型向前面形成的混合物中加入基质孔隙,分为粒内孔、粒间孔、裂缝,得到干岩石的体积模量与剪切模量;
步骤6,通过Gassmann(1951)方程进行流体替代,计算饱和流体岩石的体积模量与剪切模量;
步骤7,利用计算的体积模量与剪切模量,结合计算的总密度和岩石速度之间的关系,得到泥页岩的纵波速度和横波速度,并与实测数据进行对比,判断是否吻合,如果吻合,结束,如果不吻合,返回步骤5,调整基质孔隙中不同孔隙形态所占的体积分数,直至吻合。
实施例中,所述吻合的误差范围为10%。
实施例中,当实际测井曲线中只有纵波测井曲线时,用纵波速度进行约束。
实施例中,步骤3所述的孔隙度分为基质孔隙与干酪根孔隙,基质孔隙分为粒内孔、粒间孔、裂缝,分别用纵横比0.8、0.15、0.01来表征,干酪根孔隙为充填其中的微裂隙,设定纵横比为0.01,通过DEM模型实现干酪根与微裂隙的混合,并通过Gassmann实现微裂隙中充满油气。
实施例中,步骤4所述干酪根与油气混合物与基质矿物的混合,通过DEM模型加入干孔隙,设定纵横比为0.8,并通过Brown-Korringa方程来实现固体替代,加入干酪根与油气混合物。
以上具体实施方式仅用于说明本发明,而非用于限定本发明。
下面通过一口页岩气井纵横波速度预测与优势层段特征分析来说明:
图1是基于改进Xu-Payne理论的复杂孔隙泥页岩岩石物理模型建模流程图。
步骤1,结合测井曲线,岩心数据及地质地化报告,得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征;
步骤2,根据步骤1得到的矿物组分的体积分数及弹性参数,通过Voigt-Reuss-Hill平均公式(1952)得到岩石骨架的体积模量和剪切模量; 
步骤3,将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙,通过DEM(Differential Effective Medium)模型得到干酪根与干酪根孔隙的混合物,设定干酪根孔隙的纵横比为0.01,并通过Gassmann(1951)方程在微裂隙中充满油气;
步骤4,通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙,设定其纵横比为0.8,并通过Brown-Korringa(1975)方程实现固体替代;
步骤5,利用DEM模型向前面形成的混合物中加入基质孔隙,分为粒内孔、粒间孔、裂缝,得到干岩石的体积模量与剪切模量;
步骤6,通过Gassmann方程进行流体替代,计算饱和流体岩石的体积模量与剪切模量;
步骤7,利用计算的体积模量与剪切模量,结合计算的总密度和岩石速度之间的关系,得到泥页岩的纵波速度和横波速度,并与实测数据进行对比,判断是否吻合(允许误差范围10%),如果吻合,结束,如果不吻合,返回步骤5,调整基质孔隙中不同孔隙形态所占的体积分数,直至吻合(注意,实际测井曲线中往往只有纵波测井曲线,此时可以用纵波速度进行约束)。
图2为通过测井曲线计算,结合地化分析干酪根数据进行约束得到的干酪根体积分数分布特征,其与地化分析数据的相关系数达到0.86,证明了用来评价的干酪根分布特征的可靠性及准确性。
图3为依据测井曲线得到的各种矿物组分分布及干酪根计算数据,根据本文的模型预测得到的纵波速度与实测测井纵波速度的对比,其相关系数达到0.99,误差率在1%,证明了该模型的准确性。
图4与图5为根据计算的干酪根数据与基于本文模型得到的纵横波速度得到的含有机质泥页岩的弹性参数随干酪根体积分数的变化特征,这与前人的研究结果是一致的,侧面说明了本文模型建模思路的正确性。
图6为依据本模型得到的综合数据以及上面分析的干酪根对含有机质岩石弹性特征影响的分析,优选的页岩气井发育较好层段,从微电阻率扫描成像上看裂缝很发育,这与气测结果及实际压裂结果是一致的。 

Claims (5)

1.一种建立复杂孔隙泥页岩岩石物理模型方法,其特征在于,包括以下步骤:
步骤1,结合测井曲线,岩心数据及地质地化报告,得到组成泥页岩各组分的含量、弹性参数、孔隙度、有机质分布特征;
步骤2,根据步骤1得到的矿物组分的体积分数及弹性参数,通过Voigt-Reuss-Hill平均公式(1952)得到岩石骨架的体积模量和剪切模量; 
步骤3,将步骤1得到的孔隙度分为基质孔隙和干酪根孔隙,通过DEM(微分等效介质)模型得到干酪根与干酪根孔隙的混合物,设定干酪根孔隙的纵横比为0.01,并通过Gassmann(1951)方程在微裂隙中充满油气;
步骤4,通过DEM模型向岩石骨架中加入与干酪根与油气混合物所占体积相等的孔隙,设定其纵横比为0.8,并通过Brown-Korringa(1975)方程实现固体替代;
步骤5,利用DEM模型向前面形成的混合物中加入基质孔隙,分为粒内孔、粒间孔、裂缝,得到干岩石的体积模量与剪切模量;
步骤6,通过Gassmann(1951)方程进行流体替代,计算饱和流体岩石的体积模量与剪切模量;
步骤7,利用计算的体积模量与剪切模量,结合计算的总密度和岩石速度之间的关系,得到泥页岩的纵波速度和横波速度,并与实测数据进行对比,判断是否吻合,如果吻合,结束,如果不吻合,返回步骤5,调整基质孔隙中不同孔隙形态所占的体积分数,直至吻合。
2.根据权利要求1所述的方法,其特征在于,所述吻合的误差范围为10%。
3.根据权利要求1所述的方法,其特征在于,当实际测井曲线中只有纵波测井曲线时,用纵波速度进行约束。
4.根据权利要求1所述的方法,其特征在于,步骤3所述的孔隙度分为基质孔隙与干酪根孔隙,基质孔隙分为粒内孔、粒间孔、裂缝,分别用纵横比0.8、0.15、0.01来表征,干酪根孔隙为充填其中的微裂隙,设定纵横比为0.01,通过DEM模型实现干酪根与微裂隙的混合,并通过Gassmann实现微裂隙中充满油气。
5.根据权利要求1所述的方法,其特征在于,步骤4所述干酪根与油气混合物与基质矿物的混合,通过DEM模型加入干孔隙,设定纵横比为0.8,并通过Brown-Korringa方程来实现固体替代,加入干酪根与油气混合物。
CN201410153372.8A 2014-04-16 2014-04-16 一种建立复杂孔隙泥页岩岩石物理模型方法 Expired - Fee Related CN104007485B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410153372.8A CN104007485B (zh) 2014-04-16 2014-04-16 一种建立复杂孔隙泥页岩岩石物理模型方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410153372.8A CN104007485B (zh) 2014-04-16 2014-04-16 一种建立复杂孔隙泥页岩岩石物理模型方法

Publications (2)

Publication Number Publication Date
CN104007485A true CN104007485A (zh) 2014-08-27
CN104007485B CN104007485B (zh) 2016-01-20

Family

ID=51368228

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410153372.8A Expired - Fee Related CN104007485B (zh) 2014-04-16 2014-04-16 一种建立复杂孔隙泥页岩岩石物理模型方法

Country Status (1)

Country Link
CN (1) CN104007485B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316977A (zh) * 2014-11-12 2015-01-28 中国石油大学(华东) 一种页岩气储层的岩石模量计算方法
CN105203739A (zh) * 2015-10-29 2015-12-30 中国石油大学(北京) 一种基于岩石物理模型的页岩地层各向异性参数预测方法
CN105301642A (zh) * 2015-11-02 2016-02-03 中国石油天然气股份有限公司 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置
CN106285653A (zh) * 2015-06-05 2017-01-04 中国石油化工股份有限公司 确定页岩气储层的有机孔隙度的方法
CN106290105A (zh) * 2016-07-20 2017-01-04 中国石油大学(华东) 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法
CN108957580A (zh) * 2017-05-22 2018-12-07 中国石油化工股份有限公司 一种反演页岩地层孔隙结构的方法
CN109143397A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 碳酸盐岩储层缝洞充填物识别方法及系统
CN109655903A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 页岩层横波速度预测方法及系统
CN110189353A (zh) * 2019-06-10 2019-08-30 中国石油大学(华东) 一种泥页岩能谱矿物分布图校准方法及系统
CN110471129A (zh) * 2019-09-20 2019-11-19 西南石油大学 一种深层页岩高温高压下的各向异性岩石物理建模方法
CN110552690A (zh) * 2018-05-30 2019-12-10 中国石油化工股份有限公司 一种页岩储层脆性评价的方法
CN110836836A (zh) * 2018-08-17 2020-02-25 中国石油化工股份有限公司 干酪根弹性模量测试方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110292767A1 (en) * 2010-05-27 2011-12-01 Jianchun Dai Determining a quantity of a given material in a subterranean structure
CN103713327A (zh) * 2013-12-31 2014-04-09 孙赞东 基于测井及地震资料的页岩气储层最小闭合压力评价方法
CN103713320A (zh) * 2013-12-31 2014-04-09 孙赞东 一种富有机质泥页岩岩石物理模型的建立方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110292767A1 (en) * 2010-05-27 2011-12-01 Jianchun Dai Determining a quantity of a given material in a subterranean structure
CN103713327A (zh) * 2013-12-31 2014-04-09 孙赞东 基于测井及地震资料的页岩气储层最小闭合压力评价方法
CN103713320A (zh) * 2013-12-31 2014-04-09 孙赞东 一种富有机质泥页岩岩石物理模型的建立方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
施振飞: "江苏油田泥灰岩裂缝性储层测井解释方法研究", 《测井技术》 *
章海宁等: "适用复杂孔隙结构地层消除背景导电饱和度模型", 《测井技术》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316977A (zh) * 2014-11-12 2015-01-28 中国石油大学(华东) 一种页岩气储层的岩石模量计算方法
CN104316977B (zh) * 2014-11-12 2017-02-15 中国石油大学(华东) 一种页岩气储层的岩石模量计算方法
CN106285653A (zh) * 2015-06-05 2017-01-04 中国石油化工股份有限公司 确定页岩气储层的有机孔隙度的方法
CN105203739A (zh) * 2015-10-29 2015-12-30 中国石油大学(北京) 一种基于岩石物理模型的页岩地层各向异性参数预测方法
CN105301642A (zh) * 2015-11-02 2016-02-03 中国石油天然气股份有限公司 非均匀孔隙岩石及其固态有机质体积含量确定方法及装置
CN106290105A (zh) * 2016-07-20 2017-01-04 中国石油大学(华东) 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法
CN106290105B (zh) * 2016-07-20 2017-08-11 中国石油大学(华东) 一种碳酸盐岩储层溶蚀孔隙体积含量预测方法
CN108957580A (zh) * 2017-05-22 2018-12-07 中国石油化工股份有限公司 一种反演页岩地层孔隙结构的方法
CN109143397A (zh) * 2017-06-28 2019-01-04 中国石油化工股份有限公司 碳酸盐岩储层缝洞充填物识别方法及系统
CN109143397B (zh) * 2017-06-28 2020-05-19 中国石油化工股份有限公司 碳酸盐岩储层缝洞充填物识别方法及系统
CN109655903A (zh) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 页岩层横波速度预测方法及系统
CN110552690A (zh) * 2018-05-30 2019-12-10 中国石油化工股份有限公司 一种页岩储层脆性评价的方法
CN110836836A (zh) * 2018-08-17 2020-02-25 中国石油化工股份有限公司 干酪根弹性模量测试方法
CN110189353A (zh) * 2019-06-10 2019-08-30 中国石油大学(华东) 一种泥页岩能谱矿物分布图校准方法及系统
CN110189353B (zh) * 2019-06-10 2021-01-19 中国石油大学(华东) 一种泥页岩能谱矿物分布图校准方法及系统
CN110471129A (zh) * 2019-09-20 2019-11-19 西南石油大学 一种深层页岩高温高压下的各向异性岩石物理建模方法

Also Published As

Publication number Publication date
CN104007485B (zh) 2016-01-20

Similar Documents

Publication Publication Date Title
CN104007485A (zh) 一种建立复杂孔隙泥页岩岩石物理模型方法
CN102426390B (zh) 一种非均质泥砂岩储层储量确定方法
CN103713320B (zh) 一种富有机质泥页岩岩石物理模型的建立方法
CN104516017B (zh) 一种碳酸盐岩岩石物理参数地震反演方法
CN102445708B (zh) 三维等效富泥质砂岩速度预测模型
CN107203005B (zh) 一种定量化计算裂缝描述参数的方法
CN104007482B (zh) 一种基于各向异性有效场的泥页岩岩石物理模型方法
CN105840189A (zh) 一种页岩气岩石物理与测井方法
CN104806232B (zh) 一种确定碳酸盐岩储层孔隙度下限的方法
CN108399270A (zh) 一种确定页岩地层中各向异性泥质比例的方法
CN104181585A (zh) 一种地球物理勘探横波估算方法和系统
CN110824556A (zh) 一种非常规致密砂岩储层的岩石物理模型建立方法及应用
CN103713327B (zh) 基于测井及地震资料的页岩气储层最小闭合压力评价方法
CN103867198B (zh) 碳酸盐岩天然气层地层密度判别法
CN109509111A (zh) 探井地层压力的预测方法及系统
CN104570065B (zh) 一种利用地震波阻抗定量反演孔隙度的方法
CN102156297A (zh) 基于砂岩油藏叠后地震数据的流体替换方法
CN106202879A (zh) 二氧化碳驱油过程中纵横波速度动态估算方法
CN105317429A (zh) 一种确定页岩地层干酪根含量的方法
CN106096081A (zh) 缝洞型底水油藏的储量估算方法
CN101315428A (zh) 基于横波速度的高孔隙碎屑岩沉积地层的高分辨率勘探方法
CN108665545A (zh) 录井参数三维地质模型建立方法
CN108875115B (zh) 一种确定岩石强度的方法
CN110954942A (zh) 粘土含量识别的图版制作方法及系统
CN104198516B (zh) 应用核磁共振t2几何平均值预测横波速度的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: SUN ZANDONG

Effective date: 20150626

Owner name: SUN YONGYANG

Free format text: FORMER OWNER: SUN ZANDONG

Effective date: 20150626

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20150626

Address after: 102249, China University of Petroleum, 18, Xuefu Road, Beijing, Changping District (Beijing)

Applicant after: Sun Yongyang

Applicant after: Sun Zandong

Address before: Beijing city Changping District Road No. 18, China University of Petroleum (Beijing)

Applicant before: Sun Zandong

C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190314

Address after: 102249 18 Fu Xue Road, Changping District, Beijing

Patentee after: China University of Petroleum (Beijing)

Address before: 102249 Xuefu Road, Changping District, Beijing, China University of Petroleum (Beijing)

Co-patentee before: Sun Zandong

Patentee before: Sun Yongyang

TR01 Transfer of patent right
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: 20160120

Termination date: 20210416