《《采油工程》考试改革课堂考核环节试题第一次采油大作业答案.docx》由会员分享，可在线阅读，更多相关《《采油工程》考试改革课堂考核环节试题第一次采油大作业答案.docx（42页珍藏版）》请在三一文库上搜索。
1、第一题生产初期假设该井可以自喷生产，井筒中的流动可以分为两段。 下部分泡点压力以下为纯液流，上端低于泡点压力之后为气液两相流。忽略加速度压力梯度部分。600.92m,气液混合物高度为了简化计算，大概确定摩阻压力梯度的比例，讲井筒管流分为两部分， 纯液流和气液 两相流。以第一组数据为例，根据混合液的密度可以得到液柱高899.08m。(1 ) 纯液柱段摩阻压降和总压降计算: 油藏条件下的原油密度：根据油层物理第一章的内容，我们可以得到油藏条件下的溶解气油比33Rs = 25.94m /m于是可得油藏条件下的原油密度:3:oi -831.88kg/m因此该段原油密度可近似取原始条件下原油析出气体前可
2、忽略压力所引起的密度变化, 的原油密度。于是该段的平均密度：3：m = %i(1-10%)： w 10% = 848.7kg/m重力压力梯度:dPh / dh :g原油流速:qm - (qoBoqwBw)/86400雷诺数:NRe其中粘度为油水的体积加权平均值，原油的粘度根据油层物理学中相关公式得到。 根据雷诺数的大小，所给四组生产条件下的流动皆为水力光滑区。故有：,0.3164f = 14Ne摩擦损失梯度:v2dpf/dh = fvmD 2(2)气液共存段摩阻压降及总压降计算：为简化运算，气液共存段不分段，使用Orkiszewski方法进行计算。该段平均压力为：(9+0.1) /2=4.50
3、5MPa ,按照温度梯度计算中点温度值作为该段平均温度。气，液的就地流量:qgP0TZ(Rp-R)q。86400PT。qo Bo qw Bwql -86400按照Orkiszewski方法流型划分原则，分别计算Lb （泡流界限），Ls （段塞流界限）等等，并分析不同产量下的流型。四种产量下的流型都为泡流，其摩擦损失梯度可以按照液相 进行计算。但是液相流速需要计算持气率Hg,泡流持气率计算参见 采油工程原理与设计液相真实流速：VlhqlAp(1-Hg)按照相同的方法可得到摩擦阻力系数，从而计算该段的摩擦阻力损失。根据持气率亦可计算该段平均密度，从而得到气液共存段的重力压力梯度。（3）根据各段的压
4、力梯度和高度，可分别计算总的重力压差和摩阻压差，从而得到摩阻压 力梯度所占的比例。结果见附表产液摩阻压降/Pa重力压降/Pa摩阻压降占比/%11299.478651112635962.980.002419.81265.84591712444729.22;0.010225.42378.9115612331780.590.0193131.34039.60687912218334.340.0331比占阻摩0.03500.03000.02500.02000.01500.01000.00500.0000第二题一、计算地层压力：P=1.15x1000x9.8x1500=16.905MPa；选择合适的 产量
5、和压力（可自己选择）:流压（12.9MPa）一流量（19.8m3/d）。二、根据题目相应地层和流体参数计算某一油管直径的井筒压力梯度损失。（可利用matlab编程来实现）已知参数：油层厚度5m；中深1500m；地层压力系数1.15;原油饱和压力9.0MPa；含水率 为10%;生产气油比20m3/t ;原油相对密度 0.85;天然气相对密度 0.7; 50 c脱 气原油粘度30mPa.s处理参数：假设地层温度梯度为每100米3摄氏度，地表温度为 20摄氏度，结合地层的厚度就可知地层的温度分布。利用饱和压力将井筒分为不同的两段，根据不同气 液流型计算相应状态下的流体参数。计算方法：参考助教复印的压
6、力梯度计算方法。采用两项垂直管流的 Orkiszewski方法：用压力增量迭代法计算1）确定起始点压力P1,计算深度增量 A Z和分段数N：2）初始计算段的压力降为 Ap设，并计算下端压力:3）计算该段的平均压力 P和平均温度T :-pa+PiP一丁4）计算P和T下的流体性质参数及流动参数:（1）原油的API度141.5-131.5（2）溶解气油比：因为 r API30%(1404x10、产了56.06（3）原油的体积系数F = 5,615 P- + 2.257+40V.殳= 0.972+0000 叶(4)原油密度1000(+1,206x10)130(5)油水混合液体的密度(6)液体的粘度死油
7、的粘度7 = 3.0324-0.02023；幽+画活油的粘度月二 10715(5 6154+1001少刁二5.44(5.61物+150逆月(1000/J/_ =1000水的粘度:1D03-1.475X10 e(33 +1 fl r)L98i* 10 (J2+I.e T f1000液体的粘度:内叫(1-J(7)油、天然气的表面张力:424-0047(1,87 + 32)-0267 心城1。一%=iooo(8)水、天然气的表面张力：%3 丁1000_52J-8201BxlQ-?x加峋二7000_ 248-1.87_汀呻= 7罚 一 %37 凭)1 + %37 网(9)油、水混合物和天然气的表面张力
8、：巧。-八)+ %八(10)天然气的粘度：丹WIO“明其中：(9.4+0.02)(492+1.87)1J701 + 19AT+18TJ= 24-0211 = 3.5+986/(492+1,87)+0.01= 28.96(ii)天然气的压缩因子z=i+(4q+4 片总+(4+4口况+4)用其中的= 0.27%/阿)7； =(273,15+7)/7；p产pE条件计算公式。0.7普=92.22+1762=10&x(4.88-u.39/J% = 92.22+1”肛.以=W&x0,78-O.25J(12)天然气密度p =3.4844 xl03x_22_Z(T+273.15)5)判断流态计算流量：气体体积
9、流量：_4(0-&)(+273.15)293.15x86400x10气体质量流量：86400%液体体积流量：】一86400匕”以十一加液体体积流量：1000分+变也正跖二86400无因次气体流速：M包严Ls =1.701-0.7277v；/D25rL = 50 + 36% &Lm = 75+84”且严5）计算平均密度与摩擦损失梯度0.002361g (IO,a+DD】的+ 0.455+0.5691gD-XX=(l 巩+ 0,510O.OOlflg 口叫+ D+ 0.722+0 631gD代入可求出平均密度根据相对粗糙度和雷诺数查表的摩擦阻力系数f=0.04 ,可得摩擦梯度:fp佃+匕42D6)
10、计算压降as1 %诉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.0
11、490.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.007345
12、60.060.0075487图1不同油管半径对应的压力梯度图中的最低点对应最小的压力梯度，对应的油管尺寸：3，（即最优油管半径）第三题（思路）多相流计算只考虑重力和摩擦损失，故井筒压力损失由重力损失和摩擦损失组成。有题目1可知摩擦损失，进而可得摩擦梯度2以井口为求解节点，关系式为:P2=P l + P I2f+P 12g式中，Pl2f为摩擦损失，Pl2g为重力损失。P2=Pwf；X h ；Pi2g = ：mgh 十o（1 - fw） Ugh根据题目所给的流量与流压数据，计算出 Pi值,填入下表。液量 Q(m3/d)11.019.825.430.3流压PWf(MPa)14.112.912.491
13、2.1Pi(MPa)Pl把结果画到图2中得到曲线A （曲线A不一定是直线），与曲线B的交点所对应的产量即为该井对应的最大产量。 若产量可观，则认 为投产初期油藏能量充足，为充分利用地层能量，选择自喷生产，若 产量不可观，选择气举生产。第四题(思路)dp1、gs gsin 0 + v 成式中pl为多相混合物的密度；vN为多相混合物的流速；蚀为多相 混合物流动时的摩擦阻力系数;d为管径;p为压力;h为深度;g为重力加 速度；的井斜角的余角。fM=0.1, pN=874.14kg/m3, d=0.065m, h=1500m,(1)利用Beggs-Brill方法计算井筒内的压力梯度。设定一系列 的产量
14、，求出对应的井底流压，在每一个压力下运用压力增量迭代法 从井底向上计算，为简化计算过程，将整个井筒作为一段进行计算， 根据平均压力和平均温度计算 Rs, Bo,乙po, ug原油密度，天然 气密度，气液的就地流量、表观流速、质量流量等流体的性质参数， 根据Beggs-Brill方法的方法判断流型，再计算阻力系数，最后利用 公式计算压力梯度，若在误差范围以内，则取估计值为油压值，之后 绘制pt-Q的关系曲线B,根据公式2-5绘制每一个油嘴直径下产量与 油压的关系曲线，即油嘴特性曲线 G,曲线B、曲线G与IPR曲线构 成了综合协调曲线。(2)在曲线B上找到油压为2MPa的点，所对应的产量即为要 求
15、的最大产量(3)在协调曲线上找到产量为28t/d,确定与之对应或者相近的 油嘴直径，曲线G与B的交点确定实际的产量和油压， G与IPR曲 线的交点确定井底流压。(4)假设油压低于0.5MPa时该井停喷，向下计算该油压下的井底流压，根据初始值与终值的压差以及压力递减速度，即求得停喷期。第五题根据题目中流压与产量数据画出IPR曲线，并确定30t/d下的井底流压。油管直径：88.9mm;原油的压缩系数 Bo为:1.086;油水混合物的密度为：818kg/m3；换算单位30t/d=39.8m 3/d 对应井底流压 10.5MPa。(如图所示)图1 IPR曲线以Pwf (10.5MPa)为起点，计算压力
16、分布曲线。平衡点(710, 4.31),注气点(640, 3.78)即注气深度为640m,工作阀所在位置的油管压力P为3.78MPa。图确定平衡点和注气点不同注气量下的压力分布曲线表1不同气液比下的井口压力总气液比对应开口压力10000.96761220001.48351930001.70438140001.83146250001.916133图不同气液比下的井口压力图定注气压力及定井口压力下确定注气点深度及气液比最后结果：在1.5MPa的定井口压力下， 总气液比（TGLR）为2020m3/t,即注入气液比 为2000m3/t o定产量为 30t/d,注气量为60000m3/d。1000150
17、02000600050004000300020001000第六题第一部分：确定注气深度和产量已知注气压力和注气量，尽可能地获得最大可能的产量，确定注气点深度及产量的步骤：(1)假定一组产量(20t/d, 25 t/d, 30 t/d, 35 t/d),根据可提供的注气量(5万方/天)和地层生产气液比计算出每个产量的总气液比 TGLR (2378m3/ m3 ,1905m3/ m3 ,1591m3/ m3, 1366.2m3/ m3) c(2)以给定的地面注入压力 4MPa,计算环形空间气柱压力分布曲线B,用注入压力减去(0.5-0.7MPa)彳故B的平行线，即为注入气点深度线Co(3)以定井口
18、压力为起点，利用多想垂直管流，计算每个产量下的油管压力分布曲线D1, D2, D3, D4, D5。(步骤见第三部分)，他们与曲线C的交点即为各个产量所对应的的注气点以及注气深度。结果如下图所示：(4)从每个产量对应的注气点压力和深度开始。利用多相管流根据油层生产气液比向下计算每个产量对应的注气点以下的压力分 布曲线Ai, A2, A3, A4, A5.(方法同第三部分)。(5)在前面得到的IPR曲线上，根据上步计算出的产量与计算流压 的关系，交点所对应的的产量和压力即为该井在给定注气量和 井口油管压力下的最大产量 Q。综上，最大产量为34.8t/d,对应注气点地层深度为820m(6)估算最大
19、启动压力pe = Lg = 1500*850*9.8 12.5MPa第二部分：确定气举阀(1)第一个阀的下入深度：L1 = pmax 105 - 20 = 65.96m1:g(2)第二个阀的下入深度确定产量为35t/d时，利用油管内压力分布曲线：L2=L1+ Pa2-Pt1105 _10=ll2.93m7g(3)第三个阀的下入深度L3 =L2pa3 - pt2105 -10 = 158.63m；g(4)第四个阀的下入深度L4 = L3pa4 - pt3 105 -10 = 203,7m:g(5)第五个阀的下入深度L5 扎- pa5 -pt4 105 10 = 247.89m ;?g(6)第六个
20、阀的下入深度L6 = L5pa6 -pt5 105 10 = 291,7m:g(7)第七个阀的下入深度L7 =L6pa7 - pt6105 -10 335.76m:g(8)第八个阀的下入深度L8 =L7pa8 - pt7 105 -10 -377.67m：g(9)第九个阀的下入深度L9 = L8pa9 - pt8 105 -10 = 417.04m:g(10)第十个阀的下入深度Lii =Lio - pa11 pt10105 -10 = 490.1m:g(12)第十二个阀的下入深度L12 ML11 pa12 -pt11105 10 =524.41m用(13)第十三个阀的下入深度L13 = L12
21、 pa13 - pt12105 -10 = 556.19m二 g(14)第十四个阀的下入深度L14 =L13 pa14 -pt13105 -10 587.98m：g(15)第十五个阀的下入深度L15 =L14pa15 - pt14105 -10 = 619.76m:g(16)第十六个阀的下入深度L16二L15pa16 - pt9105 - 10 = 646.49m:g(17)第十七个阀的下入深度pa17 - p t165L17 =L16 10 70 =673.2m:g(18)第十八个阀的下入深度L18 =L17pa18 - pt17105 -10 = 697.40m:g(19)第十九个阀的下入
22、深度pa19 - p t185L19 = L18 : 10 70 = 721.578m:g(20)第二十个阀的下入深度L20=Li9 pa20 - pt19105 10 = 745.25m：g(21)第二十一个阀的下入深度L21 =L20pa21 -pt20 105 -10 770.084m7g(22)第二十二个阀的下入深度L 22 = L 21 pa22 - pt21105 -10 = 793.56m7 g(23)第二十三个阀的下入深度L23=L22 - pa23 -pt22105 -10 =823.56m；g823.56m注气点深度，共需气举阀23个。第三部分：油管压力分布的计算采用两项垂
23、直管流的Orkiszewski方法：用压力增量迭代法计算取产量q=15t/d1) 确定起始点压力Pi,计算深度增量AZ和分段数N:p1=1.5MPa, AZ=100rm N=152) 初始计算段的压力降为 Ap设，并计算下端压力：p 设=0.5MPap2设=P1 + 设=1.5 + 0.5 = 2MPaP2 设 +P1 p=F-3)计算该段的平均压力p和平均温度T：1.5 2 =1.75MPa2T 十 20C4)计算p和T下的流体性质参数及流动参数:(1)原油的API度141.5api 131.5=34.97o(2)溶解气油比：因为丫 API 30g (1.404 10p)1.1871 黑3
24、.3/ 3、Rs =父101囱书2书59.67 = 40.0341 (m /m )56.06(3)原油的体积系数VF =5.615RSJ十 2.25T+ 40= 294.0539；oBO = 0.972 0.000147 F 1.175 = 1.0889(4)原油密度1000( o 1.206 10与Rs g)130=811.661(kg/m3)(5)油水混合液体的密度一、。(1一 fw)wfw = 830.4949(6)液体的粘度死油的粘度z= 3.0324- 0.02023 API = 2.325y = 10z = 211.3489x 二 y(32 1.8T11.163 = 1.46071
25、0x -11000=0.0279活油的粘度-0.515A= 10.715(5.615Rs + 100)=0.04510.338B = 5.44 5.615Rs 150=0.7339oA(1000od)b1000=0.0063水的粘度:we1.003-1.479 10- 32 1.8T 1.982 10- 32 1.8T1000= 1.0411液体的粘度:l=o1-fwwfw=0.1098(7)油、天然气的表面张力:二 og42.4 0.047 (1.8T +32 ) 0.267，pi 父 e.015*/P =0 02501000(8)水、天然气的表面张力:3.62575 10 - p二 23.
26、33 ,。00 巾.0715137.7852.5一8.7018 或 P=0.05101000一 248 -1.8T_wg = c” ,- 23.33 - 137.78206二 137.78 =0.0717(9)油、水混合物和天然气的表面张力:二 1 =1-fw二 wgfw=0.0296(10)天然气的粘度:出 = K 父10(一俏=0.0109 g其中:1.5= 106.84979.4 0.02M 492 1.8TK 二701 19M 1.8Ty =2.4 -0.2x =1.2888x =3.5 986/ 492 1.8T0.01M =5.5559M =28.96 g =20.272(11)天
27、然气的压缩因子Z =1+(A + A2/Tr +A3/T；)Pr +(a +A/Tr )pR +A6pR/Tr3其中：R =0.27pr/ ZTrTr = 273.15 T /TcPr = P / PcAi=0.31506, A2=-1.0467, Ai=-0.5783,Ai=0.5353, Ai=-0.6123, Ai=0.6815。条件计算公式飞g 27Tc = 92.22 +176.67”, Pc = 10隈(4.88 - 0.39尸g )飞g 0.7Tc = 92.22 +176.67%, pc = 106 秋 4.78 0.25yg )综上，利用试算法可得 Z=0.9502(12)天
28、然气密度y：g =3.4844 10“ g=14.772gZ T 273.155）判断流态计算流量：气体体积流量:QgZq(Rp - Rs)(T 273.15)3 p二 0.025 m3 /s)293.15 86400 0p 10气体质量流量:1.206qo(Rp - Rs) gW =2o-L_p= 0.4815 kg/sg86400 o液体体积流量:Qi =q1Bo(1 一 fw) Bwfw= 3.6 10-4 m3/s86400 o(1 - fw)wfw液体体积流量:1.206qRs g10005 W1 二o: 0.3007 kg /s86400一 qg / P1、1/4 无因次气体流速：
29、9=7（） =57.485Ap gcLb = 1.701 0.7277vt2 /D (由于 LB因此取Pi +即卜作为新的Pi从第二步重新开始计算，知道满足需求后再 以上端的末端压力作为下段的起点压力开始计算下一段。8）利用编制的程序计算井筒的压力分布并绘图。第七题(1)已知两组S10层产能测试表。由地层梯度 0.5100.571psi/ft算得S10层位的测试点的地层压力应该在 21.57MPa24.15MPa之间。对比产能测试结果， 发现第二组数据符合要求，认为第二组数据准确。由测试点的流压大于饱和压力，产液含水率为零。得 S10 层采油指数：J = qo(test)一=226=77.93
30、m3/(d MPa) Pr - Pwf(test) 22.45-19.55当 Pwf (test)之 Pb 时，q。= J( Pr - Pwf ),C C PwfC C Pwf1Pwf(test) ( Pb 时，q。=qb+qc 1 0.2 0.8 一Pb、Pb 川13其中 qb =J(Pr - Pb) =77.93 (22.45-2.31) = 1569.51m /d =qc -JPb/1.8 =77.93 2.31 = 1.8 =100.01m3/dPwf (MPa)00.511.522.3122.45Qo (m3/d)1669.521661.441645.871622.801592.23
31、1569.510.00得S10层流入动态曲线如图:Qo (m3/d)(2)由于各小层属于同一油藏，认为油层的物性相似。由 j =空koha代入 S10 层数据得 Bo(ln 三+s)=164.87mPa s出(4-3 s)rw 4rw 4所以S2层的采油指数为2二 koha_re3JoBo(ln,一 s)rw 42二 4.578 3.7 86.4164.873 ,一、=55.75m /(d MPa)各小层温度相差不大，认为饱和压力相似。A pwf (test)之 Pb 时，qo=J(PrPwf)Pwf (test)父 Pb时，qo =qb +q。1-0.2pf -0.8Pb其中 qb =J(p
32、r - pb) =55.75 (15.72-2.31) = 747.61m3/d =qc =Jpb/1.8 =55.75 2.31，1.8 = 71.55m3/ dPwf(MPa)00.511.522.3115.72Qo(m3/d)819.15813.38802.24785.73763.86747.610.00S4(1)层的采油指数为,2二 kohaJ 二口 oBo(ln - s) rw42 二 1.015 1.31 86.4164.87=4.38m3/(d MPa)3其中 qb = J(pr - Pb) =63.16m /d3qc =Jpb/1.8 -5.621m /dPwf (MPa)00
33、.511.522.3116.73Qo (m3/d)68.7868.3367.4566.1564.4463.160.00Qo Cm3/d)S4(2)层的采油指数为2二 kohaJ =7 3oB(ln - - s)rw42 二 5.21 1.31 86.4164.873_=44.36m /(d MPa)其中 qb = J(pr - Pb) =639.67m3/d qc =Jpb/1.8 = 56.93m3/dPwf(MPa)00.511.522.3116.73Qo(m3/d)696.60691.99683.13670.00652.60639.670S4(3)层的采油指数为2二 kohaJ 二 J
34、3jBo(ln_es)rw42 二 1.25 2.314 86.4164.873=9.52m /(d MPa)其中 qb = J(pr - pb) = 137.28m3/dQc (m3/d)qc =Jp1.8 =12.22m3/dPwf(MPa)00.511.522.3116.73Qo(m3/d)149.49148.57146.66143.82140.07137.280S6(1)层的采油指数为.2 二haJ 二。B(ln e - s)rw42 二 7.77 1.019 86.4164.87=26.06m3/(d MPa)其中其中 qb = J(pr - pj = 371.36m3 /d_ _3
35、.qc =Jpb/1.833.44m /dPwf(MPa)00.511.522.3116.56Qo(m3/d)404.80402.10396.89389.18378.96371.360Qq ( m3/d)S6(2)层的采油指数为2二 kohaoBo(ln - - s)rw42 二 1.68 1.143 86.4164.873=6.32m /(d MPa)其中 qb = J( pr - Pb) = 90.06m3 / d = qc = Jpb/1.8 =8.11m3 /dPwf (MPa)00.511.522.3116.56Qo (m3/d)98.1797.5196.2594.3891.9090
36、.060Qq ( m3/d)总的流入动态曲线为:Pwf(MPa)00.511.522.3115.7216.5616.7322.45Qo(m3/d)3906.513883.323838.493772.063684.053618.65610.51468.91445.760.00Qo (rn3/d)第八题解：将K7层的产能测试数据代入达西公式：% = 2减/3二H4出3和 2+E）% 44&（1口2-?+）= 0 52细产仪得到：.将地下原油视为同一性质。由题目中“生产气油比较低，饱和压力较 低；地饱压差较大，生产过程中几乎无自由气。”得知可以将井筒中 流动视为单相流。（K7层泡点压力为2.93MP
37、a, 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)层:0123456TS产“3/d2383 06522 30.0292076 9731923. 9171770.861617. 8041464. 7401311. 6921158.635流压I*9101112_J3141515.515. ST产量1005.579352.5
38、2316997669546.4107393. 354524CL29S2ST . 2420410.T1393018K2(l)层1PR曲线0 SOO 1000 ISOO 200025003000产量m3/dj二丝二竺二M%)采油指数 ,1-K2(2)层流%0123456T6产“3/d3583.0373355.1093127. 182899. 2522671.3232443. 3942215. 466liger.5371759.609流压i*910n12IS141515.515.72产量n)3/d1531.W1303 7511075.B23E47 E943619.9657392 0372164. 1036S3.14429001000200030004000产量m%/dj =二 32=227 9斓MPa)采油指数K6(1)层采油指数K6(2)层01000200030004000产量m3/dOOi 1 1J 二曳二业巴二212川加) 采油指数 / 一K7层5000ifiUMFa01234567T9非52 43703-4的3454.9145206.352 非718042 709,249J46Q. 691?Zli_15S1郭35的MlJQEJIPii10111213141515.5161&.516.9产鼠.3/ d1715.029ME6, 41217 919W. 36铝720. 8OT44