计算井筒多相流的BeggsBrill方法.docx

1、计算井筒多相流的BeggsBrill方法作业一：计算井筒多相流的Beggs-Brill方法基本参数：套管尺寸：5 1/2, 油管：2 1/2, GOR=24/, =0.9659, =100cp, 井底温度63，井口温度25，井口压力1MPa，产液137方/天，油层深度1490m，井底流压8MPa，含水率29，泵的单级扬程4.6m。设计步骤：截图：Sub BeggsBrill(cal As Single, rw As Single, ro As Single, rg As Single, yco2 As Single, yh2s As Single, yn2 As Single, ynacl A

2、s Single, _ppc As Single, tpc As Single, d2 As Single, d1 As Single, eras As Single, lon As Single, alfa0 As Single, T1 As Single, pcal1 As Single, T2 As Single, qw As Single, qo As Single, Qsc As Single, daltal As Single, PCAL2 As Single)Do If lsum + daltal lon Thendaltal = lon - lsum End Iflmid =

3、lsum + daltal / 2 T = T1 + (T2 - T1) * lmid / lon Do daltap0 = daltap p = p1 + Sgn(cal) * daltap / 2 p利用假设的daltap计算，带入以后的plsum = lsum + daltal p2 = p1 + Sgn(cal) * daltap p1 = p2 Loop Until lsum = lon PCAL2 = p2 End Subpsc = 0.101325tsc = 293densc = 1.204 * rgdeno = 1000 * rodenw = 1000 * rwqow = qo

4、 + qw If qo 0 Thengor = Qsc / qo Call PVTO(ro, rg, gor, T, p, pbo, rso, Bo, uo, sto) End If If qw 0 Thengwr = Qsc / qw Call PVTW(T, p, ynacl, gwr, rsw, bw, uw, stw) End If If qow 0 Then GLR = Qsc / qow RS = (rso * qo + rsw * qw) / qow If RS = GLR Then RS = GLR End Ifqgsc = qow * (GLR - RS) Elseqgsc

5、= Qsc End If If qgsc 0 Thenppr = p / ppctpr = T / tpc Call dprzf(ppr, tpr, dr, z, cpr)bg = psc / tsc * z * T / p Call lgeug(rg, yn2, yco2, yh2s, p, T, z, ug)qg = qgsc * bgwg = densc * qgscdensg = wg / qg Elseqg = 0densg = 0bg = 0 End If Ql = qo * Bo + qw * bwwl = deno * qo + denw * qw + densc * RS *

