CN105822264B - 天然气水合物层钻井井筒水合物动态分解位置检测方法 - Google Patents

天然气水合物层钻井井筒水合物动态分解位置检测方法 Download PDF

Info

Publication number
CN105822264B
CN105822264B CN201610162158.8A CN201610162158A CN105822264B CN 105822264 B CN105822264 B CN 105822264B CN 201610162158 A CN201610162158 A CN 201610162158A CN 105822264 B CN105822264 B CN 105822264B
Authority
CN
China
Prior art keywords
pressure
hydrate
gas
decomposition
drilling
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
Application number
CN201610162158.8A
Other languages
English (en)
Other versions
CN105822264A (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.)
Southwest Petroleum University
China National Offshore Oil Corp CNOOC
Original Assignee
Southwest Petroleum University
China National Offshore Oil Corp CNOOC
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 Southwest Petroleum University, China National Offshore Oil Corp CNOOC filed Critical Southwest Petroleum University
Priority to CN201610162158.8A priority Critical patent/CN105822264B/zh
Publication of CN105822264A publication Critical patent/CN105822264A/zh
Application granted granted Critical
Publication of CN105822264B publication Critical patent/CN105822264B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/01Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells specially adapted for obtaining from underwater installations
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B41/00Equipment or details not covered by groups E21B15/00 - E21B40/00
    • E21B41/0099Equipment or details not covered by groups E21B15/00 - E21B40/00 specially adapted for drilling for or production of natural hydrate or clathrate gas reservoirs; Drilling through or monitoring of formations containing gas hydrates or clathrates
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/001Survey of boreholes or wells for underwater installation
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/04Measuring depth or liquid level
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/06Measuring temperature or pressure
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/10Locating fluid leaks, intrusions or movements

Abstract

本发明公开了一种天然气水合物层钻井井筒水合物动态分解位置检测方法,特别是一种应用于海洋天然气水合物的勘探钻井或者开采钻井领域的天然气水合物层钻井井筒水合物动态分解位置检测方法。本发明提供一种可以实时采集数据,检测准确度高,并且可以实时检测处动态变化的水合物分解位置的天然气水合物层钻井井筒水合物动态分解位置检测方法。包括以下几个步骤:A、按照钻井井深位置范围划分节点;B、将参数代入多相流模型计算混合泥浆密度分布和环空压耗;C、计算立管压力;D、比较计算立管压力与监测立管压力。本申请通过地面监测立压的变化,计算并实时预测水合物动态分解位置,有效降低海洋天然气水合物层钻井风险、减少损失。

Description

