课题研究答辩课件

课题名称：地层及井筒流动题课题名称：地层及井筒流动题油层厚度—5m中深—1500m

地层压力系数—1.15

地温梯度--2.90C/100m

原油饱和压力—10.0MPa含水率为—10%

生产气油比—25m3/t原油相对密度—0.85

天然气相对密度—0.7

500C脱气原油粘度—30mPa.s

该井射孔完井，孔密24孔/m，孔径14mm，油管内径62mm。该井试油资料如下:基础数据资料油层厚度—5m中深—1500m题目一题目一的任务为便于在后续的设计计算中简化进行井筒压力损失计算，试设定一般生产条件，通过不同产量下的井筒摩阻及总压力损失的计算，确定井筒摩擦压力梯度占总压力梯度的大约比例，以及该比例与产量的大致关系。生产条件：直井自喷生产题目一题目一的任务为便于在后续的设计计算中简化进行井筒压力损一、纯液柱段摩阻压降和总压降计算（1）油藏条件下溶解气油比的计算：原油的API度：地面脱气原油的有效分子量：泡点压力系数：天然气的摩尔分数：使用Lastater公式计算Rs:

由于溶解气油比大于21.25,取以第一组数据为例一、纯液柱段摩阻压降和总压降计算（1）油藏条件下溶解气油比的（2）雷诺数的计算油藏条件下原油的密度为：取纯油流的原油密度等于油藏条件下原油的密度，该段的平均密度：重力压力梯度为：混合流体粘度的计算：雷诺数的计算：

，

属于一、纯液柱段摩阻压降和总压降计算层流（2）雷诺数的计算一、纯液柱段摩阻压降和总压降计算层流（3）摩擦损失梯度的计算

一、纯液柱段摩阻压降和总压降计算（5）摩阻压降计算

（4）纯油流高度的计算

油层压力：

纯油流压力降：

纯油流高度：摩擦损失梯度为：

（3）摩擦损失梯度的计算一、纯液柱段摩阻压降和总压降计算（5二、气液两相总压降和摩阻压降计算以第一组数据为例已知：Pwf=14.1MPaQ=17.5m3/d（1）确定起始点压力

P1=10MPa

计算深度增量

△Z=(1500-608)=892

分段数N=1

假设井口压力为1.4MPa

计算,下的流体性质

参数及流动参数,由平均温度和压力计算得Rs=21.25Bo=1.068,Z=0.933

μl=2.733mPa.s

μg=0.0122mPa.s,

σ=0.0167

ρα=1.293

（Beggs-Brill方法）二、气液两相总压降和摩阻压降计算以第一组数据为例已知：Pwf二、气液两相总压降和摩阻压降计算原油密度:

天然气密度:

气液就地流量:就地的气体、液体及混合物的表观流速：

qg=0

二、气液两相总压降和摩阻压降计算原油密度:天然气密度:气液就二、气液两相总压降和摩阻压降计算•液、气及总的质量流量

•入口体积含液率（无滑脱持液率）

•弗鲁德数

•液体速度准数

二、气液两相总压降和摩阻压降计算•液、气及总的质量流量二、气液两相总压降和摩阻压降计算(2)确定流型因为

所以判断为过渡流。

二、气液两相总压降和摩阻压降计算(2)确定流型二、气液两相总压降和摩阻压降计算（3）根据流型计算

由于流型为过渡流，查表1-6和1-7得：

二、气液两相总压降和摩阻压降计算（3）根据流型计算二、气液两相总压降和摩阻压降计算(4)计算阻力系数

其中：其中，二、气液两相总压降和摩阻压降计算(4)计算阻力系数其二、气液两相总压降和摩阻压降计算(5)利用公式计算压力梯度

和压力降和

(6)计算该段末端的压力

二、气液两相总压降和摩阻压降计算(5)利用公式计算压力梯度二、气液两相总压降和摩阻压降计算(7)取所以假设值合理。(8)井下压力为14.1mPa，井口压力为1.4mPa，所以总压力损失为12.7mPa。892m处至井口的压力损失为10-1.4=8.6mPa。摩擦损失梯度为

，摩擦损失为2583pa.，

