基于数值模拟方法计算天然气地下储气库的渗漏量.docVIP

参考文献 1　陈元千,杨通佑等.石油及天然气储量计算方法.北京:石 油工业出版社,1990 2　黄炳光,刘属知.实用油藏工程与动态分析方法.北京:石 油工业出版社,1998 3　马世煜.凝析气藏开采技术.北京:石油工业出版社,19964　陈元千.油气藏工程计算方法.北京:石油工业出版社, 1990 5　何更生.油层物理.北京:石油工业出版社,1994 6　余元洲,杨广荣等.凝析气藏物质平衡方程的改进与应用.油气地质与采收率,2002;9(4:66~68 (收稿日期　2002-11-05　编辑　居维清 　　3本文为黑龙江省博士后基金资助项目(L RB -KY 01026成果。 　　作者简介:谭羽非,女,副教授,哈尔滨工业大学博士后;主要研究方向为地下储气数值模拟及优化分析,城市煤气负荷预测、控制等。地址:(150090黑龙江省哈尔滨市南岗区海河路204号哈工大二校区2622信箱。电话:(04512303019。E -mail :tanyufei @263.net 基于数值模拟方法计算天然气地下储气库的渗漏量3 谭羽非 (哈尔滨工业大学 　　谭羽非.基于数值模拟方法计算天然气地下储气库的渗漏量.天然气工业,2003;23(2:99~102 摘　要　对于建造在枯竭油气田、含水层的地下储气库,在满足城市调峰的注采动态运行过程中,确保天然气回采率是保证储气库经济合理运行的基础。采用数值模拟与参数自动拟合相结合的算法,在动态确定储气库压力随时间的变化后,通过调整实际的注采量,使计算压力与实测压力吻合,所得到的计算注采量与实际注采量的差值,即为储气库的渗漏量。模拟实例的计算结果表明:所提出的方法能为实际储气库在动态运行时估算渗漏量、确定渗漏路线以及通过调整储运方式来减少渗漏量等提供理论依据。 主题词　地下储气库　渗透　漏失量　数值模拟　计算 　　利用城市附近的枯竭油气田、含水层建造天然气地下储气库,在实际运行过程中常发生渗漏问题。其主要原因:一是由于储库内气体在采出时,压力下降,造成边底水的侵入,这就导致在下一个周期注气时,天然气要沿阻力最小的路径前进,使部分气体进入含水层成为束缚饱和气体,造成内部渗漏;二是在向储库回注气体时,必须以更高的压力将气体从井口注入到储库内,这就极可能使天然气穿过覆盖着 气藏的岩层逃逸,造成天然气的外部渗漏。这两种 原因导致了气库中的实际存气量往往小于“帐面上” 的储气量。如美国怀俄明州的Leroy 地下储气库, 在1975~1982年运行期间,由于储存的天然气大量 流失,不得不于1982年停用〔1〕。 目前国外进行气藏的储量校核主要有以下两种 方法。 (1定体积法 采用估算或计算出的膨胀系数、储库的几何形状、孔隙率及地质资料,用压力对由天然气所占据的孔隙体积进行积分,计算得出天然气量并与“帐面上”的天然气存量进行对比,由差额得出渗漏量。(2压力—储量数据图形分析以虚拟压力p/Z 对储量连续作图,如果储库的体积恒定不变,所画出的是一条直线,其斜率与截距可说明气藏的体积与储量;如果是水流驱动的储气 库,所画出的是一条曲线。 以上方法用于对储气库实际储量的校核,误差 较大。究其原因是由于储库注采气周期交替频繁, 压差变化较大,气体流速快,短时间内注入、采出气 体量大等造成的。 笔者采用数值模拟与参数自动拟合相结合的算 法,在动态确定储库压力随时间的变化后,通过调整? 99?第23卷第2期　　　　　　　　　　　　　　天　然　气　工　业　　　　　　　　　　　　　　工程建设 实际的注采量,使计算压力与实测压力吻合,所得到的计算注采量与实际注采量的差值,即为储库的渗漏量。模拟实例的计算结果表明:所提出的方法能为实际储库在动态运行时估算渗漏量、确定渗漏路线以及通过调整储运方式来减少渗漏量等提供理论依据。 数学模型的建立及求解 　　(1依据连续性方程建立流动方程,整理后的结果如下(详见文献[2]: 　　 (a l Φl +δa HQ Vl /V = 　　H φB l 5S l 5t +S l (c φ+c l 5Φl 5t (1式中:下角标l =g ,气体 w ,水,φ表示孔隙度; Φl = p l -ρl g Z ,称Φl 为势函数;a l =Hkk rl (S μl (p B l (p ,其中:k 表示渗透率;k rl 表示相对渗透率;μl ,B l 分别表示流体黏度、流体体积系数。Q Vl =Q l /ρsc 表示标准状态流体体积流量;c φ、c l 分别表示岩石、流体的压缩系数。在井点处δa =1,在非井点处δa =0。 建立如下辅助方程。 气、水两相饱和度平衡方程:S g +S w =1 (2毛细管力平衡方程:p c (S =Φg -Φw -ΔrZ (3 (2边值条件1初始条件 由于流动过程近似等温,初始条件即为储