自喷井井底压力分布

1、重庆科技学院课程设计报告院（系）：石油与天然气工程学院_ 专业班级：石油普 09-4_09-4_学生姓名：_张鹏_学号:_20094400462009440046_设计地点（单位）重庆科技学院_设计题目: :_ _自喷井井底压力分布完成日期：20122012 年 7 7 月 5 5 日指导教师评语：_成绩（五级记分制）：_指导教师（签字）:_ _一2自喷井井底压力分布张鹏（重庆科技学院，石油工程09 级，学号 2009440046）刖言：2012 年 6 月 25 日至 2012 年 7 月 5 日在徐春碧老师的指导下，我们进行了为期两周的采气工程课程设计， 主要设计内容为气水同产井井底压力分

2、布情况，在老师的指导及整个团队的合作下，我们完成了最初的寻找资料、寻求思路到编程、写报告等一系列的设计。本报告利用 Beggs-Brill 方法和 Duns-Ros 方法求解井底压力。主要目的是利用采油、采气工 程和多项管流的知识来解决现实生产中某一油、气井站具体问题，提高学生解决实际问题的 的能力。设计题目3已知水产量 为 38m /d,GLR=200,油压为 1.8Mpa（表压），油管尺寸为 62mm（21/ 2），1、计算方法1）、持液率H|（）的相关规律2）、混合速度Vm及n的相关计算油层中部深度为 1800m,w1-05,0.65,e=0.016.0.8mPa?s，天然气视临界温度为

3、 225k，视临界压力为4.7Mpa,井口温度为 20C（ 293k），井筒内温度变化为3.5 C/100m，表面张力为0.06N/M,求井底压力分布。（用两种方法做）二、设计方法Beggs-Brill法基本压力公式：P106mg sinfmVm2n/(2D)mVmVsg/ P 106|H|(0)g(1 H)H|(?H|(0)1 C sin 1.81 .sin31.8C (1|)|n(de fN|vNFR9)VmVslVsgQL/A Qg/A3fnexp(s)In y0.0523 3.18ln y 0.8725(ln y)20.01853(1 n y)4gP1/ZRT12、根据各参数的取值范围

4、，可以确定管内液体的流态，具体的标准如表1)、确定管内液体的流态分层L40.01 andNFR=0.01 andNFR=0.01 andL2NFRL1间隔0.01=i0.4 andL2NFR=0.4 andL3NFR=L4分散l=L|l=0.4 andNFRL42)、判断流形的基本公式表1d、e、f、g的确定流型上/下坡defg分离流上坡0.011-3.768：3.539-1.614 I间歇流上坡2.960.305：-0.44730.0978分散流上坡0000各种流型下坡4.7-0.36920.12440.5053、算法框图已知参数P猜测参数F2(0)64、设计程序源代码（见附页程序一）三、设

5、计方法二Duns-Ros法1、计算方法1）、压缩因子的计算Z = 1 +(0.3051 -1.0467 / Tr -0.5783 / Tr3)r+(0.5353-0.6123 / Tr + 0.6816 / Tr3)Tr = T/TcPrP/Pcr0.27 pr/ ZTr2)、气体密度的计算g3484.4gP(i)/(T z)Vsg Qsc (3.447 104Z T/P(i)/SVsl Qw / S3) 、液体粘度计算mexp1.0q0.01479.8( Z i 0.035T T)/2 32 1.9821051.8( z i 0.035T T)/21034) 、平均温度的计算T二(z i 0

6、.035 + 2T) / 2 + 273.154)、四个基本准数的计算NVIvsl4zgNvg = Vsg4|/ gND= Vsg4I/gNI=i4l/ g2、流动形态的判断由课本气液两项管流P36页流动形态分布图可知：I 区：液相是连续的，包括泡状流、弹状流和部分的沫状流；II 区：液相和气相是连续的，包括段塞流和沫状流的其余部分;III 区：气相是连续的，指雾状流。其中71)第 I 区发生在NvgL1LQNVI的范围内。此时Ns2F1F2NVIF3答1 Nvi其中FI,F2, F3和 F4作为 N 的函数在课本 P37图 3-5 图给出。L!和L2为若斯因数作为ND的函数在图 3-6 中给

7、出。2) 第 II 区发生在L1LzNviNvg50 36Nvi的范围内。此时其中Fs和 F6作为 N1的函数在课本 P38图 3-7 图给出。3) 2)第 III 区0 75发生在Nvg75 84Nvl的范围内。此时4、算法框图已知参数P猜测参数F2(0)NS(1+F5)卄(1 F7NQ3、NSP的计算v Ns(vPf(i)Ph(i)=(P(i)(0.982vgF6N(1+F5)(T F7Nvl)2(9.81/Vsl Vsg) !FrPh(i)Vsl2+ (1x 0.25l)(v Vsl Vsg)2I(1 Vsg/Vsl)-)I)9.8pf(i)/(1 (1Vsl1/24 v Vsg) )/(2 v)z/dgVsg) Vsg)/P(i)/1000000)/10685、设计程序源代码（见附页程序二） 四、程序代码1、Beggs-Brill法源代码2、Duns-Ros法源代码参考文献：1、徐春碧等.气液两相管流理论与计算，重庆科技学院，2006.12,陈家朗石油