不含水组合型水平井流入动态模型.docxVIP

不含水组合型水平井流入动态模型 目前，对水平井流量动态的研究主要集中在两个方面：一是考虑不同的表皮系数，通过水平井损失率公式计算流量动态。该方法需要大量精确的地层参数，实际应用比较困难。其次，回归生产数据的动态回归关系是真正意义上的动态流量预测。后者通常是对溶解气驱油藏水平井的生产数据或数值模拟结果进行回归,建立水平井流入动态方程,这些方程仅适用于油藏压力低于饱和压力的油藏。而对于非溶解气驱、已进入含水期的绝大多数油藏中的水平井,尚缺乏可用的流入动态模型和预测计算方法。笔者通过筛选现有的溶解气驱水平井流入动态方程,研究适用于油藏压力高于饱和压力情况的组合型IPR方程,并将其扩展到油、气、水三相的情况,建立一套适用于油、气、水三相和整个油藏压力范围的水平井IPR模型和预测计算方法。 1 刘春有不同井角的特性分析 目前共有4种溶解气驱水平井IPR方程,分别由Cheng、Bendakhlia、刘想平和黄炳光提出。Bendakhlia用两种三维三相黑油模拟器研究得到一种溶解气驱水平井的流入动态关系;黄炳光将研究溶解气驱垂直井流入动态的方法用于水平井分析,根据Bendakhlia的结果,也得到一种相似的流入动态关系方程。由于方程本身数学形式的原因,使用本文的研究方法无法利用Bendakhlia和黄炳光的方程进一步建立油、气、水三相的流入动态模型(限于篇幅,推导与证明过程从略)。因此本文中仅以其他两种方法为基础进行研究。 刘想平用黑油模拟器研究了溶解气驱油藏水平井的流入动态,采用回归方法得到如下方程: qoqomax=1-apwfpr-(1-a)2(pwfpr)2-a(1-a)(pwfpr)3.(1)qoqomax=1?apwfpr?(1?a)2(pwfpr)2?a(1?a)(pwfpr)3.(1) 式中,qo为产油量,t/d;qomax为最大产油量,t/d;pwf为井底流压,MPa;pr为油藏压力,MPa;a为待定参数,0≤a≤1。 Cheng对溶解气驱油藏中斜井和水平井进行了数值模拟,并用回归的方法得到了如下不同井斜角油井的IPR方程: qoqomax=A-Bpwfpr-C(pwfpr)2.(2)qoqomax=A?Bpwfpr?C(pwfpr)2.(2) 式中,A,B,C为井斜角的函数,对于水平井,A=0.9885,B=-0.2055,C=1.1818。 方程(1),(2)均是从溶解气驱水平井的数值模拟结果回归而来,因此仅能用于油藏压力低于饱和压力的溶解气驱油藏。 2 不含水组的水平井的ipr模型 2.1 组合型ipr曲线a点及b点的确定 图1为不含水组合型IPR曲线。在prpb的油藏中,当pwfpb时,为单相流动,采液指数为常数,IPR曲线为直线,此时的流入动态关系为 qo=J(pr-pwf),qb=J(pr-pb). (3) 式中,J为采液指数,t/(d·MPa);pb为饱和压力,MPa;qb为饱和压力pb下的产量,t/d。 当pwfpb时,为两相渗流,IPR为曲线,分别用pb,qc代替溶解气驱IPR方程(1)中的pr,qomax并进行产量叠加即可得到组合型IPR曲线上pwfpb时的关系方程: qo=qb+qc[1-apwfpb-(1-a)2(pwfpb)2-a(1-a)(pwfpb)3].(4)qo=qb+qc[1?apwfpb?(1?a)2(pwfpb)2?a(1?a)(pwfpb)3].(4) 其中 qc=qomax-qb. 由于图1中IPR曲线在A点(pwf=pb)的连续性,方程(3),(4)在A点的导数应相等。对方程(3),(4)分别求导并将pwf=pb代入使二者相等得到 -J=qc[-a1pb-2(1-a)2pwfp2b-3a(1-a)(pwfpb)21pb]|pwf=pb?J=qc[?a1pb?2(1?a)2pwfp2b?3a(1?a)(pwfpb)21pb]∣∣pwf=pb, (5) J=(2-a2)qcpb,qc=Jpb(2-a2).(6)J=(2?a2)qcpb,qc=Jpb(2?a2).(6) 根据方程(4),最后得到pwfpb时采液指数的表达式为 J=qo{(pr-pb)+pb2-a2[1-apwfpb-(1-a)2×(pwfpb)2-a(1-a)(pwfpb)3]}-1.(7)J=qo{(pr?pb)+pb2?a2[1?apwfpb?(1?a)2×(pwfpb)2?a(1?a)(pwfpb)3]}?1.(7) 2.2 产油量流压测试点 使用同样方法得到的以Cheng溶解气驱方程为基础的不含水组合型水平井IPR方程为 {qo=J(pr-pwf),pwf≥pb;qo=qb+Jpb(B+2C)[A-Bpwfpb-C(pwfpb)2],pwfpb.(8)???????qo=J(pr