tqb时，令qt110td，则pwf1同理，qt220td，Pwf213877Mpaqt330td，Pwf3120Mpa当qbqtqomzx时，令qt450td则按流压加权平均进行推导得：
wP1
P1
qtj15754Mpa
Pwf4f
qqbqt8180tqomzxqbj01251fwPb18166Mpa
4
f同理qt560td，Pwf55860Mpa
q时，
t
当qomzx
pwffwp1
qomzxqqomzx8fw9tJJ
令qt671td，Pwf62233Mpa综上，井底流压与产量的关系列表如下：PwfMpaQtd157471013873201203010040653816650586060223371
得到油井的流入动态曲线如下图：
IPR曲线20
PwfMpa
15105001020304050Q（td）60708090
图1
油井IPR曲线
22井筒多相流的计算
井筒多相流压力梯度方程井筒多相管流的压力梯度包括：因举高液体而克服重力所需的压力势能、流体因加速而增加的动能和流体沿管路的摩阻损失，其数学表达式如下：
dpdhρ
dvmfmmgsi
θρmvmdhρ
m
vmmd2
2
式中ρ
为多相混合物的密度vm为多相混合物的流速fm为多相混合物流
5
f动时的摩擦阻力系数d为管径p为压力h为深度g为重力加速度θ为井斜角的余角。井筒多相管流计算包括两部分：（1）由井底向上计算至泵入口处；（2）油管内由井口向下计算至泵出口处。1）由井底向上计算至泵入口处，计算下泵深度Lp。采用深度增量迭代方法，首先估算迭代深度。在本设计中为了减小工作量，采用只迭代一次的方法。计算井筒多相管流时，首先计算井筒温度场、流体物性参数，然后利用Orkiszewski方法判断流型，进行压力梯度计算，最后计算出深度增量和下泵深度Lp。按深度增量迭代的步骤井底流压12Mpa，假设压力降为02Mpa；估计一个对应的深度增量h40m，即深度为1960m。由井温关系式可以计算得到该处的井温为：8996℃。平均的压力和温度：T（908996）28998℃。平均压力P119Mpa。由平均压力和平均温度计算的得到流体的物性参数为：溶解油气比RS7131；原油体积系数B012584340；死油粘度μ
4
原油密度P073900；油水混合液的密度Pz
4
od
653710；活油粘度μ
4
O
4
331810；
水的粘
度μ
w
326310；
液体的粘度μ329610；天然气的压缩因子Z09567；
天然气的密度
g
9070。以上单位均是标准单位。
由以上的流体物性参数判断流型：不同流动型态下的m和f的计算方法不同，为此，计算中首先要判断流动形态。该方法的四种流动型态的划分界限如表1所示。表1流型界限流动型态泡流界
qgLBqtqgLBvgLSqt
限
段r