《采气工程作业样本.doc》由会员分享，可在线阅读，更多相关《采气工程作业样本.doc（11页珍藏版）》请在金锄头文库上搜索。
1、资料内容仅供您学习参考，如有不当或者侵权，请联系改正或者删除。第三章作业3-1解: ( 1) 当气井流动符合达西流动条件时, 由已知得, 计算绝对无阻流量AOF将代入上式, 得当时( 2) 当速度系数和表皮系数分别为时, 由公式, 得 ( 其中) 紊流系数 计算绝对无阻流量AOF当时, 则此时当时则此时3-2解: 由题中常规回压试井数据得到下表测试点井底流压( MPa) 产气量Qsc( 万方/天) (Qsc)lg( Qsc) 关井2.814512.779412.135147.2582251.22.716626.22687.48841.32.609844.0121937.0561441.42.5
2、00157.43294.761.总和139.7676066.5627695.P=Pr-Pwf( Pr-Pwf) /Qsclg( P) ( lg( P) ) lg(Qsc)lg(P)0.196345890.-0.0.-0.0.541494690.-0.0.-0.1.110354210.0.0.0.1.670910240.0.0.0.3.519105030.-0.0.-0.由上表可得到以下数据: ( 1) 指数式产能方程则气井的流入动态曲线为: ( 2) 二项式产能方程则第四章作业4-2某气井气层中部井深3000m, 井口温度为25, 气层温度为110, 气体相对密度0.65。液, 气、 水、 油
3、产量分别为10104m3/d、 10m3/d、 3m3/d, 井口油压2MPa, 油的相对密度0.85, 水的相对密度1.05, 油管内径62mm。试用Shiu & Beggs方法计算该井的流温梯度曲线。解: ( 1) 计算 ( 2) 计算松弛距离A由, 得 ( 3) 计算井筒静温梯度由温度计算公式可知, 代入上式, 得综上所述, 该井的流温梯度曲线为: 4-3已知某气井的参数: 井中部深度H3500m, 油管尺寸为2( 76.2mm) , 井筒平均温度 342K, 天然气相对密度为0.6, 地层压力为28MPa, 井口压力为6.5MPa, 气井产能方程为 。试利用节点分析方法, 计算气井产能
4、。解: 以井底为求解点, 求上述条件下的井底流压( 1) 经过气井的产能方程, 绘出该气井的流入动态曲线, 由已知得, 地层压力为28MPa, 则计算绝对无阻流量给定任意产量, 求该条件下的井底流压: 产气量Qsc( 万方/天) 井底流压Pwf( MPa) 028126.78058999225.31402773323.55419283421.42895238518.81488772615.46609194710.714476197.860绘出流入动态曲线为: ( 2) 根据已知条件, 计算流出动态曲线给定某一产量Qsc, 计算该条件下的井底流压: （一） 计算初值（二） 计算平均参数（三） 计
5、算由公式可计算出（四） 计算（五） 计算无量纲量S（六） 计算井底流压（七） 比较若, 则停止计算, 令否则, 进行第二次计算直到为止。计算程序如下: Private Sub Command1_Click()qsc = Text5pwh = Text1H = Text2D = Text3 * 10 -3Pc = 4.626Tc = 205.6f = 0.015pwf0 = pwh * (1 + 0.00008 * H)DoPm = (pwf0 + pwh) / 2Tm = Text4计算压缩因子Zz = Z1(Pm, Tm)计算粘度miugMg = 28.97 * 0.6roug = Pm *
6、 Mg / (z * 0.008314 * Tm * 1000)K = (9.4 + 0.02 * Mg) * (1.8 * Tm) 1.5) / (209 + 19 * Mg + 1.8 * Tm)X = 3.5 + 986 / (1.8 * Tm) + 0.01 * MgY = 2.4 - 0.2 * Xmiug = 10 -4 * K * Exp(X * roug Y)计算ReBg = 3.447 * 10 -4 * z * Tm / Pmv = qsc * 10 4 * Bg / (86400 * 3.14 * D 2 * 0.25)rou = 3484.4 * 0.6 * Pm /
7、(z * Tm)Re = 1000 * rou * v * D / miug计算Ss = 0.03415 * 0.6 * H / (z * Tm)计算井底压力pwf1 = (pwh 2 * Exp(2 * s) + 1.324 * 10 -18 * f * (qsc * 10000 * z * Tm) 2 * (Exp(2 * s) - 1) / D 5) 0.5com = Abs(pwf1 - pwf0)pwf0 = pwf1Loop Until com 0.001Text6 = pwf0End Sub计算压缩因子Z的函数Z1Public Function Z1!(ByVal Pm#, By
8、Val Tm#)Dim pr#, tr#Pc = 4.626Tc = 205.6tr = Tm / Tcpr = Pm / Pc z0 = 1 Doden1 = 0.27 * pr / (z0 * tr)Z1 = 1 + (0.3051 - 1.0467 / tr - 0.5783 / tr 3) * den1 + (0.5353 - 0.6123 / tr + 0.6816 / tr 3) * den1 2com = Abs(Z1 - z0)z0 = Z1Loop Until com 0.001Z1 = z0End Function经过以上程序, 计算可得到以下数据: 产气量Qsc( 万方/天) 井底流压Pwf( MPa) 08.2124336418.2389972528.2389972538.2720511148.3180857858.3768679368.4481216778.5314845988.6265795198.73298701108.85026067绘出气井的流入、 流出动态曲线为: 由此可确定该井的产能: ( 此时的井底流压为9.351MPa) 第五章作业5-1.解: ( 1) 应用Dupuit公式计算临界产量由, 得其中则( 2) 应用Schols公式计算临界产量由, 得( 3) 计算水锥突破时间则
