计算井筒多相流的Beggs-Brill方法.docxVIP

作业一：计算井筒多相流的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 As 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 = 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 + 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 = 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 * 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 = 30 Then