CN104778329B - 一种分析co2/ch4在褐煤中吸附机理的方法 - Google Patents
一种分析co2/ch4在褐煤中吸附机理的方法 Download PDFInfo
- Publication number
- CN104778329B CN104778329B CN201510202355.3A CN201510202355A CN104778329B CN 104778329 B CN104778329 B CN 104778329B CN 201510202355 A CN201510202355 A CN 201510202355A CN 104778329 B CN104778329 B CN 104778329B
- Authority
- CN
- China
- Prior art keywords
- lignite
- adsorption
- analysis
- absorption
- simulation
- 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
Landscapes
- Solid Fuels And Fuel-Associated Substances (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
Abstract
结合密度泛函理论(DFT)、分子动力学(MD)和蒙特卡洛模拟(GCMC)三种理论计算方法,在分子模拟水平上研究了CO2/CH4气体在一种新型褐煤架构中的吸附实例,归纳总结出一种分析CO2/CH4在褐煤中吸附机理的方法。本发明通过分析褐煤结构空隙特性、CO2/CH4稳定吸附构型、CO2/CH4与褐煤结构的相互作用能以及相互作用方式、CO2/CH4在褐煤结构中的吸附位,对比分析CO2/CH4单组分和混合竞争吸附的差异,利用Langmuir等温吸附线以及PSD曲线研究了温度、压强、拓扑结构、孔径尺寸、静电贡献、官能团等影响CH4和CO2在褐煤结构中吸附的本质因素,详细的流程图见说明书附图。本发明有助于更加深入地理解CO2、CH4及其混合物在褐煤结构中的竞争吸附作用,对煤层气开发和碳俘获与封存的发展和应用具有一定理论和实践指导意义。
Description
技术领域
本发明涉及煤层气的开发与利用研究领域,具体涉及一种分析CO2/CH4在褐煤中吸附机理的方法,其特征是利用密度泛函理论(DFT)、分子动力学模拟(MD)和巨正则蒙特卡洛模拟(GCMC)研究温度、压强、孔径尺寸、静电贡献、官能团等因素对褐煤中CO2/CH4吸附的影响。
背景技术
由人类活动导致的温室气体浓度增加所引起的全球气候变化问题已经成为影响世界经济秩序、政治格局和国际关系的一个重要因素以及决定世界能源前景的关键。减少温室气体排放除了提高能源利用效率、利用清洁能源和清洁燃烧技术外,另外的重要途径是二氧化碳的捕获和封存(CCS技术)。将CO2注入煤层不仅能提高煤层气的采收率,同时也能缓解温室效应,具有良好的经济效益和环境效益,这一技术涉及到CH4/CO2二元气体在煤层中的吸附过程,故研究CH4/CO2在褐煤结构中的吸附行为,压强、温度、孔径尺寸以及静电作用对吸附的影响对于煤炭开采和瓦斯防治具有重要的指导意义。分子动力学方法(MD)是一种计算机模拟实验方法,是研究凝聚态系统的有力工具。它是对理论计算和实验的有力补充,广泛应用于材料科学、生物物理和药物设计等。借助于蒙特卡洛方法(GCMC)我们可以找出能量最优点,而宏观性质可以通过系综平均得到。蒙特卡洛方法被广泛用于环境和材料科学以及化工吸附分离等行业,是平衡物系研究、物质微观结构、以及新材料制备和改性的最佳工具之一。DFT体系本身数值计算方法取得了巨大发展,使得DFT成为材料领域最先进、最高效的理论研究方法之一。通过DFT方法可以计算材料的物理化学性质,借以比较各种近似理论验证其正确性。本发明中我们结合密度泛函理论(DFT)、分子动力学模拟(MD)和蒙特卡洛模拟(GCMC)三种理论计算方法,在分子模拟水平上研究了CO2/CH4气体在褐煤结构中的吸附行为,考察了温度、压强、静电贡献、孔径尺寸、官能团等因素对煤层中CO2/CH4吸附的影响,通过对比褐煤结构中CO2/CH4的单组分与竞争吸附行为,揭示了CO2/CH4在褐煤结构中吸附差异的本质因素。通过归纳总结,形成一种分析CO2/CH4在褐煤中吸附机理的方法。
发明内容
本发明方法利用密度泛函理论(DFT)、分子动力学模拟(MD)和蒙特卡洛模拟(GCMC)三种理论计算方法,研究了CO2/CH4气体在褐煤结构中的吸附行为,通过计算CO2/CH4绝对吸附量、超额吸附量、选择吸附比,吸附热,拟合并对比分析Langmuir等温吸附线以及PSD曲线,考察温度、压强、静电贡献、拓扑结构、孔径尺寸等因素对煤层中CO2/CH4吸附的影响。结合实验研究文献,建立一种分析CO2/CH4在褐煤中吸附机理的方法。本发明方法所采用的计算机模拟软件为美国Accelrys公司研发的Materials Studio(MS)软件和美国西北大学的多功能计算程序Multipurpose Simulation Code简称MuSiC以及Gaussian 09。MS的Discover模块是分子力学计算引擎,它使用多种分子力学和动力学方法,以仔细推导的力场作为基础,可准确地计算出最低能量构型、分子体系的结构和动力学轨迹等。
本发明的内容可分为褐煤结构模型构建及孔隙特性分析和CO2/CH4在褐煤结构中吸附机理研究两部分,具体内容如下:
(1)褐煤结构模型构建及孔隙特性分析
依据收集褐煤的煤层深度、煤层煤阶、煤层密度、煤岩组分数据,在MS中合理构建褐煤大分子结构模型;
在Gaussian 09中采用DFT方法优化分子结构,保证分子的能量达到平衡状态,采用NBO电荷计算褐煤分子结构的分子电子结构、电子密度分布、势能图进一步考察褐煤分子结构的电荷分布情况,分析褐煤分子结构的稳定性;
在MS中构建无定形态褐煤结构,并采用分子动力学模拟(MD)方法优化无定形态褐煤结构;采用Sarkisov和Duren方法评估初始模型的孔隙结构,计算孔隙度、有效孔隙体积、比表面积、孔隙尺寸分布、孔隙最大直径和孔隙限制直径;
(2)CO2/CH4在褐煤结构中吸附机理研究
在Gaussian 09中采用DFT方法,研究CO2/CH4与褐煤结构的相互作用能以及相互作用方式,分析CO2/CH4在褐煤结构表面的选择吸附位;采用DFT方法,研究CO2/CH4在褐煤结构表面的吸附类型和吸附构型并根据吸附能大小比较CO2/CH4吸附构型的稳定性;根据吸附能考察CO2/CH4在褐煤结构表面的吸附能力及其强弱;
在MuSiC中采用GCMC方法,计算单组分CO2/CH4绝对吸附量和超额吸附量,拟合Langmuir等温吸附线并拟合Langmuir吸附公式,考察温度、压强、静电贡献对单组份CO2/CH4在褐煤结构中吸附的影响;计算单组份CO2/CH4的吸附热,比较CO2/CH4吸附能力的强弱;计算褐煤结构中CO2/CH4混合气体吸附量、静电作用贡献以及选择吸附比,研究竞争吸附中CO2/CH4吸附差异以及CO2/CH4的竞争吸附填充过程,判断CO2/CH4的吸附位,考察温度、压强、静电贡献、官能团对CO2/CH4混合气体竞争吸附的影响,揭示CO2/CH4的竞争吸附机理。
本发明的有益效果在于:提出了一种分析CO2/CH4在褐煤中吸附机理的简便方法。本发明方法计算快速、结果准确,能从分子水平层次认识褐煤结构褐煤结构中CH4和CO2的吸附机理,并确定温度、压强、拓扑结构、孔径尺寸、静电贡献等因素对吸附性能的影响,可以为煤炭开采、瓦斯防治、设计新型碳捕获材料提供理论指导,缩短研发周期,同时也可以节省大量的人力、物力和财力。
附图说明
下面结合附图和实施案例对本发明进一步说明。
图1:具体研究技术路线流程图
具体实施方式
实例基于周期性的无定形褐煤结构研究了CH4/CO2在褐煤中的单组份以及竞争吸附过程,通过密度泛函理论(DFT)对电荷密度分析、分子结构稳定性分析、相互作用能及相互作用方式分析、选择吸附位分析、吸附类型和吸附构型稳定性分析,通过分子动力学模拟(MD)对孔隙特性分析,通过巨正则蒙特卡洛模拟(GCMC)对吸附量分析、静电贡献分析、选择性吸附分析、吸附热曲线分析、官能团对吸附位影响分析、温度对吸附的影响分析、压强对吸附的影响分析,研究了温度、压强、拓扑结构、孔径尺寸、静电贡献、官能团等影响CH4和CO2吸附的本质因素。具体实施方式分为以下4步:
(1)无定形褐煤结构模型构建
利用MS软件构建褐煤分子结构,本实例中选用C39H37NO10S这一种分子,使用COMPASS力场和Smart Minimizer方法优化分子结构;在此基础上采用正交晶系和周期性边界条件构建无定形褐煤结构,晶胞长为狭缝孔的平均体积密度保持在1033kg/m3;CH4和CO2以及体系表面的非成键作用包括范德瓦尔斯作用和静电相互作用,截断半径为采用NPT系综进行MD模拟,设置压力为0.01Gpa,模拟时间步长为1.0fs,动力学时间为1000ps;
(2)第一性原理电荷计算
使用Gaussian09计算原子电荷,稳定构型计算涉及到吸附质(褐煤分子)以及吸附剂(CH4、CO2);采用DFT方法优化分子结构,保证分子的能量达到平衡状态,采用NBO电荷计算褐煤分子结构的分子电子结构、电子密度分布、势能图进一步考察褐煤分子结构的电荷分布情况,分析褐煤分子结构的稳定性;构型优化和能量计算采用DFT,6-31G+(d,p),WB97XD,NBO基组;
(3)巨正则蒙特卡洛(GCMC)模拟
在MuSiC中进行GCMC模拟,采用Peng–Robinson方程式来计算模拟压力,采用298K、313K、373K以及0.05-10MPa来模拟地层条件的温度和压强;模拟的迭代次数为1.0×108;用Lennard-Jones(LJ)势函数和库仑力作用来计算CH4-CH4、CH4-孔隙壁、CH4-CO2、CO2-CO2、以CO2-孔隙壁之间的相互作用。原子的LJ参数取自经典普遍力场(UFF);
(4)性质分析
计算CO2/CH4绝对吸附量,拟合Langmuir绝对等温吸附线并拟合Langmuir吸附公式采用ne=na-Vpρb公式计算CO2/CH4超额吸附量,拟合Langmuir超额等温吸附线,采用公式计算静电作用贡献,计算选择吸附比、吸附热,拟合选择吸附比曲线,考察静电作用贡献对CH4/CO2吸附的影响以及CO2/CH4的竞争吸附填充过程;研究不同温度、压力条件下CO2/CH4在无定形褐煤结构中的吸附行为,考察温度、压力对吸附性能的影响;拟合径向分布函数曲线,作CO2/CH4吸附密度分布图,研究CH4、CO2的吸附位,考察管能团对CO2/CH4吸附位的影响;拟合PSD曲线,考察拓扑结构、孔径尺寸对吸附性能的影响;对比分析CO2/CH4单组分和混合竞争吸附的差异,归纳总结影响CH4、CO2的吸附的本质因素,揭示CO2/CH4的竞争吸附机理。
以上所述,仅为本发明的一则实施例而已。并非对本发明作任何形式上的限制;凡熟悉本专业的普通技术人员均可按说明书附图及以上所述步骤而顺畅地实施本发明;但是,凡熟悉本专业的技术人员在不脱离本发明技术方案范围内,可利用以上所阐述的技术内容而做出的些许更改、修饰与演变的等同变化,均为本发明的等效实施例;同时,凡依据本发明的实施技术对以上实施例所作的任何等同变化的更动、修饰与演变等,均属于本发明的技术方案的保护范围之内。
Claims (1)
1.一种分析CO2/CH4在褐煤中吸附机理的方法,该方法是基于周期性的密度泛函理论(DFT)、分子动力学模拟(MD)和巨正则蒙特卡洛模拟(GCMC)三种理论计算方法,研究了CH4和CO2在褐煤中的单组份以及竞争吸附过程,通过电荷密度分析、分子结构稳定性分析、相互作用能及相互作用方式分析、选择吸附位分析、吸附类型和吸附构型稳定性分析、孔隙特性分析吸附量分析、静电贡献分析、选择性吸附分析、吸附热曲线分析、官能团对吸附位影响分析、温度对吸附的影响分析、压强对吸附的影响分析,研究了温度、压强、拓扑结构、孔径尺寸、静电贡献、官能团影响CH4和CO2吸附的本质因素,进一步结合实验研究文献,具体实施方式分为以下4步:
(1)无定形褐煤结构模型构建
利用Materials Studio软件构建褐煤分子结构,使用COMPASS力场和Smart Minimizer方法优化分子结构,在此基础上采用正交晶系和周期性边界条件构建无定形褐煤结构;CH4和CO2以及体系表面的非成键作用包括范德瓦尔斯作用和静电相互作用,采用NPT方法进行MD模拟;截断半径为
(2)第一性原理电荷计算
使用Gaussian 09中计算原子电荷,稳定构型计算涉及到吸附质褐煤分子以及吸附剂CH4、CO2;采用DFT方法优化分子结构,保证分子的能量达到平衡状态,采用NBO电荷计算褐煤分子结构的分子电子结构、电子密度分布、势能图进一步考察褐煤分子结构的电荷分布情况,分析褐煤分子结构的稳定性;构型优化和能量计算采用DFT,6-31G+(d,p),WB97XD,NBO基组;
(3)巨正则蒙特卡洛(GCMC)模拟
在Multipurpose Simulation Code中进行GCMC模拟,采用Peng–Robinson方程式来计算模拟压力,采用298K、313K、373K以及0.05-10MPa来模拟地层条件的温度和压强;模拟的迭代次数为1.0×108;用Lennard-Jones(LJ)势函数和库仑力作用来计算CH4-CH4、CH4-孔隙壁、CH4-CO2、CO2-CO2、以CO2-孔隙壁之间的相互作用;原子的LJ参数取自经典普遍力场(UFF);
4)性质分析
计算CH4/CO2绝对吸附量,拟合Langmuir绝对等温吸附线并拟合Langmuir吸附公式,计算CH4/CO2超额吸附量,拟合Langmuir超额等温吸附线,计算静电作用贡献,计算选择吸附比、吸附热,拟合选择吸附比曲线,考察静电作用贡献对CH4/CO2吸附的影响以及CH4/CO2的吸附填充过程;研究不同温度、压力条件下CH4/CO2在无定形褐煤结构中的吸附行为,考察温度、压力对吸附性能的影响;拟合径向分布函数曲线,作CH4/CO2吸附密度分布图,研究CH4、CO2的吸附位,考察管能团对CH4/CO2吸附位的影响;拟合PSD曲线,考察拓扑结构、孔径尺寸对吸附性能的影响;对比分析单组分和混合竞争吸附的差异,归纳总结影响CH4、CO2的吸附的本质因素,揭示CO2/CH4的竞争吸附机理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510202355.3A CN104778329B (zh) | 2015-04-24 | 2015-04-24 | 一种分析co2/ch4在褐煤中吸附机理的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510202355.3A CN104778329B (zh) | 2015-04-24 | 2015-04-24 | 一种分析co2/ch4在褐煤中吸附机理的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104778329A CN104778329A (zh) | 2015-07-15 |
CN104778329B true CN104778329B (zh) | 2017-11-28 |
Family
ID=53619790
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510202355.3A Expired - Fee Related CN104778329B (zh) | 2015-04-24 | 2015-04-24 | 一种分析co2/ch4在褐煤中吸附机理的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104778329B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507136A (zh) * | 2018-10-25 | 2019-03-22 | 南京邮电大学 | 一种二硫化钼与亚甲基蓝之间吸附机理的研究方法 |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105138809A (zh) * | 2015-10-21 | 2015-12-09 | 中国石油大学(华东) | 一种无定形态烟煤结构孔隙特性的分析方法 |
CN105834447B (zh) * | 2016-04-27 | 2017-09-22 | 吉林大学 | 基于密度泛函理论揭示纳米双金属CoNi吸附机理的方法 |
CN106951713B (zh) * | 2017-03-24 | 2020-05-19 | 中国石油大学(北京) | 一种模拟化学助剂作用下甲烷水合物分解过程的方法 |
CN107729650B (zh) * | 2017-10-13 | 2020-12-25 | 谭沛 | 一种涂层、防污材料及其分子设计方法 |
CN107918720B (zh) * | 2017-11-16 | 2021-07-23 | 中国石油大学(华东) | 定量分析分子筛中正戊烷/异戊烷分离的力场参数的方法 |
EP3614389B1 (en) * | 2018-08-23 | 2023-10-11 | Tata Consultancy Services Limited | Systems and methods for predicting structure and properties of atomic elements and alloy materials thereof |
CN109473144B (zh) * | 2018-10-30 | 2021-07-09 | 中国石油大学(华东) | 一种构建高分子与超临界二氧化碳混合体系全原子模型的方法 |
CN109920487B (zh) * | 2019-01-17 | 2021-04-02 | 中国科学院合肥物质科学研究院 | 一种辐照环境下钨中氢滞留与脱附的多尺度模拟方法 |
CN110762950B (zh) * | 2019-09-18 | 2020-08-14 | 宁波大学 | 高电压对多孔石墨烯进行除湿的方法、模拟方法和装置 |
CN112113881B (zh) * | 2020-09-01 | 2021-05-07 | 中国矿业大学(北京) | 煤中吸附性甲烷的密度测量方法 |
CN112132337A (zh) * | 2020-09-22 | 2020-12-25 | 中国矿业大学(北京) | 一种预测煤层瓦斯参数的方法 |
CN112229754B (zh) * | 2020-09-29 | 2022-01-28 | 西南石油大学 | 一种迭代平均法计算甲烷在页岩中吸附量的方法 |
CN112850843B (zh) * | 2021-01-07 | 2022-01-18 | 华南理工大学 | 一种全面揭示废水不同有机物光催化降解难易程度的方法 |
CN113591408B (zh) * | 2021-08-02 | 2022-08-30 | 西南石油大学 | 液态硫在碳酸盐岩纳米孔隙中微观吸附特征参数的确定方法 |
CN113779830B (zh) * | 2021-08-31 | 2022-12-09 | 西安交通大学 | 一种纳米材料微孔结构与气体吸附扩散性质间构效关系的获取方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4072192A (en) * | 1976-02-04 | 1978-02-07 | Texaco Inc. | Sulfomethylated lignite salt as a sacrifical agent in oil recovery processes |
CN102831256A (zh) * | 2011-06-16 | 2012-12-19 | 中国石油化工股份有限公司 | 一种利用计算机模拟计算化学物质溶解度参数的方法 |
CN103801272A (zh) * | 2014-01-27 | 2014-05-21 | 西安科技大学 | 一种煤/聚乙烯亚胺交联复合螯合吸附剂的制备方法 |
-
2015
- 2015-04-24 CN CN201510202355.3A patent/CN104778329B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4072192A (en) * | 1976-02-04 | 1978-02-07 | Texaco Inc. | Sulfomethylated lignite salt as a sacrifical agent in oil recovery processes |
CN102831256A (zh) * | 2011-06-16 | 2012-12-19 | 中国石油化工股份有限公司 | 一种利用计算机模拟计算化学物质溶解度参数的方法 |
CN103801272A (zh) * | 2014-01-27 | 2014-05-21 | 西安科技大学 | 一种煤/聚乙烯亚胺交联复合螯合吸附剂的制备方法 |
Non-Patent Citations (1)
Title |
---|
The adsorption behavior of CH4 on microporous carbon:effects of surface heterogeneity;Jin D 等;《Physical Chemistry Chemical Physics》;20141231;第16卷(第22期);第11037-11046页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507136A (zh) * | 2018-10-25 | 2019-03-22 | 南京邮电大学 | 一种二硫化钼与亚甲基蓝之间吸附机理的研究方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104778329A (zh) | 2015-07-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104778329B (zh) | 一种分析co2/ch4在褐煤中吸附机理的方法 | |
Collell et al. | Molecular simulation of bulk organic matter in type II shales in the middle of the oil formation window | |
Wu et al. | Molecular dynamics simulations about adsorption and displacement of methane in carbon nanochannels | |
Babatunde et al. | Molecular simulation study of CO2/CH4 adsorption on realistic heterogeneous shale surfaces | |
Zhang et al. | Adsorption and selectivity of CH4/CO2 in functional group rich organic shales | |
Wu et al. | Molecular simulation of flue gas and CH4 competitive adsorption in dry and wet coal | |
Niu et al. | Influencing factors and selection of CH4 and CO2 adsorption on Silurian shale in Yibin, Sichuan province of China | |
Fan et al. | Effects of N and S functionalities on binary gas co-adsorption on a coal macromolecule | |
Hao et al. | Thermodynamic analysis of CH4/CO2/N2 adsorption on anthracite coal: investigated by molecular simulation | |
Qin et al. | Time-and pressure-independent gas transport behavior in a coal matrix: model development and improvement | |
Zhou et al. | Edge-functionalized nanoporous carbons for high adsorption capacity and selectivity of CO2 over N2 | |
Zhu et al. | Molecular model construction of Danhou lignite and study on adsorption of CH 4 by oxygen functional groups | |
Lee et al. | Optimal separation of CO2/CH4/Brine with amorphous Kerogen: A thermodynamics and kinetics study | |
Peng et al. | Control mechanism of small organic molecules on methane adsorption capacity of coal | |
Long et al. | Nanomechanical properties of CH4-containing coal during CO2 storage under different injection pressures based on molecule dynamics | |
Li et al. | How gas recovery and carbon storage capacity response to dynamic deformation of kerogen upon CO2/CH4 competitive adsorption for CCUS? Evidence from molecular dynamics | |
Li et al. | High-throughput computational screening of hypothetical metal–organic frameworks with open copper sites for CO 2/H 2 separation | |
Zhang et al. | Molecular simulation of enhanced CH4 recovery and CO2 storage by CO2–N2 mixture injection in deformable organic micropores | |
Huang et al. | Molecular simulation of CO2 sequestration and enhanced gas recovery in gas rich shale: An insight based on realistic kerogen model | |
Tao et al. | Displacement of shale gas confined in illite shale by flue gas: A molecular simulation study | |
Meng et al. | CO2 utilization and sequestration in organic-rich shale from the nanoscale perspective | |
Zhang et al. | Molecular insights into recovery of shale gas by injecting CO2-rich industrial waste gas in kerogen slit nanopores | |
Wang et al. | Feasibility study on using CO2-rich industrial waste gas replacement for shale gas exploration based on adsorption characteristics | |
Lin et al. | The influencing factors of gas adsorption behaviors in shale gas reservoirs | |
Zhang et al. | Insights into influence of CH4 on CO2/N2 injection into bituminous coal matrix with biaxial compression strain in molecular terms |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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: 20171128 Termination date: 20200424 |
|
CF01 | Termination of patent right due to non-payment of annual fee |