天然气水合物层钻井井筒水合物动态分解位置检测方法
技术领域
[0001]本发明涉及一种天然气水合物层钻井井筒水合物动态分解位置检测方法,特别是 一种应用于海洋天然气水合物的勘探钻井或者开采钻井领域的天然气水合物层钻井井筒 水合物动态分解位置检测方法。
背景技术
[0002] 天然气水合物又称“可燃冰”,是一种高密度、高热值的非常规能源(lm3天然气水 合物可以释放出164m3甲烷气和0 • 8m37j〇,其主要以层的形式存在于海底岩层之间,且资源 量极为丰富。海洋天然气水合物层钻井过程中,水合物颗粒随钻井液上返到钻井平台,此过 程中的管流不同于常规油气钻井过程:天然气水合物在上返过程中,压力不断降低、温度不 断升高,至一定位置后,天然气水合物发生分解,使井筒压力变化明显,进一步造成井控安 全的一些列问题,对海上钻井安全产生巨大风险。
[0003] 目前,国内外对海洋天然气水合物层钻井中天然气水合物在井筒上返过程中的分 解位置有一定研究,能在一定程度上实现海洋天然气水合物层钻井中水合物分解位置预 测。现有方法中有学者通过天然气水合物相平衡曲线,结合环空井筒中温度、压力变化,预 测天然气水合物分解位置变化,但是现有天然气水合物相平衡曲线主要应用为稳态条件 下,并且水合物分解位置是理论计算结果,不能完全反映实际海洋天然气水合物层钻井中 水合物动态分解实时位置的变化。专利CN102865066A公开了一种含天然气水合物相变的深 水井筒多相流动实验装置及方法,可以实现含深水井筒内天然气水合物分解对井筒规律影 响的模拟,但是其并未提出一种实时有效的天然气水合物分解位置预测方法。专利 CN104198674A公开了一种管输天然气水合物形成在线预警装置及方法,可以计算出监测点 处水相中水合物抑制剂和盐分浓度,结合采集到的气相压力、组分和温度参数实现预警,但 是其需要复杂的数据在线采集模块,在海洋钻井中无法有效应用。因此现有技术中还没有 一种可以实时采集数据,操作简单,检测准确度高,并且可以实时检测处动态变化的水合物 分解位置的天然气水合物层钻井井筒水合物动态分解位置检测方法。
发明内容
[0004] 本发明所要解决的技术问题是提供一种可以实时采集数据,操作简单,检测准确 度高,并且可以实时检测处动态变化的水合物分解位置的天然气水合物层钻井井筒水合物 动态分解位置检测方法。
[0005] 为解决上述技术问题本发明采用的天然气水合物层钻井井筒水合物动态分解位 置检测方法,包括以下几个步骤:
[0006] A、按照钻井井深位置范围划分节点,根据钻井工况和钻井返出参数计算出钻井的 气相流量、液相流量和固相流量,并选择其中一个节点作为下一步骤的计算节点。
[0007] B、将气相流量、液相流量和固相流量以及所选计算节点的井深代入多相流模型中 得出混合泥浆密度分布和环空压耗。 y_8] C、将前述步骤的混合泥衆密度分布轉空压耗代入井筒压力模型中得出计算立 官压力; _9] D、将计算立管压力与難立籠力肺隨,如果计歡髓力与监测立管压力 相等则天然气水合物分解位置为所选计算节点的井深位置,贼计算立管压力与监测立管 压力不相_选择下-个_节点㈣计算节点并難频晅步狐 、
[0010]进一步的是,将D步骤中得出的不同时间的天然气水合物分解位置拟合成曲线,通 过曲线预测下一时刻天然气水合物分解位置。
[0011]本发明的有益效果是:本申请通过地面监测立压的变化,结合井筒返出参数,根据 井筒压力模型、井筒多相流模型,计算并实时删水合物动态分解位置,其需要采集的数据 少,并且能够及时为水合物分解所造成的井控安全等一系列问题的处理提供指导,有效降 低海洋天然气水合物层钻井风险、减少损失。
附图说明
[0012]图1是本发明中的海洋天然气水合物层钻井中井筒压力关系示意图;
[0013]图2是本发明中的根据气、液相表观流速判定多相流流型的模板示意图;
[00M]图3是本发明中的通过立压实时预测水合物分解位置的方法示意图;
[0015]图4是本发明中的立压变化和水合物分解位置趋势示例图;
具体实施方式
[0016]下面结合附图对本发明作进一步说明。
[0017]如图2所示,本发明的天然气水合物层钻井井筒水合物动态分解位置检测方法,包 括以下几个步骤:
[0018] A、按照钻井井深位置范围划分节点,根据钻井工况和钻井返出参数计算出钻井的 气相流量、液相流量和固相流量,并选择其中一个节点作为下一步骤的计算节点;
[0019] B、将气相流量、液相流量和固相流量以及所选计算节点的井深代入多相流模型中 得出混合泥浆密度分布和环空压耗;
[0020] C、将前述步骤的混合泥浆密度分布和环空压耗代入井筒压力模型中得出计算立 管压力;
[0021] D、将计算立管压力与监测立管压力进行比较,如果计算立管压力与监测立管压力 相等则天然气水合物分解位置为所选计算节点的井深位置,如果计算立管压力与监测立管 压力不相等则选择下一个相邻节点作为计算节点并重复步骤B至步骤D。
[0022]为了使计算的天然气水合物分解位置与钻井的深度位置对应,将钻井深度划分成 若干个长度区间,每个长度区间对应一个节点,具体实施时可以给划分好的节点编号,例如 井深最深处的节点编为1号节点,在1号"P点上方并且紧邻1号节点的为2号节点,并依次将 剩下的节点编号。在计算时可以选择最深处的1号节点作为初始节点,进行第一次计算,如 果1号节点不是天然气水合物的分解位置,则可以选择与其相邻的2号节点进行计算,依次 类推直至算出准确的天然气水合物的分解位置。根据钻井工况和钻井返出参数计算出钻井 的气相流量、液相流量和固相流量可采用现有的计算方法进行。
[0023]海洋天然气水合物层钻井中,水合物发生分解时,井筒中钻井液的物理性质会发 生相应变化,如钻井液密度、流量等,从而影响到钻柱内及环空井筒压力的平衡状态,进而 使地面上监测到的立管压力等参数发生相应变化。其中,立压指作业过程中立管上的压力 值,反映的是钻井液在钻柱、钻头水眼及环空中的压力损失,在数值上近似于等于整个循环 系统压力损失,也近似等于循环栗压,可以通过监测得到。
[0024] 在常规正常钻进情况下,井筒内的流体循环满足质量守恒定律,即钻井液的总量 不会发生变化;而海洋天然气水合物层钻井过程中,天然气水合物随钻井液在环空井筒中 上返至一定位置会发生分解,分解产生天然气会使环空井筒中压力发生变化,破坏原有的 压力平衡,引起压力波动,从而表征到立压监测的变化上。如图1,海洋天然气水合物层钻井 过程中,井筒压力关系可以表示为:
[0025] Pd+Pmd = A Pdf + A Paf + A Pdb+Pma
[0026] 式中:pd为立管压力,MPa;pmd、pma分别为钻柱内、环空中液柱压力,单位分别为: MPa; A pdf、A paf、A pdb分别为钻柱内、环空中、钻头压耗,单位均为MPa。
[0027]天然气水合物在环空中随钻井液上返至一定位置发生分解后,环空中液柱压力和 摩阻压耗可以分别表示为:
[0028] Pma — Pmixghd+Pmg (h—hd) (2)
[0029] A paf = A pafl+ A paf2 (3)
[0030]式中:Pm为钻井液与水合物石肩混合密、度,单位为g/cm3;Pmix为水合物分解产生天 然气与钻井液和剩余水合物岩肩混合密度,单位为g/cm3; h、hd分别为井深、天然气水合物分 解位置,单位无km; A pafl、A paf2分别为天然气水合物分解位置上部、下部的环空压耗,单位 为 MPa。
[0031]当天然气水合物发生分解后,钻柱内液柱压力Pmd、钻井液在钻杆中压耗A pdf、钻 井液钻头压耗A pdb基本不发生变化,而环空中液柱压力pma会由于水合物分解位置上部混 合钻井液密度Pmix的降低而减小,且钻井液在环空中的压耗会由于井筒流动特征的变化而 发生变化,根据井筒压力关系,反映到监测立压Pd的变化上。
[0032]根据监测到的立压变化计算水合物分解位置时,根据钻井工况和返出参数监测, 得到返出气、液、固相流量,作为初始条件,并结合多相流模型如图2计算得到某一节点内混 合泥楽密度分布和环空压耗;进而根据井筒压力关系计算立压,并与监测立压比较是否相 等;如不相等则循环迭代计算下一节点内混合泥浆密度分布和环空压耗,并计算立压再次 比较;以此循环至计算所得立压与监测立压相等,则说明天然气水合物分解位置即为此计 算节点处。同时,随着分解的进行,立压将会不断变化,采用上述方法不断迭代,实时计算并 预测出水合物动态分解位置的变化。
[0033]本发明的通过立压实时预测水合物分解位置的方法如图3所示。在海洋天然气水 合物层钻井过程中,根据钻井工况和返出参数监测,得到返出气、液、固相流量,通过多相流 模型计算第一个节点内混合泥浆密度分布和环空压耗,并根据井筒压力关系计算立压,与 井口监测立压相比较;如果不相等则循环迭代下一节点计算,直至计算立压与井口监测立 压相等则说明水合物分解位置即为节点位置。通过系统不断实时迭代计算从而预测水合物 分解位置。 '
[0034]首先假设天然气水合物不分解,通过井筒压力关系计算得到立压:如果计算得到 立压与实测立压相一致,说明此时天然气水合物在井筒中未发生分解或还未钻遇天然气水 合物层;如果计算得到立压比实测立压大,说明此时天然气水合物在井筒中已经发生分解, 根据井口返出气、液、固相流量,结合多相流模型和井筒压力关系迭代计算得到天然气水合 物分解位置;如果监测的立压比计算得到立压小且不发生变化,说明天然气水合物动态分 解位置己经达到稳定状态;如果监测的立压随时间仍在不断减小,则说明天然气水合物分 解位置在不断下移,此时根据监测立压并结合返出参数实时计算,对天然气水合物动态分 解位置实时预测,直到监测的立压降低至某一稳定状态。
[0035]立压监测数据可以用来计算并得到天然气水合物分解的实时位置,从而通过立压 的变化来表征海洋天然气水合物层钻井水合物分解位置的变化,为施工人员进一步采取措 施提供实时数据支撑。
[0036]将D步骤中得出的不同时间的天然气水合物分解位置拟合成曲线,通过曲线预测 下一时刻天然气水合物分解位置。,随着分解的进行,立压将会不断变化,采用上述方法不 断迭代,实时计算并预测出水合物动态分解位置的变化。
[0037]结合图4作进一步说明,该图表示的是立压随时间的变化以及实时反映出的天然 气水合物分解位置随时间变化情况,其中横坐标是时间,纵坐标是立压和天然气水合物分 解位置。A、B、C、D、E分别代表五处不同的钻井时间,时间A〜B处,立压未发生变化,且与假设 天然气水合物不分解时计算得到的立压Pd〇相等,表明天然气水合物未发生分解;时间B〜C 处,立压降低较快,说明天然气水合物刚开始分解,环空井筒压力降低,其分解位置迅速向 下移动;时间C〜D处,立压降低减慢,说明随着天然气水合物分解位置的下移,井筒温度的 降低占主要因素,其下移速度变缓;时间D〜E处,立压保持稳定,说明环空井筒压力降低与 井筒温度降低效果相当,天然气水合物持续分解,但其分解位置基本稳定。天然气水合物分 解位置的变化过程中,应及时结合井筒压力变化采取必要措施,以保证海洋天然气水合物 钻井的安全进行。 ~

