• / 6
  • 下载费用:1 下载币  

盐穴地下储气库溶腔形态变化数值模拟

关 键 词:
地下 储气库溶腔 形态 变化 数值 模拟
资源描述:
· 88· 天然气工业 2014年7月盐穴地下储气库溶腔形态变化数值模拟梁光川1 王梦秋11.西南石油大学石油与天然气工程学院彭星煜1 张硕22.中国石油西气东输管道公司苏北管理处梁光川等.盐穴地下储气库溶腔形态变化数值模拟.天然气工业,2014,34(7):88—92.摘 要 盐岩蠕变对盐穴地下储气库的溶腔形态有着较大的影响,易造成储气库失效、泄漏、破裂等安全问题。为此,利用美国速拉格朗日差分软件,建立了盐穴储气库模型,对储气库溶腔蠕变规律进行了数值模拟,分析了在不同蠕变时间内、不同内压下储气库溶腔在不同直径处的变形量,得到在盐岩蠕变下,盐穴储气库溶腔形态的变化规律,并给出了相应操作建议:①盐穴储气库溶腔不同直径处产生的变形量是不同的,溶腔最大直径处变形量最大,随着直径的增大变形量增大,且呈非均匀性增大,因此在建设盐穴储气库时,应对储气库溶腔的高径比进行严格地计算;②溶腔内压对储气库溶腔变形量有直接作用,较高压力下的盐岩蠕变对储气库溶腔形态变化的影响比低压时小,因此建议在盐穴储气库实际运行过程中,尽量避免低压运行,以保证储气库溶腔形态的完整性。该成果为盐穴储气库的安全运行与维护提供了可行的依据。关键词 盐穴地下储气库 盐岩蠕变 溶腔形态 变化规律 数值模拟 溶腔内压 完整性 操作压力0.3787/j.000—0976.2014.07.015N of of .10500,.25009,4,,8—92,7/25/2014.(000—0976;of to of in n of to of in of is n of in of he a on of of is of of to of be in of of be in of 光川,1972年生,教授,博士;长期从事油气地面集输系统节能改造、天然气地下储气库优化、油气管道输送理论及技术研究工作。地址:(610500)四话:13981827236。E—63.输与加工盐穴地下储气库(以下简称盐穴储气库)通过常规钻井方法钻穿盐岩层,注入淡水进行冲蚀,抽出卤水后利用形成的溶腔进行气体储存。由于盐岩具有孑透性小、塑性大等优点,运用盐穴储气库储存气体安全可靠、不易渗漏[1]。盐岩具有较大塑性,造成了储气库溶腔注采气时的可收缩性。同时,盐岩具有蠕变特性,盐岩蠕变对盐穴地下储气库的溶腔形态有较大影响,造成储气库失效、泄漏、破裂等安全问题心]。鉴于此,利用美国陕速拉格朗日差分软件,建立盐穴储气库模型,对储气库溶腔蠕变规律进行数值模拟,分析在不同蠕变时间内、不同内压下储气库溶腔在不同直径处的变形量,得到在盐岩蠕变下,盐穴储气库溶腔形态的变化规律,为盐穴储气库的安全运行与维护提供可行的依据。1 建模原理利用数值模拟方法对盐穴储气库溶腔蠕变规律及溶腔变形量开展研究,基本原理是采用拉格朗日法研究每个流体质点随时间变化的情况,着眼于某个流体质点,研究任一段时间内质点运动的轨迹以及质点的速度和压力。求解物体离散为一系列如图1所示的4四面体单元,并采用插值函数d。,一>:占:。N”进行求H=1解,函数中“,、73,分别代表四面体中节点的位移、速度。C¨’数值模拟基本原理是[3]:①通过对三维介质的离散,使所有外力与内力集中于三维网络节点上,进而将连续介质运动定律转化为离散节点上的牛顿定律;②时间与空间的导数采用沿有限空间与时间间隔线性变化的有限差分来近似地表示;③将静力问题当做动力问题来求解,运动方程中惯性项用来作为达到所求静力平衡的一种手段。节t’三维离散四面体图2溶腔模型的建立我国盐岩矿藏以多夹层盐岩为主,由于主要分析盐穴储气库溶腔在盐岩蠕变下的变形规律,因此在建立储气库溶腔模型时将复杂的盐岩夹层简化为简单的单层盐岩层。借鉴江苏金坛盐穴储气库的建库参数¨·,建立了如图2所示的腔体几何模型。建立溶腔及围岩地层模型后,将地层岩石力学参数、盐层的蠕变关系、溶腔周围的地应力参数及溶腔运行参数带人模型中进行模拟计算。图2腔体几何模型图2.1 网格划分及参数设置利用。,并接人到3分后的网格如图3所示。图3 网格划分图图3中,l/4圆柱体模拟储气库溶腔周围的地质环境,其直径是储气库溶腔模型最大直径的20倍,在网格中间位置可以看到储气库溶腔轮廓的1/4。在网万方数据天 然 气 T 川,格划分时将储气库溶腔网格一起划分,有助于储气库溶腔的形成。在此网格中,定义了2个组,分别是者对储气库溶腔进行网格划分,后者对储气库溶腔周围地质环境进行网格划分。选择采用各向同性弹性模型,其材料特性常数见表1。表1 材料特性常数参数表参数 参数值体积模量/kg· 1095 1082 5002.2模型求解及验证数值分析采用迭代方法计算,在迭代过程中监控一些变量或参数的变化,以判断分析结果是否正确,模型是否与实际相符¨。7]。监控模型最大不平衡力,所谓体系最大不平衡力,是指每一个计算循环中,外力通过网格节点传递分配到体系各节点时,所有节点的外力与内力之差中的最大值。如图4~6监控了最大不平衡力,并记录了位移等值线图和应力等值线图。4 54 0互3 5兰3.0滚2 57 2 0二<1 5O 5O l O 1.5 2 01图4最大不平衡力图由图4可以得出,该模型中最大不平衡力逐渐收敛至零,表示所建储气库溶腔模型具有收敛性。由图4、5可以看出,在建立的土体模型中,由于重力场的作用,最大位移和最大应力分别出现在模型的最底部。2.3腔体的形成及导入蠕变模型利用空模型形成腔体模型,通过设置挖后的模型如图7、8所示。形成腔体模型后,需要对溶腔施加内压,通过施加法向应力将储气库内压分别设置为10 5 7 0 3 5 择经典黏弹性蠕变模型[8。91计算在不同内压条件下溶腔不图5位移等值线图图6应力等值线图移/0~00 10~0O 40~00 60~00 80~11 00~0、5.:一一r—r~一一一78505050544332,一一一一一一一一一■—典黏弹性蠕变模型材料参数如表2所示。表2经典黏弹性模型的材料参数表参数 参数值弹性体积模量/Pa·s)1 1095 10811 010” 2 4 6 8 10 1 5 20蠕变叫M/压为10 变计算及结果分析 图l,。时储气库溶腔不同直径处的变形量图3.1蠕变计算根据一些学者的研究[1“12j,进行蠕变计算时,分别计算了储气库溶腔内压为10 5 7 0 3 5 0 变时间分别为2a、4 a、6 a、8 a、10 a、15 a、20 气库溶腔直径为20 m、30 m、40 m、50 .2结果分析根据模拟计算所得结果,可以得到不同压力下储气库溶腔不同直径处变形量图(图9~11),经总结分析,仅列出具有代表性的分析图。图9~11分别表示溶腔内压为10 0 0 同蠕变时间盐穴储气库溶腔不同直径处的变形量。对比分析图9~11,可以得出:①溶腔不同直径处变形量随时间变化规律比较一致,溶腔直径为50 蠕变时间为6 同压力下溶腔不同直径处变形量曲线均出现一个明显的转折点。o 2 4 6 8 10 15 20蠕变压为30 腔变形量曲线斜率较大,表明在此时间内溶腔变形速率较大,变形比较显著;蠕变时∞"如笛加博j。斟蜜567890●,一34_)599999~~一一一一~一7056789O●2345;999999一,,,,■■—32092 · 天然气工业 2014年7月间大于6 腔变形量曲线较平缓,溶腔变形量随时间变化不大;③内压为10 变20 腔直径为50 压为20 变20 腔直径为50 1压为30 变20 腔直径为50据表明,溶腔内压小于20腔变形量随压力增大显著增加;溶腔内压大于20 腔变形量变化不大。4 结论1)建立了盐穴储气库溶腔模型,计算分析了在盐岩蠕变下盐穴储气库溶腔形态的变形量,得到了盐穴储气库在不同压力下溶腔不同直径处的变化规律。由于受盐岩蠕变的影响,在盐穴储气库溶腔不同直径上产生的变形量是不同的,溶腔最大直径处变形量最大,随着直径的增大变形量增大,且呈非均匀增大,因此在建设盐穴储气库时,应对储气库溶腔的高径比进行严格的计算。2)盐穴储气库在不同的运行压力下变形量是不同的。当压力小于20 着内压增加,溶腔整体变形量增加;当内压大于20 着内压增加,溶腔整体变形量基本保持不变。说明在较高压力下,盐岩蠕变对溶腔形态变化的影响比低压时小。因此,建议在盐穴储气库实际运行过程中,尽量避免低压运行,以保证储气库溶腔形态完整性。参 考 文 献国地下储气库现状与发展展望[然气工业,2006,26(6):111].006,26(6):111—113.[2]王同涛.多夹层盐岩体中储气库围岩变形规律及安全性研究[D].北京:中国石油大学,2011.on of in ].01 1.[31丁国生.盐穴储气库溶腔收缩规律分析[J].天然气工业,2007,27(11):94—96.he of ].007,27(11):94—96.[41尹雪英,杨春和,陈剑文.金坛盐矿老腔储气库长期稳定性分析数值模拟[土力学,2006,27(6):869—874.on of ].006,27(6):8695]谭羽非.天然气地下储气库技术及数值模拟[M].北京:石油工业出版社,2007.nu—].007.[6]王新胜.盐岩储气库运营期稳定性研究[D].重庆:重庆大学,2009.on of of ].009.[7]李后荣.盐岩地下储气库稳定性分析[J].贵州大学学报,2011,28(1):122—125.he of ].011,28(1):122 125.[8]梁卫国,徐素国.金坛储气库盐岩蠕变特性及其实用本构研究[J].辽宁工程技术大学学报,2007,26(3):354 356.U on of ].007,26(3):354—356.[9]杨春和,陈锋,曾义金.盐岩蠕变损伤关系研究[J].岩石力学与工程学报,2002,14(11):1602—1604.of ].002,14(11):1602 1604.[101 1.008.[11] N,.of in ].013(33):171—185.[12]姜德义.岩盐溶腔稳定性及失稳控制研究[D].重庆:重庆大学,2003.on of ].003.(修改回稿日期 2014—05—13编辑何 明)万方数据盐穴地下储气库溶腔形态变化数值模拟作者: 梁光川, 王梦秋, 彭星煜, 张硕, 梁光川,王梦秋,彭星煜,南石油大学石油与天然气工程学院), 张硕,国石油西气东输管道公司苏北管理处)刊名: 天然气工业英文刊名: (期): 2014,34(7)引用本文格式:穴地下储气库溶腔形态变化数值模拟[期刊论文]014(7)
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:盐穴地下储气库溶腔形态变化数值模拟
链接地址:http://www.oilwenku.com/p-50456.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开