二、气液两相总压降和摩阻压降计算(7)取所以假设值合理。(8二、气液两相总压降和摩阻压降计算(9)同理①,

时,井底压力为12.9mPa，井口假设为0.3mPa，所以总压力损失为12.6mPa。1141m至井口的压力损失为10-0.3=9.7mPa。摩擦损失梯度为，摩擦损失为13945Pa.度为

,

时,井底压力为12.5mPa，井口假设为0.2mPa，所以总压力损失为12.3mPa。，

二、气液两相总压降和摩阻压降计算(9)同理①二、气液两相总压降和摩阻压降计算，

时，井底压力为12.1mPa，井口假设为0.1mpa，所以总压力损失为12mPa。1240至井口的压力损失为10-0.1=9.9mPa。摩擦损失梯度为，摩擦损失为26642Pa.1190m至井口的压力损失为10-0.2=9.8mPa。摩擦损失梯度为

，摩擦损失为20304Pa.二、气液两相总压降和摩阻压降计算，总结（上下段）由Beggs-Brill方法按压力增量迭代法通过编程计算得到的结果如下表：总结（上下段）由Beggs-Brill方法按压力增量迭代法通总结（上下段）绘制比例与产量的关系如下：总结（上下段）绘制比例与产量的关系如下：程序文本#include<stdio.h>voidmain(){ floatpg,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);return0;}程序文本#include<stdio.h>程序文本#include<stdio.h>#include<math.h>voidmain(){ floatpg,pl,F,Yg,pa,P,To,Z,T,Po,Rs,Rp,Qo,Qg,Ql,Bo,Ap,Vsl,Vsg,Vm,El,Nfr,Nvl,L1,L2,L3,L4; scanf("%f%f%f%f%f%f%f%f%f%f%f",&Po,&T,&Z,&Rp,&Rs,&Qo,&P,&To,&Bo,&pl,&F); Qg=(Po*T*Z*Qo)*(Rp-Rs)/(86400*P*To); Ql=Qo*Bo/86400; Vsl=Ql/0.003; Vsg=Qg/0.003; Vm=Vsl+Vsg; El=Ql/(Ql+Qg); Nvl=Vsl*pow(pl/(9.81*F),0.25); Nfr=Vm*Vm/(9.81*0.062); L1=316*pow(El,0.302);L2=0.0009252*pow(El,-2.4684); L3=0.1*pow(El,-1.4516);L4=0.5*pow(El,-6.738); printf("%f%f%f%f%f%f%f\n",El,Nfr,Nvl,L1,L2,L3,L4); return0;}程序文本#include<stdio.h>程序文本#include<stdio.h>#include<math.h>voidmain(){ floata,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);return0;}程序文本#include<stdio.h>程序文本#include<stdio.h>#include<math.h>voidmain(){ floatVm,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))/(ndl*El+ndg*(1-El)); rmna1=0.0056+0.5/pow(Nre,0.32); rmna=rmna1*pow(2.718,S); printf("%f\n",rmna);return0;}程序文本#include<stdio.h>程序文本#include<stdio.h>voidmain(){ floatpg,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);return0;}程序文本#include<stdio.h>课题研究答辩课件课题名称：地层及井筒流动题课题名称：地层及井筒流动题油层厚度—5m中深—1500m

地层压力系数—1.15

地温梯度--2.90C/100m

原油饱和压力—10.0MPa含水率为—10%

生产气油比—25m3/t原油相对密度—0.85

天然气相对密度—0.7

500C脱气原油粘度—30mPa.s

该井射孔完井，孔密24孔/m，孔径14mm，油管内径62mm。该井试油资料如下:基础数据资料油层厚度—5m中深—1500m题目一题目一的任务为便于在后续的设计计算中简化进行井筒压力损失计算，试设定一般生产条件，通过不同产量下的井筒摩阻及总压力损失的计算，确定井筒摩擦压力梯度占总压力梯度的大约比例，以及该比例与产量的大致关系。生产条件：直井自喷生产题目一题目一的任务为便于在后续的设计计算中简化进行井筒压力损一、纯液柱段摩阻压降和总压降计算（1）油藏条件下溶解气油比的计算：原油的API度：地面脱气原油的有效分子量：泡点压力系数：天然气的摩尔分数：使用Lastater公式计算Rs:

由于溶解气油比大于21.25,取以第一组数据为例一、纯液柱段摩阻压降和总压降计算（1）油藏条件下溶解气油比的（2）雷诺数的计算油藏条件下原油的密度为：取纯油流的原油密度等于油藏条件下原油的密度，该段的平均密度：重力压力梯度为：混合流体粘度的计算：雷诺数的计算：

，

属于一、纯液柱段摩阻压降和总压降计算层流（2）雷诺数的计算一、纯液柱段摩阻压降和总压降计算层流（3）摩擦损失梯度的计算

一、纯液柱段摩阻压降和总压降计算（5）摩阻压降计算

（4）纯油流高度的计算

油层压力：

纯油流压力降：

纯油流高度：摩擦损失梯度为：

（3）摩擦损失梯度的计算一、纯液柱段摩阻压降和总压降计算（5二、气液两相总压降和摩阻压降计算以第一组数据为例已知：Pwf=14.1MPaQ=17.5m3/d（1）确定起始点压力

P1=10MPa

计算深度增量

△Z=(1500-608)=892

分段数N=1

假设井口压力为1.4MPa

计算,下的流体性质

参数及流动参数,由平均温度和压力计算得Rs=21.25Bo=1.068,Z=0.933

μl=2.733mPa.s

μg=0.0122mPa.s,

σ=0.0167

ρα=1.293

（Beggs-Brill方法）二、气液两相总压降和摩阻压降计算以第一组数据为例已知：Pwf二、气液两相总压降和摩阻压降计算原油密度:

天然气密度:

气液就地流量:就地的气体、液体及混合物的表观流速：

qg=0

二、气液两相总压降和摩阻压降计算原油密度:天然气密度:气液就二、气液两相总压降和摩阻压降计算•液、气及总的质量流量

•入口体积含液率（无滑脱持液率）

•弗鲁德数

•液体速度准数

二、气液两相总压降和摩阻压降计算•液、气及总的质量流量二、气液两相总压降和摩阻压降计算(2)确定流型因为

所以判断为过渡流。

二、气液两相总压降和摩阻压降计算(2)确定流型二、气液两相总压降和摩阻压降计算（3）根据流型计算

由于流型为过渡流，查表1-6和1-7得：

二、气液两相总压降和摩阻压降计算（3）根据流型计算二、气液两相总压降和摩阻压降计算(4)计算阻力系数

其中：其中，二、气液两相总压降和摩阻压降计算(4)计算阻力系数其二、气液两相总压降和摩阻压降计算(5)利用公式计算压力梯度

和压力降和

(6)计算该段末端的压力

二、气液两相总压降和摩阻压降计算(5)利用公式计算压力梯度二、气液两相总压降和摩阻压降计算(7)取所以假设值合理。(8)井下压力为14.1mPa，井口压力为1.4mPa，所以总压力损失为12.7mPa。892m处至井口的压力损失为10-1.4=8.6mPa。摩擦损失梯度为

，摩擦损失为2583pa.，

二、气液两相总压降和摩阻压降计算(7)取所以假设值合理。(8二、气液两相总压降和摩阻压降计算(9)同理①,

时,井底压力为12.9mPa，井口假设为0.3mPa，所以总压力损失为12.6mPa。1141m至井口的压力损失为10-0.3=9.7mPa。摩擦损失梯度为，摩擦损失为13945Pa.度为

,

时,井底压力为12.5mPa，井口假设为0.2mPa，所以总压力损失为12.3mPa。，

二、气液两相总压降和摩阻压降计算(9)同理①二、气液两相总压降和摩阻压降计算，

时，井底压力为12.1mPa，井口假设为0.1mpa，所以总压力损失为12mPa。1240至井口的压力损失为10-0.1=9.9mPa。摩擦损失梯度为，摩擦损失为26642Pa.1190m至井口的压力损失为10-0.2=9.8mPa。摩擦损失梯度为

，摩擦损失为20304Pa.二、气液两相总压降和摩阻压降计算，总结（上下段）由Beggs-Brill方法按压力增量迭代法通过编程计算得到的结果如下表：总结（上下段）由Beggs-Brill方法按压力增量迭代法通总结（上下段）绘制比例与产量的关系如下：总结（上下段）绘制比例与产量的关系如下：程序文本#include<stdio.h>voidmain(){ floatpg,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);return0;}程序文本#include<stdio.h>程序文本#include<stdio.h>#include<math.h>voidmain(){ floatpg,pl,F,Yg,pa,P,To,Z,T,Po,Rs,Rp,Qo,Qg,Ql,Bo,Ap,Vsl,Vsg,Vm,El,Nfr,Nvl,L1,L2,L3,L4; scanf("%f%f%f%f%f%f%f%f%f%f%f",&Po,&T,&Z,&Rp,&Rs,&Qo,&P,&To,&Bo,&pl,&F); Qg=(Po*T*Z*Qo)*(Rp-Rs)/(86400*P*To); Ql=Qo*Bo/86400; Vsl=Ql/0.003; Vsg=Qg/0.003; Vm=Vsl+Vsg; El=Ql/(Ql+Qg); Nvl=Vsl*pow(pl/(9.81*F),0.25); Nfr=Vm*Vm/(9.81*0.062); L1=316*pow(El,0.302);L2=0.0009252*pow(El,-2.4684); L3=0.1*pow(El,-1.4516);L4=0.5*pow(El,-6.738); printf("%f%f%f%f%f%f%f\n",El,Nfr,Nvl,L1,L2,L3,L4); ret