延安大学：《采油工程》课程教学资源（电子教案）第一章 油井流入动态与井筒多相流动计算

第一章油井流入动态 目的要求 讲授流体流入动态的基本内容。主要包括垂直井流入动态和水平井流入动态两个方面。 要求学生了解单相及多相油气流入动态曲线的特征及规律。 课时:4学时 重点授课内容提要 、垂直井流入动态 (一)单相液体的流入动态 达西滲流 圆形地层 定压边界(稳定流):%0.543k从(D-) HoBo(In e 1 封闭边界(拟稳定流):q。= 0.543k从(p-Dr) bo(ln 3 非圆形地层 0.543kh(p-p) HB(nX-0.75+s) 0.543kh(P,-pn) q B(nX-0.75+s) 其中: 0.543kh B(nx-0.75+s 采油指数:是一个反应油层性质、厚度、流体参数、泄油面积、完井条件等的综合指 标 数值上等于单位生产压差下的油井产量活产油量与生产压差之比:J。 IPR曲线斜率的负倒数

1 第一章 油井流入动态 目的要求 讲授流体流入动态的基本内容。主要包括垂直井流入动态和水平井流入动态两个方面。 要求学生了解单相及多相油气流入动态曲线的特征及规律。 课时：4 学时 重点授课内容提要  一、垂直井流入动态 （一）单相液体的流入动态 1 达西渗流 圆形地层： 定压边界（稳定流）： 封闭边界（拟稳定流）： 非圆形地层： 其中： 采油指数：是一个反应油层性质、厚度、流体参数、泄油面积、完井条件等的综合指 标。 数值上等于单位生产压差下的油井产量活产油量与生产压差之比： IPR 曲线斜率的负倒数 ) 2 1 (ln 0.543 ( ) s r r B kh p p q w e o o r wf o      ) 4 3 (ln 0.543 ( ) s r r B kh p p q w e o o r wf o      B  X s k h p p q o o r w f o     ln 0.75 0.543 ( )  B  X s k h p p q o o r w f o     ln 0.75 0.543 ( )  qo  J o(p r  pwf ) B  X s k h J o o o    ln 0.75 0.543  r wf o o p p q J   qo  J o(p r  pwf ) wf o o dp dq J  

2、非达西渗流 条件:油井产量很高时,在井底附近不再符合线性滲流,呈现高速非线性渗流。 Pr-p 线的斜率 Pr-p =C+D q。 (二)溶解气驱油藏油井的流入动态 条件:地层压力小于饱和压力 0.543kh 0.75+ u, B Vogel方程 1968年,沃格尔对不同流体性质、气油比、相对渗透率、井距、压裂井、污染井等各种 情况下的21个溶解气驱油藏进行了计算 无因次IPR曲线 都有类似的形状 只是高粘度油藏 及油层损害严重 时差别较大 0.2040.60.8 0102030405060 产量(M3ld) ogel方程: Vogel型方程: 1-0. 0. 2非完善井 Vogel方程的修正 油井的非完善性 ●打开性质不完善,如射孔完成

2 2、非达西渗流 条件：油井产量很高时，在井底附近不再符合线性渗流，呈现高速非线性渗流。 （二）溶解气驱油藏油井的流入动态 条件：地层压力小于饱和压力 1 Vogel 方程 1968 年，沃格尔对不同流体性质、气油比、相对渗透率、井距、压裂井、污染井等各种 情况下的 21 个溶解气驱油藏进行了计算。 Vogel 方程： Vogel 型方程： 2 非完善井 Vogel 方程的修正 油井的非完善性：  打开性质不完善，如射孔完成 2 pr  pwf  Cqo  Dqo o o r wf C Dq q p p                r wf p p o o r o w e o dp B k s r r k h q  ln 0.75 0.543 无因次 IPR 曲线 都有类似的形状， 只是高粘度油藏 及油层损害严重 时差别较大。 2 max 1 0.2 0.8                    r w f r w f o o p p p p q q 2 max 1 (1 )                     r wf r wf o o p p a p p a q q

●打开程度不完善,如未全部钻穿油层 打开程度和打开性质双重不完善 酸化、压裂等措施 油层受到损害 流动效率:在相同产量下的理想生产压差与实际生产压差之比 Prf.i p,- p in--0.75 S=0,FE=1,完善井 FE s>0,FE<1,不完善井 n-0.75+s s<0,FE>1,超完善井 Standing 方法 曲线 (三)组合型流入动态曲线 条件:P>Pb P<P中 靠近井筒部分为油气两相流动,远离井筒部分为单相流动

3  打开程度不完善，如未全部钻穿油层  打开程度和打开性质双重不完善  酸化、压裂等措施  油层受到损害 流动效率：在相同产量下的理想生产压差与实际生产压差之比。 Standing 方法： （三）组合型流入动态曲线 条件： 靠近井筒部分为油气两相流动，远离井筒部分为单相流动。 r wf a r wf i p p p p FE , ,    pwf ,i  p r  (p r  pwf ,a)FE s r r r r FE w e w e     ln 0.75 ln 0.75 0， 1，超完善井 0， 1，不完善井 0， 1，完善井       s FE s FE s FE 2 , , max( 1) 1 0.2 0.8                     r w f i r w f i o F E o p p p p q q pwf ,i  p r  (p r  pwf ,a)FE IPR 曲线 pr  pb Pwf  Pb

P P 时 16 q。=J0(p P8 q=4+q1-02/ P gb=J(pr-P) =J(-B)→ qb +q|1-0. 4人-8)37 Pb 吻。=-02 pb p J。Pb 4=J(P,-p2)+2P 1-022|-08 1.8 Pu Pb (四)综合IPR曲线 Petrobras提出 油气两相的PPR曲线 ◆实质是按含水率 油气水三相的IPR曲线 纯水的IPR曲线 取纯油IPR曲线 和水IPR曲线的 加权平均值 ◆加权的物理量可 产液量

4 时 （四）综合 IPR 曲线 pwf  Pb                             2 1 0.2 0.8 b wf b wf o b c p p p p q q q pwf  Pb qo  J o(p r  pwf ) 时 时 ( ) b o pr pb q  J  qo  J o(p r  pwf ) o wf o J dp dq                               2 1 0.2 0.8 b wf b wf o b c p p p p q q q 2 0.2 1.6 o c wf c wf b b dq q p q dp p p    pwf  pb 1.8 o b c J p q                               2 1 0.2 0.8 1.8 ( ) b wf b o b wf o o r b p p p J p p q J p p  Petrobras 提出  实质是按含水率 取纯油 IPR 曲线 和水 IPR 曲线的 加权平均值  加权的物理量可 以是产量，也可以 是流压

1、按产量加权平均→q1=f,qn+(1-f)。 Pu> Pb q,=J(P,-Pu) 4=(1-fJ(p,-P2)+ 1-0. Put 2、按流压加权平均→p=fmP+(1-fm)P 油气两相的PR曲线 油气水三相的IPR曲线 纯水的IPR曲线 产液量/m3d (五)多层油藏油井流入动态 1、层间无干扰 全井|PR=各层PR的叠加

5 1、按产量加权平均 2、按流压加权平均 （五）多层油藏油井流入动态 1、层间无干扰 全井 IPR=各层 IPR 的叠加 t w w w qo q  f q  (1 f ) pwf  pb ( ) t L pr pwf q  J  pwf  pb 1 0.2 0.8 ( ) 1.8 (1 ) ( ) 2 w L r wf b wf b L b wf t w L r b f J p p p p p J p p q f J p p                                         wf w wf w w pwf o p f p f , ,   (1 )