CN111310339B - 一种等效时间修正的变储量系统典型曲线产生方法 - Google Patents
一种等效时间修正的变储量系统典型曲线产生方法 Download PDFInfo
- Publication number
- CN111310339B CN111310339B CN202010102422.5A CN202010102422A CN111310339B CN 111310339 B CN111310339 B CN 111310339B CN 202010102422 A CN202010102422 A CN 202010102422A CN 111310339 B CN111310339 B CN 111310339B
- Authority
- CN
- China
- Prior art keywords
- reserve
- pressure
- time
- typical curve
- stage
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000004519 manufacturing process Methods 0.000 claims abstract description 20
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 28
- 238000012937 correction Methods 0.000 claims description 16
- 239000012530 fluid Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 14
- 230000000694 effects Effects 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/08—Measuring diameters or related dimensions at the borehole
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Business, Economics & Management (AREA)
- Marine Sciences & Fisheries (AREA)
- Marketing (AREA)
- Agronomy & Crop Science (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Animal Husbandry (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种等效时间修正的变储量系统典型曲线产生方法,包括以下步骤:调用封闭边界油藏单井渗流模型;确定第k阶段开始时的地层压力;获得第k阶段的虚拟生产时间;修正渗流数学模型的生产时间,产生理论图版;将理论图版既典型曲线进行拟合,检验解释模型及解释参数的合理性。本发明提出的修正不同控制储量阶段等效时间的变储量系统典型曲线产生方法,沿用了单一固定储量系统的不稳定流动压力计算方式,产生多阶段控制储量变异的典型曲线,在变储量交界点上计算压力连续,展现和保持了阶段储量的动态特征,为利用现有的封闭边界油藏单井渗流模型分析变储量系统提供了简便方法。
Description
技术领域
本发明涉及油气井生产动态分析与储量评价技术领域,具体涉及一种等效时间修正的变储量系统典型曲线产生方法。
背景技术
传统的单井动态分析(例如现代产能递减分析)理论模型和解释方法中,均是假设单井的控制储量固定,将流压的下降视为两个组成部分:一是流体向井运移的流动压降,二是能量衰竭的地层压力下降;流动压降与瞬时产量相关,地层压力下降与累积产量相关。实际的油气藏生产系统中单井控制储量并非固定,将单井流动系统的控制储量视为可变量:一是有新储量贡献影响,油气藏的非均质性导致不同储集体之间存在一定的启动压差,当压差大于启动压差后,该储集体才开始有贡献,呈现出储量增加趋势;二是有井间干扰影响,在连通性强的系统中新井投产及周边邻井的产量变化,导致油气井控制储量的变化。
油气井的生产动态中控制储量的信息特征展现不突出,在动态分析中需要将原始的产量、压力转换为特征分析曲线,以突出控制储量及流动能力的影响,最常用的特征分析曲线为Blasingame诊断曲线和双对数压降诊断曲线,通过渗流数学模型产生理论图版既典型曲线对实测数据进行拟合,检验解释模型及解释参数的合理性。
对于含有变储量效应的实际生产动态,其Blasingame诊断曲线和双对数压降诊断曲线中不同于常规形态,表现出储量变异的显著影响。当储量增大时,Blasingame曲线中的-45°直线段右移,双对数压降曲线中的晚期45°压力导数直线段右移;当储量减少时,Blasingame曲线中的-45°直线段左移,双对数压降曲线中的晚期45°压力导数直线段左移。
直接开发新的变储量系统渗流数学模型的解析解,需要针对不同井型、不同边界、不同储量段数,分别建立渗流模型进行解析求解,这种方式既复杂又低效,无法利用现有的渗流模型成果,某些问题还难以求解;如何将现有的常规定产压降典型曲线改造为变储量系统的图版,将是一种简洁高效的计算途径。
发明内容
一种等效时间修正的变储量系统典型曲线产生方法,基于物质平衡原理,通过修正不同控制储量阶段的等效时间,保持沿用常规的单一固定储量系统不稳定流动压力计算方式,产生多阶段控制储量变异的典型曲线。
本发明目的在于为现有的封闭边界油藏单井渗流模型分析变储量系统提供简便方法。本发明的目的是通过以下技术方案来实现的:
一种等效时间修正的变储量系统典型曲线产生方法,包括如下步骤:
S1,根据压降模型计算出变储量系统在第k阶段储量开始时的地层压力Pr,k-1;
S2,假设孤立的k阶段储量系统从原始地层压力Pi开始,将按照虚拟的生产时间tp,k-1'在储量变化交界处产生的地层压力,等于所述步骤S1中的地层压力Pr,k-1,对比得出虚拟生产时间tp,k-1',获得相对真实开始时间tp,k-1的变储量影响的时间修正量;
S3,进入第k阶段后,计算变储量的时间修正量Δtk-1,生产时间t的典型曲线理论流压Pwf(t)变换为等效时间t+Δtk-1计算,在本阶段内变储量的时间修正量Δtk-1保持不变,进一步计算出典型曲线的压力,获得渗流数学模型的理论图版;
S4,通过所述步骤S3的渗流数学模型产生理论图版既典型曲线对实测数据进行拟合,检验解释模型及解释参数的合理性。
进一步的,所述步骤S1中的地层压力Pr,k-1,由前期k-1个储量系统的总压降表示:
其中,Pi为原始地层压力,q为定产量,Nk为第k阶段的单井控制储量,tp,k-1为第k阶段的开始时间,tp,k为第k阶段的结束时间,Ct为综合压缩系数,取tp,0=0。
进一步的,所述步骤S2中按照虚拟的生产时间tp,k-1'在储量变化交界处产生的地层压力由下式表示:
同步骤S1中的地层压力对比可得虚拟的生产时间tp,k-1'为:
相对真实开始时间tp,k-1的变储量影响的时间修正量为:
Δtk-1=tp,k-1'-tp,k-1
其中,k=1时Δtk-1=0。
进一步的,所述步骤S3中的典型曲线的压力按具体井型用不稳定试井的完整渗流数学模型计算,如下式所示:
Pwf,k=Ptrans(q,Pi,Rw,h,φ,Ct,μ,B,Re,k,...,t+Δtk-1)
其中,Pwf,k为第k阶段的井底流压,Ptrans()为不稳定渗流的压力解式;q为产量,Rw为井筒半径,h为储层厚度,φ为孔隙度,Ct为综合压缩系数,μ为流体粘度,B为流体体积系数;Re,k为第k阶段的油藏控制半径,对应第k阶段的控制储量;t为典型曲线的计算时间。
本发明的有益效果:本发明提出的修正不同控制储量阶段等效时间的变储量系统典型曲线产生方法,沿用了单一固定储量系统的不稳定流动压力计算方式,产生多阶段控制储量变异的典型曲线,在变储量交界点上计算压力连续,展现和保持了阶段储量的动态特征;本方法为利用现有的封闭边界油藏单井渗流模型分析变储量系统提供了简便方法。
附图说明
图1是本发明的方法流程图。
图2是本发明的多阶段变控制储量的地层压力下降的等效时间分析图。
图3是本发明的多阶段变控制储量等效时间的修正量变化图。
图4是本发明应用分析实例的xc10井压力指数分析曲线图。
图5是本发明应用分析实例的xc10井阶段变储量Blasingame诊断拟合曲线图。
图6是本发明应用分析实例的xc10井阶段变储量双对数压力拟合曲线图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图说明本发明的具体实施方式。
本实施例中,如图1所示,一种等效时间修正的变储量系统典型曲线产生方法如下:
第一步,根据压降模型计算出储量系统在第k阶段储量开始时的地层压力Pr,k-1。其中,多阶段变控制储量等效时间分析如图2所示,系统的原始地层压力为Pi,第k阶段的单井控制储量为Nk,设第k阶段的结束时间为tp,k,开始时间为tp,k-1。在恒定产量条件下当流动进入拟稳定流动期后,地层压力随生产时间线性下降,控制储量控制了地层压力的下降速率。对定产量q、恒定的综合压缩系数Ct的变储量系统,在第k阶段开始时的地层压力由前期k-1个储量系统的总压降表示:
其中,Pi为原始地层压力,取tp,0=0。
第二步,假设孤立的k阶段储量系统从原始地层压力Pi开始,按照虚拟的生产时间tp,k-1',在储量变化的交界处产生相同的地层压力:
同第一步中的计算出的地层压力进行对比可得,虚拟的生产时间tp,k-1'为:
相对真实开始时间tp,k-1的变储量影响的时间修正量为:
Δtk-1=tp,k-1'-tp,k-1
其中,k=1时Δtk-1=0。
第三步,进入第k阶段后,计算变储量的时间修正量Δtk-1,生产时间t的典型曲线理论流压Pwf(t)变换为等效时间t+Δtk-1计算,在本阶段内变储量的时间修正量Δtk-1保持不变,典型曲线的压力按具体的井型如直井、水平井等用不稳定试井的完整渗流数学模型计算:
Pwf,k=Ptrans(q,Pi,Rw,h,φ,Ct,μ,B,Re,k,...,t+Δtk-1)
其中,Pwf,k为第k阶段的井底流压,Ptrans()为不稳定渗流的压力解式;q为产量,Rw为井筒半径,h为储层厚度,φ为孔隙度,Ct为综合压缩系数,μ为流体粘度,B为流体体积系数;Re,k为第k阶段的油藏控制半径,对应第k阶段的控制储量;t为典型曲线的计算时间。
第四步,根据渗流数学模型产生的理论图版即典型曲线对实测数据进行拟合,检验解释模型及解释参数的合理性。
本实施例中,变储量修正时间可为正或负,如图3所示,以2阶段控制储量的系统为例,当控制储量增长时,修正时间Δt>0,等效时间延长;当控制储量下降时,修正时间Δt<0,等效时间缩短。
应用分析实例:
某储气库气井xc10井在第1阶段采气期表现出变控制储量现象,如图4所示压力指数曲线呈现2段斜率,反映出初期控制储量小、后期控制储量大,阶段控制储量分析结果见表1。变控制储量的定产压降的Blasingame图版拟合曲线见图5、双对数压降图版拟合曲线图6。
表1
本发明提出的修正不同控制储量阶段等效时间的变储量系统典型曲线产生方法,沿用了单一固定储量系统的不稳定流动压力计算方式,产生多阶段控制储量变异的典型曲线,在变储量交界点上计算压力连续,展现和保持了阶段储量的动态特征;本方法为利用现有的封闭边界油藏单井渗流模型分析变储量系统提供了简便方法。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种等效时间修正的变储量系统典型曲线产生方法,其特征在于,包括以下步骤:
S1,根据压降模型计算出变储量系统在第k阶段储量开始时的地层压力Pr,k-1;
S2,假设孤立的k阶段储量系统从原始地层压力Pi开始,按照虚拟的生产时间tp,k-1'在储量变化交界处产生的地层压力,等于所述步骤S1中的地层压力Pr,k-1,对比得出虚拟生产时间tp,k-1',获得相对真实开始时间tp,k-1的变储量影响的时间修正量;
S3,进入第k阶段后,计算变储量的时间修正量Δtk-1,生产时间t的典型曲线理论流压Pwf(t)变换为等效时间t+Δtk-1计算,在本阶段内变储量的时间修正量Δtk-1保持不变,进一步计算出典型曲线的压力,获得渗流数学模型的理论图版;
S4,对所述步骤S3的渗流数学模型产生理论图版既典型曲线对实测数据进行拟合,检验解释模型及解释参数的合理性;
所述步骤S2中在储量变化交界处产生的地层压力由下式表示:
其中,tp,k-1'为虚拟的生产时间,同步骤S2中的地层压力对比可得虚拟的生产时间tp,k-1'为:
相对真实开始时间tp,k-1的变储量影响的时间修正量为:
Δtk-1=tp,k-1'-tp,k-1
其中k=1时Δtk-1=0。
3.根据权利要求1所述的一种等效时间修正的变储量系统典型曲线产生方法,其特征在于,所述步骤S3中的典型曲线的压力按具体井型用不稳定试井的完整渗流数学模型计算,如下式所示:
Pwf,k=Ptrans(q,Pi,Rw,h,φ,Ct,μ,B,Re,k,...,t+Δtk-1)
其中,Pwf,k为第k阶段的井底流压,Ptrans()为不稳定渗流的压力解式;q为产量,Rw为井筒半径,h为储层厚度,φ为孔隙度,Ct为综合压缩系数,μ为流体粘度,B为流体体积系数;Re,k为第k阶段的油藏控制半径,对应第k阶段的控制储量;t为典型曲线的计算时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010102422.5A CN111310339B (zh) | 2020-02-19 | 2020-02-19 | 一种等效时间修正的变储量系统典型曲线产生方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010102422.5A CN111310339B (zh) | 2020-02-19 | 2020-02-19 | 一种等效时间修正的变储量系统典型曲线产生方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111310339A CN111310339A (zh) | 2020-06-19 |
CN111310339B true CN111310339B (zh) | 2020-12-04 |
Family
ID=71158380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010102422.5A Expired - Fee Related CN111310339B (zh) | 2020-02-19 | 2020-02-19 | 一种等效时间修正的变储量系统典型曲线产生方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111310339B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6842700B2 (en) * | 2002-05-31 | 2005-01-11 | Schlumberger Technology Corporation | Method and apparatus for effective well and reservoir evaluation without the need for well pressure history |
CN103995961A (zh) * | 2014-05-06 | 2014-08-20 | 中国石油天然气股份有限公司 | 一种测试井的生产数据分析方法 |
CN104500030A (zh) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | 超高压气藏生产动态异常数据诊断及修正方法 |
CN104632187A (zh) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | 一种确定产水煤层气井动态储量的方法 |
CN105913332A (zh) * | 2016-04-22 | 2016-08-31 | 中国石油天然气股份有限公司 | 一种确定致密气藏开发指标的方法及装置 |
CN108804819A (zh) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | 一种低渗气藏动态储量评价方法 |
CN109931056A (zh) * | 2019-02-21 | 2019-06-25 | 西南石油大学 | 一种变控制储量气井动态分析方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413030B (zh) * | 2013-07-24 | 2017-07-07 | 中国石油天然气股份有限公司 | 一种缝洞型碳酸盐岩气藏动态分析方法及系统 |
CN104141490B (zh) * | 2014-07-17 | 2016-08-31 | 中国石油天然气股份有限公司 | 有水气藏单井水侵情况判断及气井产量控制方法和装置 |
CN105260543B (zh) * | 2015-10-19 | 2018-06-01 | 中国石油天然气股份有限公司 | 基于双孔模型的多重介质油气流动模拟方法及装置 |
US20180210977A1 (en) * | 2017-01-23 | 2018-07-26 | Saudi Arabian Oil Company | Systems and Methods for Automatic Reservoir Simulation Model Performance Tuning |
US10439433B2 (en) * | 2017-03-01 | 2019-10-08 | University Of Central Florida Research Foundation, Inc. | Adaptive power grid restoration |
-
2020
- 2020-02-19 CN CN202010102422.5A patent/CN111310339B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6842700B2 (en) * | 2002-05-31 | 2005-01-11 | Schlumberger Technology Corporation | Method and apparatus for effective well and reservoir evaluation without the need for well pressure history |
CN104632187A (zh) * | 2013-11-14 | 2015-05-20 | 中国石油化工股份有限公司 | 一种确定产水煤层气井动态储量的方法 |
CN103995961A (zh) * | 2014-05-06 | 2014-08-20 | 中国石油天然气股份有限公司 | 一种测试井的生产数据分析方法 |
CN104500030A (zh) * | 2014-11-21 | 2015-04-08 | 中国石油天然气股份有限公司 | 超高压气藏生产动态异常数据诊断及修正方法 |
CN105913332A (zh) * | 2016-04-22 | 2016-08-31 | 中国石油天然气股份有限公司 | 一种确定致密气藏开发指标的方法及装置 |
CN108804819A (zh) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | 一种低渗气藏动态储量评价方法 |
CN109931056A (zh) * | 2019-02-21 | 2019-06-25 | 西南石油大学 | 一种变控制储量气井动态分析方法 |
Non-Patent Citations (2)
Title |
---|
"致密低渗气藏气井动态分析方法";陈伟 等;《天然气勘探与开发》;20091025;参见正文第3章 * |
"页岩气藏压裂水平井Blasingame产量递减分析方法建立与应用";姜瑞忠 等;《石油学报》;20191215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111310339A (zh) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108612525B (zh) | 一种气藏动态储量计算方法 | |
CN109281663B (zh) | 一种确定煤层气储层渗透率的方法 | |
CN107506948B (zh) | 一种基于动态泄流体积的页岩油气综合产量分析方法 | |
CN105569646B (zh) | 一种油气井技术可采储量预测方法 | |
CN108804819A (zh) | 一种低渗气藏动态储量评价方法 | |
CN108756868B (zh) | 一种注气开发油藏开发中后期可采储量的评估方法 | |
CN104832156A (zh) | 一种预估气井产量的方法 | |
CN110878690B (zh) | 基于碳同位素演化确定页岩气井生产递减曲线的方法 | |
CN105631184A (zh) | 一种油气井产量递减分析方法和系统 | |
CN111401595B (zh) | 一种低粘度原油油藏水驱波及系数预测方法 | |
CN109931056B (zh) | 一种变控制储量气井动态分析方法 | |
CN114997080B (zh) | 一种缝洞型气藏未稳定产能测试资料解释方法 | |
CN111502652A (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN115293063A (zh) | 一种气水两相渗流预测模型及其建立方法和应用 | |
CN111310338B (zh) | 一种原始地层压力修正的变储量系统典型曲线产生方法 | |
CN107130955B (zh) | 井底流压的确定方法及储集体天然能量的确定方法 | |
CN111310339B (zh) | 一种等效时间修正的变储量系统典型曲线产生方法 | |
CN113445988B (zh) | 一种低渗碳酸盐岩气藏气井产能评价方法 | |
CN108520143B (zh) | 一种注气开发油藏气油比上升率表征方法 | |
CN116127874B (zh) | 一种凝析井地层压力及单井控制储量联合一体的评价方法 | |
CN108266184A (zh) | 一种基于生产数据的气井单井控制储量估算方法 | |
CN110895637A (zh) | 一种产水气井单井控制储量的方法 | |
CN111305825B (zh) | 一种变控制储量的气井压力动态模拟方法 | |
CN115345090A (zh) | 一种欠饱和煤层气储层压降漏斗动态传播的计算方法 | |
CN112228053A (zh) | 致密气藏间歇井开井时间优化方法 |
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 |
Granted publication date: 20201204 Termination date: 20220219 |
|
CF01 | Termination of patent right due to non-payment of annual fee |