6、 qow If Ql 0 Thendensl = wl / Qlfo = qo * Bo / Ql End Iful = uo * fo + uw * (1 - fo)stal = sto * fo + stw * (1 - fo)End Subgor = gorin * 5.615 T = tin * 1.8 - 460 p = 145 * pinapi = 141.5 / ro - 131.5 If api= pbo Or rso = gor Thenrso = gor End If If api= pbo Or rso = gor Then Bo = Bo * Exp(co * (pbo

7、 - p)zz = 3.0324 - 0.0203 * apiyy = 10 zzxx = yy * T (-1.163)uod = 10 xx - 1aa = 10.715 * (rso + 100) (-0.515)bb = 5.44 * (rso + 150) (-0.338) uo = aa * uod bbuo = 100 If pbo 0 And (p = pbo Or rso = gor) Thenuo = uo * (p / pbo) (2.6 * p 1.187 * Exp(-0.0000898 * p - 11.513) End Ifsto = (42.4 - 0.047

8、* T - 0.267 * api) * Exp(-0.0007 * p) / 1000pbo = pbo / 145rso = rso / 5.615End Sub T = tin * 1.8 - 460 F p = 145 * pin a1 = 2.12 + 0.00345 * T - 0.0000359 * T 2 b1 = 0.0107 - 0.0000526 * T + 1.48 * 10 (-7) * T 2 c1 = -8.75 * 10 (-7) + 3.9 * 10 (-9) * T - 1.02 * 10 (-11) * T 2scnacl = 1 - (0.0753 -

9、0.000173 * T) * ynaclrws = (a1 + b1 * p + c1 * p 2) * scnacl / 5.615 If rws gwr Thenrws = gwr End If If rws 0 Thennv = Vsl * (densl / (ga * stal) 0.25 End Ifdenn = na * densl + (1 - na) * densg un = (na * ul + (1 - na) * ug) / 1000 ren = denn * Vm * don / un fmn = (1 / (1.14 - 2 * Log(eras / don + 2

10、1.25 / ren 0.9) / Log(10) 2 End Sub If na 0 Then ls1 = 316 * na 0.302 ls2 = 0.0009252 * na (-2.4684) ls3 = 0.1 * na (-1.4516) ls4 = 0.5 * na (-6.738) If na 0.99999 Then ident = 6 ElseIf na 0.01 And nr = 0.01 And nr = 0.01 And na ls3 And nr = 0.4 And (nr ls3 And nr = ls4) Then ident = 2 ElseIf na= ls

11、1 Or na = 0.4 And nr = ls4 Then ident = 3 ElseIf na= 0.01 And (nr = ls2 And nr = 0 Then Select Case ident Case 1dd = 0.011ee = -3.768ff = 3.539gg = -1.614 Case 2dd = 2.96ee = 0.305ff = -0.4473gg = 0.0978 Case 3dd = 1ee = 0ff = 0gg = 0 End Select Elsedd = 4.7ee = -0.3692ff = 0.1244gg = -0.5056 End If

12、 cc = (1 - na) * Log(dd * na ee * nv ff * nr gg) If cc 0 Then cc = 0 End Iftm = Sin(1.8 * alfa)hui = 1 + cc * (tm - tm 3 / 3) If hui 0 Thenhui = 0 End If Select Case ident Case 1dd = 0.98ee = 0.4846 cc = 0.0868 Case 2dd = 0.845ee = 0.5351 cc = 0.0173 Case 3dd = 1.065ee = 0.5824 cc = 0.0609 End Selec

13、t hl0 = dd * na ee * nr (-cc) If hl0 na Then hl0 = na End Ifhlsta = 0.918 * hl0 * hui If hlsta 1 Thenhlsta = 1 End IfEnd Subga = 9.807 Select Case ident Case 4 ls2 = 0.0009252 * na (-2.4684) ls3 = 0.1 * na (-1.4516)ident = 1 Call HOLDUP(alfa, ident, na, nv, nr, hlsta) hlsta1 = hlstaident = 2 Call HO

14、LDUP(alfa, ident, na, nv, nr, hlsta) hlsta2 = hlstatm = (ls3 - nr) / (ls3 - ls2)hlsta = tm * hlsta1 + (1 - tm) * hlsta2 Case 5hlsta = 0 Case 6hlsta = 1 Case Else Call HOLDUP(alfa, ident, na, nv, nr, hlsta) End Select If hlsta 1 And yb 1.2 Thensb = Log(2.2 * yb - 1.2) Elsesb = yr / (-0.0523 + 3.182 *

15、 yr - 0.8725 * yr 2 + 0.01853 * yr 4) =16 End If End Iffm = fmn * Exp(sb) Tf = fm * denn * Vm 2 / (2 * don) coef = 1 - denm * Vm * vsg / (1000000 * p) T1 = (tm + Tf) / (1000000 * coef) daltap = daltal * T1 End Sub a1 = 0.31506237 a2 = -1.0467099 a3 = -0.57832729 a4 = 0.53530771 a5 = -0.61232032 a6 =

16、 -0.10488813 a7 = 0.681570001 a8 = 0.68446549 a = a1 + a2 / tr + a3 / tr 3 B = a4 + a5 / tr C = a5 * a6 * tr D = a7 / tr 3 e = a8 G = -0.27 * pr / trdr = -G Doed = e * dr 2exa = Exp(-ed) f = G / dr + 1 + a * dr + B * dr 2 + C * dr 5 f = f + D * dr 2 * (1 + ed) * exafp = -G / dr 2 + a + 2 * B * dr +

17、5 * C * dr 4fp = fp + 2 * D * dr * (1 + ed - ed 2) * exadr = dr - f / fp If dr = 0 Or dr = 1 Then Exit Sub End If Loop Until Abs(f / fp) 0.00001 z = 0.27 * pr / (dr * tr)ed = e * dr 2exa = Exp(-ed)dzdr = a + 2 * B * dr + 5 * C * dr 4dzdr = dzdr + 2 * D * dr * (1 + ed - ed 2) * exacr = 1 / pr - 0.27

18、/ (z 2 * tr) * dzdr / (1 + dr * dzdr / z)End Submg = 28.97 * rg R = 0.008314 g1 = yh2s * (5.7 * rg - 1.7) * 0.00001 g2 = yco2 * (5 * rg + 1.7) * 0.00001 g3 = yn2 * (5 * rg + 4.7) * 0.00001tr = 1.8 * T G = (9.399999 + 0.02 * mg) * tr 1.5 G = 0.0001 * G / (209 + 19 * mg + tr) + g1 + g2 + g3 j = 1.7 -

19、197.2 / tr - 0.002 * mg h = (3.5 + 986 / tr + 0.01 * mg) * 0.001 jden = mg * p / (z * R * T)ug = G * Exp(h * den j)End Sub作业二抽油杆柱设计基本参数配产：30m3/d,泵深：2300，动液面：950m，冲程：3.6，冲次6次/min，泵径：56mm，油管内径：73mm，地层水粘度：1,8mPa.s，地层水密度：1050kg/m3，抽油杆密度：7850kg/m3，最小杆径：19mm，杆使用系数：0.95设计步骤：截图：Dim hc As Single 沉没度，mDim SF

20、As Single 杆使用系数Dim s As Single 冲程，mDim ql As Single 产液量，m3/dDim dp As Single 泵径，mmDim de As Single 柱塞与泵筒间隙平均值，mDim uw As Single 地层水粘度，mPa.sDim n As Single 冲次，次/minDim rouw As Single 地层水密度，kg/m3Dim rous As Single 钢的密度，kg/m3Dim t As Single 抽油杆最小抗张强度，MPaDim Pl As Single, PLi As Single, deltaPL As Singl

21、e 输入应力范围比、计算应力范围比、两级抽油杆顶端面应力范围比的最大允许差值Dim dt As Single 油管内径，mmDim fp As Single, fr As Single, fo As Single 柱塞截面积，抽油杆截面积，阀孔截面积，m2Dim nrod As Integer 杆级数Dim L0 As Single, L1 As Single 计算所得杆柱长度，m 已设计出的杆柱总长度，mDim dfo As Single 阀孔径，mDim qr() As Single 每米抽油杆质量，kg/mDim Li() As Single 每级杆长度，mDim dr() As Sin

22、gle 抽油杆直径,mmDim b As Single 失重系数Dim m As Single 油管内径与抽油杆直径之比Dim vmax As Single 抽油杆柱最大下行速度，m/sDim vp As Single, vf As Single 柱塞最大运动速度、阀孔液流速度，m/sDim NRe As Single, ksi As Single 雷诺数，阀流量系数Dim Pmax0 As Single, Pmin0 As Single, PL0 As Single 抽油杆上下冲程底端载荷，N，底端应力范围比Dim () As Single, dL() As Single, d() As Single, d() As Single 井斜角，井斜深、井斜角、方位角变化量Dim Pmaxi1() As Single, Pmini1() As Single, PLi1() As Single 第一级杆节点最大、最小载荷，N, 最大应力范围比定义动态数组Dim Pmaxi2() As Single, Pmini2() As Single, PL