- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
e 1.7315 , f 1.5459 , g 0.7581
H1 (90 ) H1 (0) aE1b N Fr
c
1.17
二、气液两相总压降和摩阻压降计算 (4)计算阻力系数
y
其中：
H (90 )
1
E1
2
0.73
S
ln y 0.28 2 4 0.0523 3.182 ln y 0.8725 (ln y ) 0.01853 (ln y )
总结（上下段） 绘制比例与产量的关系如下：
程序文本
#include<stdio.h> void main() { float pg,pl,po,Yg,pa,P,To,Z,T,Po,Rs,Rp,Qo,Qg,Ql,Bo,Vm,Gm; scanf("%f%f%f%f%f%f%f%f%f%f%f%f",&po,&pa,&Rs,&Yg,&Bo,&P,&To,&Z,&T, &Po,&Rp,&Qo); pl=(po+pa*Rs*Yg)/Bo; pg=(Yg*pa*P*To)/(Z*T*Po); Qg=(Po*T*Z*Qo)*(Rp-Rs)/(86400*P*To); Ql=Qo*Bo/86400; Vm=Qg/0.003+Ql/0.003; Gm=pl*Ql+pg*Qg; printf("%f %f %f %f %f %f\n",Qg,Ql,Vm,Gm,pl,pg); return 0; }
至井口的压力损失为10-0.1=9.9mPa。摩擦损失梯度
dP v 2 21 .49 pa / m ，摩擦损失为26642Pa. 为 dh 2D
总结（上下段） 由Beggs-Brill方法按压力增量迭代法通 过编程计算得到的结果如下表：
数据 井筒摩阻损 失/pa 1 2 3 4 3709.09 16984.21 23967.97 30792.19 压力损失 /mpa 12.45 12.40 12.35 12.28 摩阻损失梯 压力损失梯 度pa/m 2.47 11.32 15.98 20.53 度mpa/m 8300.00 8266.67 8230.00 8186.67 0.03% 0.14% 0.19% 0.25% 梯度比
程序文本
#include<stdio.h> #include <math.h> void main() { float a,b,c,d,e,f,g,El,Nfr,Nvl,H1,H2,C,five; scanf("%f%f%f%f%f%f%f%f%f%f",&a,&b,&c,&d,&e,&f,&g,&El,&Nvl,&Nfr); H2=a*pow(El,b)/pow(Nfr,c); C=(1-El)*log(d*pow(El,e)*pow(Nvl,f)*pow(Nfr,g)); five=1+0.3*C; H1=H2*five; printf("%f %f \n",H1,C); return 0; }
二、气液两相总压降和摩阻压降计算
② 1190m至井口的压力损失为10-0.2=9.8mPa。摩擦损失 dP v 2 17 .06 pa / m ，摩擦损失为20304Pa. 梯度为
dh 2D
③ Pwf 12 .1Mpa ，Q 59 .8m 3 / d 时，井底压力为12.1mPa， 井口假设为0.1mpa，所以总压力损失为12mPa。1240
二、气液两相总压降和摩阻压降计算
• 原油密度: • 天然气密度: • 气液就地流量:
l o a Rs g
BO 813 .89 kg / m3
g
g a PTO
ZTPO
48 .48 kg / m 3
qg=0
• 就地的气体、液体及混合物的表观流速：
D2 Ap 0.003 m 2 4
dP P计 Z 8.36 Mpa dz
(6)计算该段末端的压力
P t 10 8.36 1.64Mpa
二、气液两相总压降和摩阻压降计算
，
(7)取 0.4Mpa
Pt 1.4 0.24 0.4
所以假设值合理。
(8)井下压力为14.1mPa，井口压力为1.4mPa，所以总
vsg 0
vm vsl vsg 0.072 m / s
ql vsl 0.072 m / s Ap
二、气液两相总压降和摩阻压降计算
•液、气及总的质量流量
Gl l ql 0.176
Gg 0
Gm Gl Gg 0.176
•入口体积含液率（无滑脱持液率）
ql El 1 ql q g
程序文本
#include<stdio.h> #include <math.h> void main() { float Vm,pl,El,pg,ndl,ndg,Nre,rmna1,rmna,S; scanf("%f%f%f%f%f%f%f",&Vm,&pl,&pg,&El,&ndl,&ndg,&S); Nre=0.062*Vm*(pl*El+pg*(1-El))/(ndl056+0.5/pow(Nre,0.32); rmna=rmna1*pow(2.718,S); printf("%f\n",rmna); return 0; }
流量（m3/d） 17.5 41.6 50.4 59.8
流压（MPa）
14.1
12.9
12.5
12.1
题目一
为便于在后续的设计计算中简化进行井筒压力损失 计算，试设定一般生产条件，通过不同产量下的井 筒摩阻及总压力损失的计算，确定井筒摩擦压力梯 度占总压力梯度的大约比例，以及该比例与产量的 大致关系。
取纯油流的原油密度等于油藏条件下原油的密度，该
段的平均密度：
重力压力梯度为：
混合流体粘度的计算：
雷诺数的计算：
， 属于
一、纯液柱段摩阻压降和总压降计算 （3）摩擦损失梯度的计算
摩擦损失梯度为：
（4）纯油流高度的计算
油层压力： 纯油流压力降： 纯油流高度：
（5）摩阻压降计算
二、气液两相总压降和摩阻压降计算
程序文本
#include<stdio.h> void main() { float pg,pl,H1,rmna,Gm,Vm,Vsg,P,gradient,V1,V2,V3; scanf("%f%f%f%f%f%f%f%f",&pl,&pg,&H1,&rmna,&Gm,&Vm,&Vsg,&P); V1=(pl*H1+pg*(1-H1))*9.81; V2=rmna*Gm*Vm/(2*0.062*0.003); V3=(pl*H1+pg*(1-H1))*Vm*Vsg/P; gradient=0-(V1+V2)/(1-V3); printf("%f\n",gradient); return 0; }
'es 0.084
其中， 0.0056
'
0.5 ( N Re ' ) 0.32
二、气液两相总压降和摩阻压降计算
和
(5)利用公式计算压力梯度
dP dz
和压力降 P计
m vm 1H1 g (1 H1 ) g 2G DA p dP 9.37 10 3 Mpa / m dz 1 1H1 g (1 H1 )vm vsg / P
生产条件：直井自喷生产
一、纯液柱段摩阻压降和总压降计算
（1）油藏条件下溶解气油比的计算：
原油的API度：
地面脱气原油的有效分子量：
泡点压力系数：
天然气的摩尔分数： 使用Lastater公式计算Rs:
由于溶解气油比大于21.25,取
一、纯液柱段摩阻压降和总压降计算 （2）雷诺数的计算
油藏条件下原油的密度为：
压力损失为12.7mPa。892m处至井口的压力损失为 10-1.4=8.6mPa。
dP v 2 摩擦损失梯度为 2.90 pa / m ，摩擦损失为 dh 2D 2583pa.
度为
二、气液两相总压降和摩阻压降计算
，
(9)同理
3 ① Pwf 12 .9 Mpa , Q 41 .6m / d 时,井底压力为12.9mPa，
•弗鲁德数
N Fr
vm 2 0.009 gD
l N vl vsl g
0.25
•液体速度准数
0.6
二、气液两相总压降和摩阻压降计算 (2)确定流型
L1 316 E10.302 316
L2 92 .52 10 5 E12.4684 0.001
（Beggs-Brill 方法）
已知：Pwf=14.1MPa Q=17.5m3/d
（1）确定起始点压力
P1=10MPa 计算深度增量 △Z=(1500-608)=892 分段数 N=1 假设井口压力为1.4MPa 计算 P 5.7 mpa , T 307.79K 下的流体性质 参数及流动参数,由平均温度和压力计算得 Rs=21.25 Bo=1.068,Z=0.933 μ l=2.733mPa.s μ g=0.0122mPa.s, σ =0.0167 ρ α =1.293
课题名称：地层及井筒流动第一题
基础数据资料
