油井流入动态及多相流动计算.ppt

第一章 油井流入动态与井筒多相流动计算,Inflow Performance Relationship & Wellbore Multiphase Flow,油井流入动态 Well Inflow Performance 井筒气液两相流基本概念 Basic Conception of Wellbore Two Phase Flow 计算气液两相垂直管流的Orkiszewski方法 Orkiszewski Correlation 计算井筒多相管流的Beggs-Brill 方法 Beggs-Brill Correlation,主要内容,自喷井生产系统,分离器 地面油嘴 井口 安全阀（海上油井） 节流器（海上油井） 井底流压Pwf 井底油层面上的压力Pwfs 平均地层压力Pr 集气管网 油罐,井筒设备：油管、封隔器、配产器； 地面设备：井口装置（又称采油树），内含有油嘴。,油井生产过程,井口到分离器的流动（ PB Psep ） 近似水平管流,油气从油藏流到井底（PrPwf） 地层中的渗流,四个基本流动过程：,通过油嘴的流动（Pt PB） 嘴流,从井底流到井口（Pwf Pt） 多相管流（泡流、段塞流、环流、雾流）,第一节 油井流入动态（IPR曲线）,油气井流入动态：在一定的油层压力下，流体（油，气，水）产量与相应的井底流压的关系，反映了油藏向该井供油气的能力。,表示产量与井底流压关系的曲线(Inflow Performance Relationship Curve),称为流入动态曲线，简称IPR曲线。,复合型：PwfPb Pe ， 单两相渗流,直线型：PwfPb,单相渗流， 牛顿流体（刚性水驱）,曲线型：a.非牛顿流体单相驱， b.PwfPb,两相流， 溶解气驱，粘弹流体。,IPR曲线的影响因素： 油藏驱动类型；完井状况；油藏及流体物性。,一、单相流体流入动态,供给边缘压力不变圆形地层中心一口井的产量公式为：,圆形封闭油藏，拟稳态条件下的油井产量公式为：,1.符合线性渗流规律时的流入动态,（一）单相液体流入动态,非圆形封闭地层, 油井产量公式可根据泄油面积形状和油井位置进行校正。,单相流动时，油层物性及流体性质K、o、Bo、re、rw、S基本不随压力变化:,分析：qo与Pwf的关系？,其中,单位生产压差下的油井产油(液)量，反映油层物性、流体物性、完井条件及泄油面积等与生产有关的综合指标，m3/(d. Mpa)。,采油(液)指数的定义,采油指数J的确定：,根据试井资料绘制,根据油藏参数计算：,改变油井工作制度，当油井稳定生产后，测得一系列产量与流压，绘制IPR曲线，IPR曲线斜率的负倒数即为采油指数。,注意问题：,对于非直线型IPR曲线，由于斜率不是定值，按不同定义求得的采油指数不同。 在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。,反映了油藏向井的供油能力； 反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响； 通过油井流入动态研究为油藏工程提供检验资料； 为采油工程的下一步工作提供依据； 检查钻井、固井、完井和各项工艺措施等技术水平的优劣。,IPR曲线的应用：,例：A井位于正方形泄油面积的中心，,根据系统试井，计算,例题,1）绘制IPR曲线,解：,2）计算采油指数,3）查表得,4）直线外推至q=0,求得,5）Pwf=8.8MPa时,Qo=20(12-8.8)=64m3,2.符合非线性渗流规律时的流入动态,条件：油井产量很高时，在井底附近不再符合线性渗流，呈现高速非线性渗流。,胶结地层的紊流速度系数：,非胶结地层紊流速度系数：,如果试井资料在单相渗流呈现非线性渗流，可绘制 与 的关系曲线 。,由此可以看出， 与 呈线性关系。绘制的直线的斜率为D，其截距则为C。,（二）单相气体流入动态,1.符合线性渗流规律时的流入动态,条件：定压边界、圆形气层中心有一口气井稳定生产时，距井轴r处的流量为：,根据气体连续方程和状态方程，将半径r处的流量折算为标准状态下的气井产量qg。,引用假(拟)压力的概念:,所以,用数值积分法或其它方法求得拟压力 后，再求得气井产量。,在工程中常近似地用平均压力 求和Z,即,积分得,令,则,1)绘制 与 的关系曲线；,应用：,2)求出斜率D；,3)预测稳定线性渗流条件下的气井流入动态。,2.符合非线性渗流规律时的流入动态,(1)二项式方程,矿场上将 与 的关系曲线称为二项式特征曲线。,利用上式求得一组不同 下的 ，即可绘出气井的流入动态曲线。,将 = =1.03105Pa时的产量称为气井的绝对无阻流量,用途： 衡量气井生产能力的大小； 进行气井间生产能力的比较； 是气井配产的依据。,C 与气层及流体性质有关的系数； n渗流指数，0.51。,利用试井资料求出C，n值，求得其绝对无阻流量,与液体相比，气体具有更大的压缩性，应注意问题：,气体的产量均指标准状态下的产量,气井的压力均采用绝对压力而非表压力,（2）指数式方程,二、油气两相渗流时的流入动态 Two Phase Flow IPR,o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制IPR曲线。,平面径向流，直井油气两相渗流时油井产量公式为：,(一)垂直井油气两相渗流时的流入动态,1968年，沃格尔对不同流体性质、油气比、相对渗透率、井距、压裂井、污染井等各种情况下的21个溶解气驱油藏进行了计算。,1.Vogel method(沃格尔方法),除高粘度原油及严重污染的油井外，绘制的 与 具有相似性，这一参考曲线被称为沃格尔曲线。,Vogel曲线及方程,对某一油田而言，在一定时间内是常数，可通过测试数据确定。,若测得两种不同工作制度下的油井产量，则有：,a.计算,c.根据给定的流压及计算的相应产量绘制IPR曲线。,b.给定不同流压，计算相应的产量：,已知地层压力和一个测试点：,利用Vogel方程绘制IPR曲线的步骤,油藏压力未知，已知两个测试点,a.确定油藏平均压力,b.计算,d.根据给定的流压及计算的相应产量绘制IPR曲线,c.给定不同流压，计算相应的产量,2.费特柯维奇方法Fetkovichs Approximation,油气两相渗流：,假设 与压力 成直线关系，则：,令：,当 时：,所以：,3.非完善井Vogel方程的修正,油水井的非完善性： 打开性质不完善；如射孔完成 打开程度不完善；如未全部钻穿油层 打开程度和打开性质双重不完善 油层受到损害 酸化、压裂等措施,非完善井:,非完善井附加压力降:,油井的流动效率Flow Efficiency FE：,油井的理想生产压差与实际生产压差之比。,油井不完善性的确定,利用流动效率计算非完善直井流入动态,无因次IPR曲线, Standing方法(FE=0.51.5),关键是求理想井底流压,Standing方法不完善井IPR曲线的步骤：,b.预测不同流压下的产量,c.根据计算结果绘制IPR曲线,Harrison无因次IPR曲线(FE1), Harrison方法,查FE曲线,IPR曲线 （FE一定）,当FE1.5时，应用Standing方法计算结果不合理，可应用Harrison方法：,Pwf=0,三种图版对应的油井流动效率范围不同 Harrison方法和Standing方法图版中的最大产液量是FE=1时的最大产液量，不是油井实际的最大产液量。 Harrison方法的图版可以获得高流动效率井和低流压井的最大产液量，而Standing方法不能。,VogelStandingHarrison的区别,(二)溶解气驱油藏斜井和水平井的IPR曲线,1.Cheng方法：,P=Pwf/Pr； q=qo/qomax； A、B、C为取决于井斜角的系数。,2.Bendakhlia方法,其它水平井产能计算模型：,Borisov模型：,Giger模型：,Joshi模型：,Renard & Dupuy模型：,组合型IPR曲线,三、组合型IPR曲线,油藏中同时存在单相流、油气两相流。,3.当 后，油藏中出现两相流动。,直线段 采油指数,？,，直线与曲线的斜率相等：,组合型IPR曲线的应用,当测试压力 ，求 时的产量并绘制IPR曲线,求不同qwf下的产量qo：,绘制IPR曲线,时,时,当测试压力 ，求不同 下的产量并绘制IPR曲线,求J：由,求qc,求产量qo：,绘制IPR曲线,单相流采油指数：,时,时,已知D井平均地层压力为16MPa，Pb为13MPa,Pwf为8MPa时的产量为71.45m3/d，试计算Pwf为14MPa和7MPa时的产量并绘制该井的 IPR曲线。,a.计算 (注意工作点的位置),例 题,b.计算,c.计算指定流压下的产量 (注意Pwf与Pb的关系),e.根据给定的流压及计算的相应产量绘制IPR曲线。,d.不同流压下，计算相应的产量为：,四、油气水三相流入动态(自学),Petrobras提出了计算三相流动IPR曲线的方法。,综合IPR曲线的实质: 按含水率取纯油IPR曲线和水IPR曲线的加权平均值。,当已知测试点计算采液指数时，是按产量加权平均；,当预测产量或流压时,是按流压加权平均。,五、多层油藏油井流入动态,1.层间无干扰，全井IPR曲线为各层IPR曲线的迭加,多层油藏油井流入动态,随着流压的降低，供液的小层增加、产量增加、采液指数增加。,2.层间有干扰,油层压力小于水层压力,a.全井IPR曲线 b.油层单独测试的IPR曲线 c.水层单独测试的IPR曲线,A油层静压Po B全井静压Pe C水层静压Pw,当井底流压大于全井静压(PwfPe)时，井口不出液，而此时水层出水，出水量为CB曲线。,水去哪了？,当PePwfPo时，井口出水，但水层出水量高于井口产水量。,水又去哪了？,当井底流压小于油层静压(PwfPo)时，井口同时产出水和油，此时全井的IPR曲线为油、水层IPR曲线的迭加。,在PwfPo之前，井口全部产水，油井含水率为100。,当PwfPo时，井口开始产油，含水率开始下降。,结论:,1,2,3,时，发生高压层向低压层的倒灌层间干扰；,时，无层间干扰,井口含水率随流压的降低而降低，即降低流压既可以提高油井产量，又可以降低含水率（但产水量也增加）。,油层压力大于水层压力,1,2,3,油井含水率随流压的降低而上升;,时，只产油；,时，同时产出水和油;,本 节 小 结,基本概念：油井流入动态；采液指数；表皮系数；流动效率,基本方法：IPR曲线求取,基本方法的扩展与修正：,液相,气液两相,液相+气液,油气水三相,采液指数,Vogel,组合式,Petrobras,Standing Harrison 水平井、定向井,表皮效应:由于钻井、完井、作业或采取增产措施，使井底附近地层的渗透率变差或变好，从而引起附加流动压力的效应。,流入动态研究方法,作 业,第二节 井筒气液两相流基本概念,气-液混合物在垂直管中的流动为第二个流动过程，其多相流动规律研究是进行自喷、气举、抽油井生产分析的基础。,研究目的：确定沿程压力损失的流动规律,研究途径： 基本流动方程 实验资料相关因次分析 近似关系,Basic Conception of Wellbore Two Phase Flow,一、井筒气液两相流动的特性,(一)气液两相流动与单相液流的比较,a.出现条件,单相,单相与多相共存,全井多相,多相流,油（或+水）+气,单相流,油（或+水）,(二)气液混合物在垂直管中的流动结构变化,1.油气混合物在垂直管中的流动特征,流压：从油层流到井底后具有的压力。,油压：流压作用下，克服静液柱压力和流动 阻力后的压力。,多相,垂直管流,+,混合物的密度、粘度、流速等会随温度、压力变化而变化 流动形态会发生变化,能量损失主要包括： 重力损失 摩擦损失 滑脱损失,油气混合物在垂直管中的流动特征,2.气液混合物在垂直管中的流动结构变化,从井底到井口，自喷井内可能出现的流态包括纯油(液)流、泡流、段塞流、环流和雾流。 实际上，在同一口井内，一般不会出现完整的流态变化。,原油从井底流到井口，主要克服重力损失和摩擦损失，压力越来越小，当井筒压力低于饱和压力时，溶解气从油中分离出来，气体体积不断膨胀，导致油气分布状态发生变化。, 纯液流 Liquid Flow,井筒压力低于饱和压力，溶解气从油中分离出来，气体以小气泡分散在液相中。,泡流 Bubble Flow,井筒压力大于饱和压力,单相液流，重力损失为主。,特点：气体是分散相，液体是连续相；重力损失为主，滑脱现象比较严重。,各流态的主要特征,泡流,混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡合并成大气泡，井筒内将形成一段液一段气的结构。,特点： 气体呈分散相，液体呈连续相； 一段气一段液交替出现； 气体膨胀能得到较好的利用； 滑脱损失变小； 摩擦损失变大。,段塞流 Slug Flow,段塞流,环流 Annular Flow,特点： 气液两相都是连续相； 气体举油作用主要是靠摩擦携带； 滑脱损失变小； 摩擦损失变大。,压力继续下降，气泡从中间突破，油管中心是连续的气流，管壁为油环。,环流,压力进一步下降，气体体积流量足够大，占据整个油管，液体以液滴形式分散在气流中。,特点： 气体是连续相，液体是分散相； 气体以很高的速度携带液滴喷出井口； 气、液之间的相对运动速度很小； 气相是整个流动的控制因素。,雾流 Mist Flow,雾流,(三)滑脱现象及特性参数,在垂直管流中，由于流体间密度差异，产生气体超越液体流动的现象。 (Slip Effect),1.滑脱现象,井筒气液两相流动中，混合物的重力损失一般大于其它能量损失； 重力损失的大小直接取决于井深和混合物密度； 混合物的密度与滑脱现象有关。,气液两相过流断面示意图,因滑脱而产生的附加压力损失称为滑脱损失。,通常用有、无滑脱时混合物的密度之差 来表示单位管长上的滑脱损失：,滑脱损失的实质： 液相的流动断面增大引起混合物密度的增加。,在气液两相流动状态下，液相所占单位管段容积的份额，称为持液率，也称液相存容比。 (Liquid Holdup) 气相所占单位管段容积的份额，称为空隙率，也称气相存容比。 (Gas Holdup),2.求取混合物密度,无滑脱持液率（体积含液率）：,无滑脱混合物密度：,无滑脱时的混合物密度,滑脱时的混合物密度,滑脱损失,总结,多相垂直管流特性,流态变化,压力损失,纯油流,泡流,段塞流,环流,雾流,重力损失：,摩擦损失：,滑脱损失：,计算难度增大,关键是计算,倾斜管流能量平衡关系示意图,二、井筒气液两相流能量平衡方程 及压力分布计算步骤,1.能量平衡方程推导,质量为m的流体，在一定状态下具有的能量包括：内能、位能、动能和压缩（或膨胀）能。,断面1和2流体能量平衡关系：,写成微分形式：,由于,（摩擦产生的热量消耗的功）,m取单位质量1,写成压力梯度形式,取z轴方向为自上而下：,水平管流:,垂直管流:,多相混合物流动:,研究流动过程中混合物密度、速度、摩擦系数的变化规律和计算方法是研究多相管流的中心问题。,以计算段下端压力为起点，重复步。,2.垂直多相管流压力分布计算步骤,重复的计算，直至 。,按深度增量迭代的步骤,已知任一点(井口或井底)的压力作为起点，任选一个合适的压力间隔p。,估计一个对应的深度增量h 。,计算该管段的平均温度及平均压力，并确定流体性质参数。,并计算该段的压力梯度dp/dh。,计算对应于的该段管长(深度差)h。,计算该段下端对应的深度及压力。,P0=Pwf,h,P1=P0+P,h,Pt,以计算段下端压力为起点，重复步。,重复的计算，直至 。,按压力增量迭代的步骤（自学）,已知任一点(井口或井底)的压力作为起点，以固定的h 将井筒分为n段。,估计一个对应的压力增量p。,计算该管段的平均温度及平均压力，并确定流体性质参数。,并计算该段的压力梯度dp/dh。,计算对应于的该段的压降p。,计算该段下端对应的深度及压力。,P0=Pwf,p,P1=P0+P,p,Pt,注意的问题： a.计算压力分布过程中，温度和压力是相关的； b.流体物性参数计算至关重要； c.不同的多相流计算方法差别较大，因此在实际应用中有必要根据油井的实际情况筛选精度相对高的方法。,第三节 Orkiszewski方法,综合了Griffith & Wallis 和 Duns & Ros 等方法；,处理过渡性流型时，采用Ros方法(内插法)；,针对每种流动型态提出存容比及摩擦损失的计算方法；,提出了四种流型，即泡流、段塞流、过渡流及环雾流。,把Griffith段塞流相关式改进后推广到了高流速区；,1967年奥其斯泽斯基提出了适用于多相垂直管流计算方法，其主要特点包括：,过渡流 Transition Flow,只有在雾流条件，气体体积流量远大于液体体积流量，近似等于总的体积流量，其它条件动能变化近似为0。,一、压力降公式及流动型态划分界限Pressure Gradient Correlation and Flow pattern Transitions,压力降是摩擦能量损失、势能变化和动能变化之和：,根据气体定律：,又由于,且,所以,把梯度变为计算管段，则：,计算管段的压降公式为,Orkiszewski方法流型划分界限,不同流动型态下 和 的计算方法不同。,二、平均密度及摩擦损失梯度的计算,气相存容比(含气率)Hg ：管段中气相体积与管段容积之比值。 液相存容比(持液率)HL ：管段中液相体积与管段容积之比值。,1.泡流,平均密度:,可由滑脱速度来确定,滑脱速度：气相流速与液相流速之差。,则：,泡流摩擦损失梯度按液相进行计算：,一般油管的绝对粗糙度为,2.段塞流,平均密度:,段塞流的摩擦梯度：,段塞流计算中，关键是滑脱速度vs的计算。 目前，vs的计算方法有两种：查图迭代法和经验公式法。,流体分布系数,泡流雷诺数：,C1Nb曲线,雷诺数：,滑脱速度的计算迭代法,滑脱速度：,C2NRe曲线,（假设Vs）,滑脱速度的计算经验公式法,过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。,(3)过渡流,雾流混合物平均密度计算公式与泡流相同：,由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。,雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。,摩擦梯度:,(4)雾流,所以：,Orkiszewski方法计算流程框图,第四节 Beggs & Brill 方法,水和空气、聚丙烯管实验基础上总结的方法,建立流型分布图，将七种流型归为三类，增加了过渡流,计算时先按水平管流计算，然后采用倾斜校正系数校 正成相应的倾斜管流,倾斜度 -90 +90，分上坡和下坡流动,1973年提出，适用于水平、垂直和任意倾斜管流计算,Beggs & Brill 两相水