Claims (2)

1.天然气水合物层钻井井筒水合物动态分解位置检测方法,其特征在于:包括以下几 个步骤: A、 按照钻井井深位置范围划分节点,根据钻井工况和钻井返出参数计算出钻井的气相 流量、液相流量和固相流量,并选择其中一个节点作为下一步骤的计算节点; B、 将气相流量、液相流量和固相流量以及所选计算节点的井深代入多相流模型中得出 混合泥浆密度分布和环空压耗; C、 将前述步骤的混合泥浆密度分布和环空压耗代入井筒压力模型中得出计算立管压 力; 4 D、比较计算立管压力与监测立管压力是否相等,如果不相等则循环迭代计算下一相邻 节点=混合泥浆密度分布和环空压耗,并按照步骤(^计算立管压力再次比较;以此循环直至 计•胃A管压力与井口监测立管压力相等则说明天然气水合物分解位置即为此计算节点处; 同时•,随着分解的进行,立管压力将会不断变化,采用上述方法不断迭代,实时计算并预测 出水合物动态分解位置的变化。 ,2.¾权利要求i所述的天然气水合物层钻井井筒水合物动态分解位置检测方法,其特 于:将D步骤中得出的不同时间的天然气水合物分解位置拟合成曲线,通过曲线预测下 一时刻天然气水合物分解位置。
CN201610162158.8A 2016-03-18 2016-03-18 天然气水合物层钻井井筒水合物动态分解位置检测方法 Expired - Fee Related CN105822264B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610162158.8A CN105822264B (zh) 2016-03-18 2016-03-18 天然气水合物层钻井井筒水合物动态分解位置检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610162158.8A CN105822264B (zh) 2016-03-18 2016-03-18 天然气水合物层钻井井筒水合物动态分解位置检测方法

Publications (2)

Publication Number Publication Date
CN105822264A CN105822264A (zh) 2016-08-03
CN105822264B true CN105822264B (zh) 2018-08-17

Family

ID=56524201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610162158.8A Expired - Fee Related CN105822264B (zh) 2016-03-18 2016-03-18 天然气水合物层钻井井筒水合物动态分解位置检测方法

Country Status (1)

Country Link
CN (1) CN105822264B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108804722B (zh) * 2017-04-26 2021-07-27 中国石油化工股份有限公司 一种用于钻井仿真的参数计算方法及装置
CN111827909B (zh) * 2020-08-04 2021-06-01 中国石油大学(华东) 一种海域天然气水合物开路循环钻井中井筒压力的主动控制方法及控制装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006096779A (ja) * 2004-09-28 2006-04-13 National Institute Of Advanced Industrial & Technology 窒素によるメタンハイドレートの分解方法及び分解装置
CN101709639A (zh) * 2009-11-20 2010-05-19 中国石油大学(华东) 模拟深水油气开采的井筒多相流动装置
CN102865066A (zh) * 2012-10-16 2013-01-09 中国石油大学(华东) 含天然气水合物相变的深水井筒多相流动实验装置及方法
CN102943620A (zh) * 2012-08-27 2013-02-27 中国石油大学(华东) 基于钻井环空井筒多相流动计算的控压钻井方法
JP5229703B2 (ja) * 2008-12-12 2013-07-03 独立行政法人産業技術総合研究所 ガスハイドレートの生成・分解状況の検知方法
CN104198674A (zh) * 2014-08-29 2014-12-10 西南石油大学 管输天然气水合物形成在线预警装置及方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006096779A (ja) * 2004-09-28 2006-04-13 National Institute Of Advanced Industrial & Technology 窒素によるメタンハイドレートの分解方法及び分解装置
JP5229703B2 (ja) * 2008-12-12 2013-07-03 独立行政法人産業技術総合研究所 ガスハイドレートの生成・分解状況の検知方法
CN101709639A (zh) * 2009-11-20 2010-05-19 中国石油大学(华东) 模拟深水油气开采的井筒多相流动装置
CN102943620A (zh) * 2012-08-27 2013-02-27 中国石油大学(华东) 基于钻井环空井筒多相流动计算的控压钻井方法
CN102865066A (zh) * 2012-10-16 2013-01-09 中国石油大学(华东) 含天然气水合物相变的深水井筒多相流动实验装置及方法
CN104198674A (zh) * 2014-08-29 2014-12-10 西南石油大学 管输天然气水合物形成在线预警装置及方法

Also Published As

Publication number Publication date
CN105822264A (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
US9790743B2 (en) Natural gas hydrate formation drilling simulation device
US9103201B2 (en) Method and system for configuring crude oil displacement system
US20120330553A1 (en) Method for generating a general enhanced oil recovery and waterflood forecasting model
CN102943620A (zh) 基于钻井环空井筒多相流动计算的控压钻井方法
EP3500725B1 (en) Fluid production network leak detection
CN105822264B (zh) 天然气水合物层钻井井筒水合物动态分解位置检测方法
Zhang et al. Real-time estimation and management of hydrate plugging risk during deepwater gas well testing
Lin et al. Calculation analysis of pressure wave velocity in gas and drilling mud two-phase fluid in annulus during drilling operations
Sleiti et al. Early gas kick detection in vertical wells via transient multiphase flow modelling: A review
Zhou Research on heat transfer in geothermal wellbore and surroundings
Mao et al. Dynamical Simulation of High-Pressure Gas Kick in Ultra-Deepwater Riserless Drilling
Li et al. Reservoir characterization during underbalanced drilling of horizontal wells based on real-time data monitoring
Li et al. Model for calculating the gas volume fraction of a gas-cut wellbore through natural gamma-ray logging
He et al. Real-Time Interpretation Model of Reservoir Characteristics While Underbalanced Drilling Based on UKF
Biru Analysis of well testing, temperature and pressure in high-temperature wells of Aluto-Langano, Ethiopia
CN112627733A (zh) 深水控压钻井水力参数实时优化方法及设备
Kaya et al. Investigation of pressure transient analysis methods for single-phase and CO2-rich geothermal reservoirs
Yang et al. Development of a Dual-Measured-Points Early Gas Kick Detection Method based on the pressure responses of two PWD tools
Hill Research in multiphase flow developed into new global discipline of flow assurance-an oil industry perspective
Leng et al. Rational Allocation of Multilayer Gas Wells with Producing Water in Tight Gas Reservoirs
Gao et al. A Comprehensive Model for Simulating Supercritical-Water Flow in a Vertical Heavy-Oil Well
EP2397649A1 (en) Method and system for determining relative mass fluxes
Gao et al. The Influence of Dissolution of Hydrate in Blocked BOP on Gas Volume Fraction in Deepwater Drilling
Se et al. Novel Methodology to Assess Miscible Gas Flood Performance Combining Time-Lapse Logging and Simulation Model Output
Tareke Evaluation of wellbore parameters for well HE-4, Hellisheidi, SW-Iceland

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180817

Termination date: 20190318

CF01 Termination of patent right due to non-payment of annual fee