[发明专利]一种基于雷诺数的环空压力计算方法有效
申请号: | 201910493699.2 | 申请日: | 2019-06-06 |
公开(公告)号: | CN110197040B | 公开(公告)日: | 2023-04-07 |
发明(设计)人: | 冯福平;丛子渊;严茂森;胡小康;刘圣源;张祯 | 申请(专利权)人: | 东北石油大学 |
主分类号: | G06F30/20 | 分类号: | G06F30/20 |
代理公司: | 哈尔滨东方专利事务所 23118 | 代理人: | 曹爱华 |
地址: | 163319 黑龙江*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: |
本发明涉及的是一种基于雷诺数的环空压力计算方法,它包括:一、通过现场测试或实验收集注采井的基本参数;二、计算注采井环空保护液下部泄露流量;三、计算窜流气体脱离水泥环表面时的初生气泡半径;四、判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;五、计算环空保护液中气体上升速度;步骤六:计算环空带压值,环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力后,若 |
||
搜索关键词: | 一种 基于 雷诺数 压力 计算方法 | ||
【主权项】:
1.一种基于雷诺数的环空压力计算方法,其特征在于:步骤一:通过现场测试或实验收集注采井的基本参数,包括:地层的温度压力、水泥环孔隙度、水泥环渗透率、井口和井筒温度、环空保护液的相关物性参数和窜流气体的相关物性参数;步骤二:计算注采井环空保护液下部泄露流量qg:
式中:K为水泥环渗透率,m2;A为环空截面积,m2;Tsc为标况温度,K;μg为气相粘度,Pa·s;Tf为地层温度,K;Lf为水泥环长度,m;Psc为标况压力,Pa;Pf为地层压力,Pa;Pc为环空保护液下部面压力,Pa;步骤三:计算窜流气体脱离水泥环表面时的初生气泡半径rb:
式中:ρl为环空保护液密度,kg/m3;ρg为窜流气体密度,kg/m3;g为重力加速度,9.8m/s2;μl为环空保护液粘度,Pa·s;σ为环空保护液与气体的表面张力系数,N/m;φ为水泥环孔隙度,无量纲;
为水泥颗粒直径,m;步骤四:判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;气泡雷诺数公式为:
①当气泡雷诺数Re<1时,选用小雷诺数气泡上升速度计算公式:
式中:CD为阻力系数,其值为24/Re。②当气泡雷诺数Re>1时,选用大雷诺数气泡上升速度计算公式:
步骤五:计算环空保护液中气体上升速度vt,w:
其中:![]()
![]()
a3=‑(a1+a2+a4)
式中:α为含气率,无量纲;n为流性指数,无量纲;步骤六:计算环空带压值Pt:
式中:Vm为环空保护液体积,m3;Vt为气顶体积,m3;Cm为环空保护液压缩系数,Pa‑1;Twh为井口温度,K;A为环空截面积,m2;Twb为井筒温度,K;Lf为环空保护液高度,m;环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力Pt1后,若Pf=Pc,则计算结束;若Pf≠Pc,则i=i+1,返回步骤二重新计算,直至地层压力Pf与环空保护液下部压力Pc相等。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于东北石油大学,未经东北石油大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201910493699.2/,转载请声明来源钻瓜专利网。