《采油工程》考试改革课堂考核环节试题第一次采油大作业答案

第一题

生产初期假设该井可以自喷生产，井筒中的流动可以分为两段。下部分泡点压力以下为纯液流，上端低于泡点压力之后为气液两相流。忽略加速度压力梯度部分。

为了简化计算，大概确定摩阻压力梯度的比例，讲井筒管流分为两部分，纯液流和气液两相流。以第一组数据为例，根据混合液的密度可以得到液柱高600.92m，气液混合物高度899.08m。

（1） 纯液柱段摩阻压降和总压降计算： 油藏条件下的原油密度：

?oi??o??aRs?gBo

根据油层物理第一章的内容，我们可以得到油藏条件下的溶解气油比

Rs?25.94m3/m3

于是可得油藏条件下的原油密度：

?oi?831.88kg/m3

原油析出气体前可忽略压力所引起的密度变化，因此该段原油密度可近似取原始条件下的原油密度。于是该段的平均密度：

?m??oi(1?10%)??w?10%?848.7kg/m3

重力压力梯度：

dph/dh??m?g

原油流速：

qm?(qoBo?qwBw)/86400

雷诺数：

NRe?Dvm?m?m

其中粘度为油水的体积加权平均值，原油的粘度根据油层物理学中相关公式得到。 根据雷诺数的大小，所给四组生产条件下的流动皆为水力光滑区。故有：

f?摩擦损失梯度：

0.316414NRe

dpf/dh?f2?mvmD2

（2） 气液共存段摩阻压降及总压降计算： 为简化运算，气液共存段不分段，使用Orkiszewski方法进行计算。该段平均压力为：（9+0.1）/2=4.505MPa，按照温度梯度计算中点温度值作为该段平均温度。

气，液的就地流量：

qg?p0TZ(Rp?Rs)qo86400pT0qoBo?qwBw86400

ql?按照Orkiszewski方法流型划分原则，分别计算LB（泡流界限），LS（段塞流界限）等等，并分析不同产量下的流型。四种产量下的流型都为泡流，其摩擦损失梯度可以按照液相进行计算。但是液相流速需要计算持气率Hg，泡流持气率计算参见《采油工程原理与设计》。 液相真实流速：

vlh?qlAp(1?Hg)

按照相同的方法可得到摩擦阻力系数，从而计算该段的摩擦阻力损失。 根据持气率亦可计算该段平均密度，从而得到气液共存段的重力压力梯度。

（3）根据各段的压力梯度和高度，可分别计算总的重力压差和摩阻压差，从而得到摩阻压力梯度所占的比例。结果见附表 产液 11 19.8 25.4 31.3

摩阻压降/Pa 299.4786511 1265.845917 2378.91156 4039.606879 重力压降/Pa 摩阻压降占比/% 12635962.98 0.0024 12444729.22 0.0102 12331780.59 0.0193 12218334.34 0.0331 摩阻与产量关系0.03500.03000.0250摩阻占比0.02000.01500.01000.00500.0000051015产液量

20253035第二题

一、 计算地层压力：P=1.15x1000x9.8x1500=16.905MPa；选择合适的产量和压力(可自己选择):流压(12.9MPa)—流量(19.8m/d)。

二、 根据题目相应地层和流体参数计算某一油管直径的井筒压力梯度损失。（可利用matlab编程来实现）

已知参数：油层厚度5m；中深1500m；地层压力系数1.15；原油饱和压力9.0MPa；含水率

为10%；生产气油比20m3/t；原油相对密度0.85；天然气相对密度0.7；50℃脱气原油粘度30mPa.s

处理参数：假设地层温度梯度为每100米3摄氏度，地表温度为20摄氏度，结合地层的厚

度就可知地层的温度分布。利用饱和压力将井筒分为不同的两段，根据不同气液流型计算相应状态下的流体参数。

计算方法：参考助教复印的压力梯度计算方法。

采用两项垂直管流的Orkiszewski方法：用压力增量迭代法计算

1） 确定起始点压力p1，计算深度增量ΔZ和分段数N： 2） 初始计算段的压力降为Δp设，并计算下端压力：

3)计算该段的平均压力

和平均温度

：

3

4)计算

和

下的流体性质参数及流动参数：

（1）原油的API度

（2）溶解气油比：因为γ

API

＞30

（3）原油的体积系数

（4）原油密度

（5）油水混合液体的密度

（6）液体的粘度 死油的粘度

活油的粘度

水的粘度：

液体的粘度：

（7）油、天然气的表面张力：

（8）水、天然气的表面张力：

（9）油、水混合物和天然气的表面张力：

（10）天然气的粘度：

其中：

（11）天然气的压缩因子

其中

条件 （12）天然气密度

5）判断流态 计算流量： 气体体积流量：

气体质量流量：

液体体积流量：

液体体积流量：

无因次气体流速：

计算公式

5）计算平均密度与摩擦损失梯度

代入可求出平均密度

根据相对粗糙度和雷诺数查表的摩擦阻力系数f=0.04，可得摩擦梯度：

6）计算压降

7）计算误差

8）利用编制的程序计算井筒的压力分布并绘图。

三、 计算不同油管直径下井筒压力梯度，选出最优解（对应现有油管尺寸）。 四、 Matlab编程产生的计算结果。

油管半径 0.025 0.026 0.027 0.028 0.029 0.03 0.031 0.032 0.033 0.034 0.035 0.036 0.037 0.038 0.039 0.04 0.041 0.042 表1 不同油管半径下的压力梯度 压力梯度 油管半径 0.0088662 0.0083937 0.0080762 0.0079108 0.0078051 0.0077286 0.00767 0.0076229 0.0075835 0.0075503 0.0075204 0.0074939 0.0074688 0.0074442 0.0074207 0.0073964 0.0073712 0.0073456

0.043 0.044 0.045 0.046 0.047 0.048 0.049 0.05 0.051 0.052 0.053 0.054 0.055 0.056 0.057 0.058 0.059 0.06 压力梯度 0.007337 0.0073486 0.0073608 0.007373 0.0073855 0.0073983 0.007411 0.0074238 0.0074366 0.0074494 0.0074623 0.007475 0.0074876 0.0075001 0.0075124 0.0075247 0.0075367 0.0075487

图1 不同油管半径对应的压力梯度

图中的最低点对应最小的压力梯度，对应的油管尺寸:31（即最优油管半径）。

2

第三题（思路）

多相流计算只考虑重力和摩擦损失，故井筒压力损失由重力损失和摩擦损失组成。有题目1可知摩擦损失，进而可得摩擦梯度?

2

?dp? ?。dL??f1

以井口为求解节点，关系式为： P2=P1+P12f+P12g

式中，P12f为摩擦损失，P12g为重力损失。 P2=Pwf；

dp?P12f=????dL?f?h；

P12g??mgh?[?o(1?fw)??wfw]gh

根据题目所给的流量与流压数据，计算出P1值

，填入下表。 液量Q(m3/d) 11.0 流压14.1 19.8 12.9 25.4 12.49 30.3 12.1 Pwf(MPa) P1(MPa)

P1（MPa）

A 1MPa

Q(m3/d) B

把结果画到图2中得到曲线A（曲线A不一定是直线），与曲线B的交点所对应的产量即为该井对应的最大产量。若产量可观，则认为投产初期油藏能量充足，为充分利用地层能量，选择自喷生产，若产量不可观，选择气举生产。

本文来源：https://www.bwwdw.com/article/vwov.html