采油工程考试改革课堂考核环节试题第一次采油大作业复习资料要点.docx

第一题生产初期假设该井可以自喷生产，井筒中的流动可以分为两段。下部分泡点压力以下为纯液流，上端低于泡点压力之后为气液两相流。忽略加速度压力梯度部分。为了简化计算，大概确定摩阻压力梯度的比例，讲井筒管流分为两部分，纯液流和气液两相流。以第一组数据为例，根据混合液的密度可以得到液柱高600.92m，气液混合物高度899.08m。（1） 纯液柱段摩阻压降和总压降计算：油藏条件下的原油密度：根据油层物理第一章的内容，我们可以得到油藏条件下的溶解气油比于是可得油藏条件下的原油密度：原油析出气体前可忽略压力所引起的密度变化，因此该段原油密度可近似取原始条件下的原油密度。于是该段的平均密度：重力压力梯度：原油流速：雷诺数：其中粘度为油水的体积加权平均值，原油的粘度根据油层物理学中相关公式得到。根据雷诺数的大小，所给四组生产条件下的流动皆为水力光滑区。故有：摩擦损失梯度：（2） 气液共存段摩阻压降及总压降计算：为简化运算，气液共存段不分段，使用Orkiszewski方法进行计算。该段平均压力为：（9+0.1）/2=4.505MPa，按照温度梯度计算中点温度值作为该段平均温度。气，液的就地流量：按照Orkiszewski方法流型划分原则，分别计算LB（泡流界限），LS（段塞流界限）等等，并分析不同产量下的流型。四种产量下的流型都为泡流，其摩擦损失梯度可以按照液相进行计算。但是液相流速需要计算持气率Hg，泡流持气率计算参见采油工程原理及设计。液相真实流速：按照相同的方法可得到摩擦阻力系数，从而计算该段的摩擦阻力损失。根据持气率亦可计算该段平均密度，从而得到气液共存段的重力压力梯度。（3）根据各段的压力梯度和高度，可分别计算总的重力压差和摩阻压差，从而得到摩阻压力梯度所占的比例。结果见附表产液摩阻压降/Pa重力压降/Pa摩阻压降占比/%11299.478651112635962.980.0024 19.81265.84591712444729.220.0102 25.42378.9115612331780.590.0193 31.34039.60687912218334.340.0331 第二题一、 计算地层压力：P=1.15x1000x9.8x1500=16.905MPa；选择合适的产量和压力(可自己选择):流压(12.9MPa)流量(19.8m3/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)计算该段的平均压力和平均温度：4)计算和下的流体性质参数及流动参数：（1）原油的API度（2）溶解气油比：因为API30（3）原油的体积系数（4）原油密度（5）油水混合液体的密度（6）液体的粘度死油的粘度活油的粘度水的粘度：液体的粘度：（7）油、天然气的表面张力：（8）水、天然气的表面张力：（9）油、水混合物和天然气的表面张力：（10）天然气的粘度：其中：（11）天然气的压缩因子其中条件计算公式（12）天然气密度5）判断流态计算流量：气体体积流量：气体质量流量：液体体积流量：液体体积流量：无因次气体流速：5）计算平均密度及摩擦损失梯度代入可求出平均密度根据相对粗糙度和雷诺数查表的摩擦阻力系数f=0.04，可得摩擦梯度：6）计算压降7）计算误差8）利用编制的程序计算井筒的压力分布并绘图。三、 计算不同油管直径下井筒压力梯度，选出最优解（对应现有油管尺寸）。四、 Matlab编程产生的计算结果。表1 不同油管半径下的压力梯度油管半径压力梯度油管半径压力梯度0.0250.00886620.0430.0073370.0260.00839370.0440.00734860.0270.00807620.0450.00736080.0280.00791080.0460.0073730.0290.00780510.0470.00738550.030.00772860.0480.00739830.0310.007670.0490.0074110.0320.00762290.050.00742380.0330.00758350.0510.00743660.0340.00755030.0520.00744940.0350.00752040.0530.00746230.0360.00749390.0540.0074750.0370.00746880.0550.00748760.0380.00744420.0560.00750010.0390.00742070.0570.00751240.040.00739640.0580.00752470.0410.00737120.0590.00753670.0420.00734560.060.0075487图1 不同油管半径对应的压力梯度图中的最低点对应最小的压力梯度，对应的油管尺寸:（即最优油管半径）。第三题（思路）多相流计算只考虑重力和摩擦损失，故井筒压力损失由重力损失和摩擦损失组成。有题目1可知摩擦损失，进而可得摩擦梯度。12以井口为求解节点，关系式为： P2=P1+P12f+P12g式中，P12f为摩擦损失，P12g为重力损失。P2=Pwf；P12f=；根据题目所给的流量及流压数据，计算出P1值，填入下表。液量Q(m3/d)11.0 19.825.430.3流压Pwf(MPa)14.1 12.912.4912.1P1(MPa)P1（MPa）A1MPaBQ(m3/d)把结果画到图2中得到曲线A（曲线A不一定是直线），及曲线B的交点所对应的产量即为该井对应的最大产量。若产量可观，则认为投产初期油藏能量充足，为充分利用地层能量，选择自喷生产，若产量不可观，选择气举生产。第四题（思路）1、gsin+v/d*式中为多相混合物的密度;v为多相混合物的流速;f为多相混合物流动时的摩擦阻力系数;d为管径;p为压力;h为深度;g为重力加速度; 为井斜角的余角。f=0.1，=874.14kg/m3，d=0.065m，h=1500m，（1）利用Beggs-Brill方法计算井筒内的压力梯度。设定一系列的产量，求出对应的井底流压，在每一个压力下运用压力增量迭代法从井底向上计算，为简化计算过程，将整个井筒作为一段进行计算，根据平均压力和平均温度计算Rs，Bo，Z，o,g，原油密度，天然气密度，气液的就地流量、表观流速、质量流量等流体的性质参数，根据Beggs-Brill方法的方法判断流型，再计算阻力系数，最后利用公式计算压力梯度，若在误差范围以内，则取估计值为油压值，之后绘制pt-Q的关系曲线B，根据公式2-5绘制每一个油嘴直径下产量及油压的关系曲线，即油嘴特性曲线G，曲线B、曲线及IPR曲线构成了综合协调曲线。（2）在曲线B上找到油压为2MPa的点，所对应的产量即为要求的最大产量（3）在协调曲线上找到产量为28t/d，确定及之对应或者相近的油嘴直径，曲线及的交点确定实际的产量和油压，及曲线的交点确定井底流压。（）假设油压低于0.5MPa时该井停喷，向下计算该油压下的井底流压，根据初始值及终值的压差以及压力递减速度，即求得停喷期。第五题一、 根据题目中流压及产量数据画出IPR曲线，并确定30t/d下的井底流压。油管直径:88.9mm;原油的压缩系数Bo为:1.086;油水混合物的密度为：818kg/m3 ;换算单位30t/d=39.8m3/d对应井底流压10.5MPa。（如图所示）图1 IPR曲线二、 以Pwf（10.5MPa）为起点，计算压力分布曲线。平衡点（710，4.31），注气点（640，3.78）即注气深度为640m，工作阀所在位置的油管压力P为3.78MPa。图 确定平衡点和注气点三、 不同注气量下的压力分布曲线表1 不同气液比下的井口压力总气液比对应井口压力10000.96761220001.48351930001.70438140001.83146250001.916133图 不同气液比下的井口压力图 定注气压力及定井口压力下确定注气点深度及气液比最后结果：在1.5MPa的定井口压力下，总气液比(TGLR)为2020m3/t，即注入气液比为2000m3/t。定产量为30t/d,注气量为60000m3/d。第六题第一部分：确定注气深度和产量已知注气压力和注气量，尽可能地获得最大可能的产量，确定注气点深度及产量的步骤：（1） 假定一组产量（20t/d，25 t/d，30 t/d，35 t/d），根据可提供的注气量（5万方/天）和地层生产气液比计算出每个产量的总气液比TGLR（2378m3/ m3 ,1905m3/ m3 ,1591m3/ m3，1366.2m3/ m3）。（2） 以给定的地面注入压力4MPa，计算环形空间气柱压力分布曲线B，用注入压力减去（0.5-0.7MPa）做B的平行线，即为注入气点深度线C。（3） 以定井口压力为起点，利用多想垂直管流，计算每个产量下的油管压力分布曲线D1，D2，D3，D4，D5。（步骤见第三部分），他们及曲线C的交点即为各个产量所对应的的注气点以及注气深度。结果如下图所示：（4） 从每个产量对应的注气点压力和深度开始。利用多相管流根据油层生产气液比向下计算每个产量对应的注气点以下的压力分布曲线A1，A2，A3，A4，A5.（方法同第三部分）。（5） 在前面得到的IPR曲线上，根据上步计算出的产量及计算流压的关系，交点所对应的的产量和压力即为该井在给定注气量和井口油管压力下的最大产量Q。综上，最大产量为34.8t/d，对应注气点地层深度为820m。（6） 估算最大启动压力第二部分：确定气举阀（1） 第一个阀的下入深度：（2）第二个阀的下入深度确定产量为35t/d时，利用油管内压力分布曲线：=112.93m（3）第三个阀的下入深度158.63m（4）第四个阀的下入深度（5）第五个阀的下入深度（6）第六个阀的下入深度（7）第七个阀的下入深度（8）第八个阀的下入深度（9）第九个阀的下入深度（10）第十个阀的下入深度（11）第十一个阀的下入深度（12）第十二个阀的下入深度（13）第十三个阀的下入深度（14）第十四个阀的下入深度（15）第十五个阀的下入深度（16）第十六个阀的下入深度（17）第十七个阀的下入深度（18）第十八个阀的下入深度（19）第十九个阀的下入深度（20）第二十个阀的下入深度（21）第二十一个阀的下入深度（22）第二十二个阀的下入深度（23）第二十三个阀的下入深度823.56m注气点深度，共需气举阀23个。第三部分：油管压力分布的计算采用两项垂直管流的Orkiszewski方法：用压力增量迭代法计算取产量q=15t/d1） 确定起始点压力p1，计算深度增量Z和分段数N：p1=1.5MPa，Z=100m，N=152） 初始计算段的压力降为p设，并计算下端压力：3)计算该段的平均压力和平均温度：4)计算和下的流体性质参数及流动参数：（1）原油的API度（2）溶解气油比：因为API30（m3/m3）（3）原油的体积系数（4）原油密度（5）油水混合液体的密度（6）液体的粘度死油的粘度活油的粘度=0.0451=0.7339=0.0063水的粘度：=1.0411液体的粘度：=0.1098（7）油、天然气的表面张力：=0.0250（8）水、天然气的表面张力：=0.0715=0.0510=0.0717（9）油、水混合物和天然气的表面张力：=0.0296（10）天然气的粘度：=0.0109其中：=106.8497=1.2888=5.5559=20.272（11）天然气的压缩因子其中A1=0.31506，A2=-1.0467，A1=-0.5783，A1=0.5353，A1=-0.6123，A1=0.6815。条件计算公式综上，利用试算法可得Z=0.9502（12）天然气密度=14.7725）判断流态计算流量：气体体积流量：气体质量流量：液体体积流量：液体体积流量：无因次气体流速：=57.485（由于LB0.13，取0.13）=79.254经判断，该段属于断塞流5）计算平均密度及摩擦损失梯度根据段塞流相关公式，可求出相关参数vs=0.及：代入可求出平均密度=98.55kg/m3根据相对粗糙度和雷诺数查表的摩擦阻力系数f=0.04，可得摩擦梯度：=2019.286）计算压降：=0.7746Mpa7）计算误差，取因此取作为新的从第二步重新开始计算，知道满足需求后再以上端的末端压力作为下段的起点压力开始计算下一段。8）利用编制的程序计算井筒的压力分布并绘图。第七题(1) 已知两组S10层产能测试表。由地层梯度0.5100.571psi/ft算得S10层位的测试点的地层压力应该在21.57MPa24.15MPa之间。对比产能测试结果，发现第二组数据符合要求，认为第二组数据准确。由测试点的流压大于饱和压力，产液含水率为零。得S10层采油指数：当时，时，其中=Pwf（MPa）00.511.522.3122.45Qo（m3/d）1669.52 1661.44 1645.87 1622.80 1592.23 1569.51 0.00 得S10层流入动态曲线如图： (2) 由于各小层属于同一油藏，认为油层的物性相似。由代入S10层数据得所以S2层的采油指数为各小层温度相差不大，认为饱和压力相似。当时，时，其中=Pwf（MPa）00.511.522.3115.72Qo（m3/d）819.15 813.38 802.24 785.73 763.86 747.61 0.00 S4(1)层的采油指数为其中Pwf（MPa）00.511.522.3116.73Qo（m3/d）68.78 68.33 67.45 66.15 64.44 63.16 0.00 S4(2)层的采油指数为其中Pwf（MPa）00.511.522.3116.73Qo（m3/d）696.60 691.99 683.13 670.00 652.60 639.67 0S4(3)层的采油指数为其中Pwf（MPa）00.511.522.3116.73Qo（m3/d）149.49 148.57 146.66 143.82 140.07 137.28 0S6(1)层的采油指数为其中其中Pwf（MPa）00.511.522.3116.56Qo（m3/d）404.80 402.10 396.89 389.18 378.96 371.36 0S6(2)层的采油指数为其中=Pwf（MPa）00.511.522.3116.56Qo（m3/d）98.17 97.51 96.25 94.38 91.90 90.06 0总的流入动态曲线为：Pwf（MPa）00.511.522.3115.7216.5616.7322.45Qo（m3/d）3906.51 3883.32 3838.49 3772.06 3684.05 3618.65 610.51 468.91 445.76 0.00 第八题解：将K7层的产能测试数据代入达西公式：得到：将地下原油视为同一性质。由题目中“生产气油比较低，饱和压力较低；地饱压差较大，生产过程中几乎无自由气。”得知可以将井筒中流动视为单相流。（K7层泡点压力为2.93MPa，K8层为1.8MPa）小层数据如下层位厚度渗透率压力K2(1)4.882.61215.57K2(2)5.513.44515.72K6(1)2.233.10616.64K6(2)6.072.9216.56K79.322.22116.9K872.61517.1将上表数据代入达西公式可以得到产量及井底流压之间的关系：K2(1)层：采油指数K2(2)层采油指数K6(1)层采油指数K6(2)层采油指数K7层采油指数K8层采油指数合采数据：第 16 页