复试采油课件第二章

第2章油井流入动态与井筒多相流动2.1油井流入动态（1）单相液体流入动态；（2）油气两相流入动态；（3）组合型流入动态；（4）油气水三相流入动态；（5）斜井和水平井流入动态；（6）多层油藏流入动态2.2井筒气液多相流（1）Orkiszewski方法；（2）Beggs-Brill方法2．1油井流入动态油井流入动态：油井产量与井底流动压力的关系。物理意义：反映了油藏向井的供油能力，反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响，是采油工程与油藏工程的衔接点。作用：为油藏工程提供检验资料为采油工程的下一步工作提供依据检查钻井、完井和各项工艺措施等技术水平的优劣2．1油井流入动态油井流入动态曲线（IPR曲线）表示产量与井底流压关系的曲线，简称IPR曲线（InflowPerformanceRelationship

）。水驱与溶解气驱聚合物驱流入动态曲线参考文献：（1）李福军等，采用聚合物驱油井流入动态研究，钻采工艺，1997年第20卷第2期（2）郑俊德等，聚合物驱采油井ipr曲线确定方法，石油钻采工艺，1998年第20卷第5期（3）夏惠芬等，聚合物驱油井地层流体流入动态关系研究，钻采工艺，1998年第21卷第3期2．1油井流入动态2．1油井流入动态气井流入动态曲线参考文献：（1）陈元千，气井新的无因次ipr方程及应用，油气井测试，1998年第7卷第4期（2）卢德唐等，气井的瞬时ipr曲线，天然气工业，2003年第23卷第6期（3）雷振中，气井无因次ipr关系式的推导及图版绘制，钻采工艺，1996年第19卷第1期2．1油井流入动态一种新型的油井流入动态曲线参考文献：（1）苏永新，一种新型的油井流入动态曲线，油气井测试，2001年第10卷第5期（2）郭冀义等，油井ipr曲线倒转现象分析，油气井测试，1997年第6卷第3期（3）王俊魁等，油井流入动态曲线与合理井底压力的确定，新疆石油地质，1999年第20卷第5期（4）夏慧芬等，油井ipr曲线拐点分析及计算，大庆石油地质与开发，1996年第15卷第2期2．1油井流入动态IPR发展历程

IPR最初只是经验地描述了油井产量与给定平均地层压力、井底流压之间的相互作用和影响。常规IPR曲线是基于Darcy线性定律,其合理应用的前提是采油指数保持不变。（1）1942年，Evinger和Muskat通过对渗流方程研究指出,当在油藏中存在两相渗流时产量与压力将不会像期望的那样存在直线关系,而是一种曲线关系。早期诸多研究油井工作情况的研究者中的两位。2．1油井流入动态IPR发展历程

（2）1968年，Vogel选用21个油田的实例数据(油藏岩石和流体性质有较大的变化范围)进行数值模拟得到一系列IPR关系数据。分析这些数据时,Vogel首先注意到这些实例的生产—压力关系曲线非常相似。他将每一个点的压力除以油藏平均压力、将每个点的产量除以油井最大产量进行无量纲化,发现这些无量纲化的IPR数据点最后落在一个狭小的范围内,经回归得到了后来称为Vogel方程的IPR曲线。2．1油井流入动态IPR发展历程

（3）1973年，Fetkovich曾经建议用油井等时试井数据来评价其生产能力,他在气井产能经验方程基础上,根据对6个油田、40口不同的油井生产数据分析结果,提出了后来称为Fetkovich方程的IPR关系式。（4）1976年，Jones、Blount和Glaze通过研究用多流量短时测试预测油井流入动态,考虑到非达西流动的影响,根据Forchheimer方程得到一种二项式IPR方程。2．1油井流入动态IPR发展历程

（5）1992年，Wiggins完成了一项非常有意义的工作,他对油气两相渗流拟稳态解式进行Tailor展开,解析得到了IPR方程一般形式。（6）1995年，Sukarno在数值模拟基础上得到了一种IPR曲线方程,试着考虑当井底流压变化时由于表皮变化(受产量变化影响)而引起的流动效率的变化。2．1油井流入动态IPR发展历程

2．1油井流入动态IPR研究方法

①利用井底压力计测量测试井的地面产量q和井底压力pwf,并建立其之间关系的井底流入动态关系式,主要有适用于未饱和油藏的直线IPR方程和描述饱和油藏中油井流入动态的Vogel(1968)方程;②利用试井方法确定油藏流入动态;③利用油藏数值计算方法来研究油藏.

上述方法有很多不同的假设条件和适用范围条件限制,最根本的假设是将油气藏视为线性动态系统,且需建立描述油气藏动态的数学模型,用算法辨识来进行其参数辨识,从而达到建立动态模型的目的。

2．1油井流入动态④利用BP神经网络技术进行油井流入动态分析方法研究将油井流入动态视为非线性动态系统,用神经网络预测油井产量随井底流压的变化情况,建立油井流入动态的神经网络模型,从而进行油井的流入动态分析;该方法不需建立描述油井流入动态的复杂数学模型,其变量的内在关系归结为相应的网络参数,只需对其动态系统的输入输出进行网络训练就可建立神经网络预测模型。参考文献：陈军斌等，利用BP网络技术进行油井流入动态分析方法研究，西安石油学院学报(自然科学版)，2002年11月，第17卷第6期IPR研究方法

2．1油井流入动态IPR研究方法

⑤利用模拟退火算法进行油井流入动态研究

Vogel曲线仅限于产水少或不产水的井，而且提出通用方程时有很多假设条件；Standing方法由于要求知道油层的体积系数、原油粘度和相对渗透率，难于应用；陈元千推导的曲线通式虽然考虑了采出程度和油井不完善性的影响，但也仅适用于低含水率的油藏；近年来,越来越多的稳定试井资料证实,注水保持压力开发的油田,当井底流压低于饱和压力以后,由于井底附近油层渗流条件发生了变化,指示曲线向压力轴偏转,产量出现最大点,此时就不能用达西公式和Vogel方程来计算油井的产量。因此,需找到一种非数值的方法进行油井的流入动态研究。2．1油井流入动态IPR研究方法

⑤利用模拟退火算法进行油井流入动态研究模拟退火算法源于对固体退火过程的模拟,固体退火是先将固体加热至熔化,再徐徐冷却使之凝固成规整晶体的热力学过程。利用模拟退火方法进行油井流入动态研究时不需要把单相流和多相流分开处理,适用范围较广;该方法无任何前提假设条件,比Vogel方法优越。参考文献：陈军斌等，利用模拟退火算法进行油井流入动态研究，数学的实践与认识，2003年7月，第33卷第7期2．1油井流入动态IPR的应用范围（作用）

①确定油气井合理产能②预测产量③制定油气井最优工作制度④分析、评价油气井动态⑤确定井底合理流压⑥确定启动压力⑦确定地层压力⑧确定地层污染情况⑨确定油井的流动效率和沃格参数及饱和压力2．1油井流入动态2.1.1单相液体渗流时的流入动态

供给边缘压力不变圆形地层中心一口井的产量公式为：圆形封闭油藏，拟稳态条件下的油井产量公式为：2．1油井流入动态2.1.1单相液体渗流时的流入动态

对于非圆形封闭泄油面积的油井产量公式，可根据泄油面积和油井位置进行校正。2．1油井流入动态2.1.1单相液体渗流时的流入动态

单相流动时，假设油层物性及流体性质基本不随压力变化，产量公式可表示为:采油（液）指数2．1油井流入动态2.1.1单相液体渗流时的流入动态

采油(液)指数：

单位生产压差下的油井产油(液)量，反映油层性质、厚度、流体物性、完井条件及泄油面积等与产量之间关系的综合指标。2．1油井流入动态2.1.1单相液体渗流时的流入动态

采油指数J的获得：

油藏参数计算试井资料：测得3～5个稳定工作制度下的产量及其流压，便可绘制该井的实测IPR曲线2．1油井流入动态2.1.1单相液体渗流时的流入动态

例：X井位于面积A=3897m2的等边三角形泄油面积的中心，井眼半径rw=0.1m，根据高压物性资料Bo=1.2，uo=3mPa.s，由压力恢复试井资料获得s=+3。测试资料如表所示。试绘制IPR曲线，并求采油指数J、油层静压Pr及油层参数koh。流压（MPa）11.1510.269.749.15产量（m3/d）17.434.145.656.82．1油井流入动态2.1.1单相液体渗流时的流入动态

（1）绘制IPR曲线2．1油井流入动态2.1.1单相液体渗流时的流入动态

（2）求采油指数及油层参数由图2-2查得（3）直线外推至q＝0处得油藏压力为12MPa2．1油井流入动态2.1.1单相液体渗流时的流入动态

对于单相液体流动的直线型IPR，采油指数可定义为产油量与生产压差之比，或者单位生产压差下的油井产油量；也可定义为每增加单位生产压差时，油井产量的增加值，或油井IPR曲线斜率的负倒数。对于多相流体流动的非直线型IPR，采油指数准确的定义为油井IPR曲线斜率的负倒数。注意事项：2．1油井流入动态2.1.1单相液体渗流时的流入动态

对于非直线型IPR曲线，由于其斜率不是定值，按上述几种定义所求得的采油指数则不同。所以，对于具有非直线型IPR曲线的油井，在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。2．1油井流入动态2.1.1单相液体渗流时的流入动态

当油井产量很高时，井底附近将出现非达西渗流：胶结地层的紊流速度系数：非胶结地层紊流速度系数：2．1油井流入动态2.1.1单相液体渗流时的流入动态

如果在单相流动条件出现非达西渗滤，也可利用试井所得的产量和压力资料求得C和D值。由试井资料绘制的～直线的斜率为D，其截距则为C。2．1油井流入动态2.1.2油气两相渗流时的流入动态

o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制IPR曲线。平面径向流，直井油气两相渗流时油井产量公式为：2．1油井流入动态2.1.2油气两相渗流时的流入动态

（1）Vogel方法①假设条件：a.圆形封闭油藏，油井位于中心；b.均质油层，含水饱和度恒定；c.忽略重力影响；d.忽略岩石和水的压缩性；e.油、气组成及平衡不变；f.油、气两相的压力相同；g.拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。2．1油井流入动态2.1.2油气两相渗流时的流入动态

（1）Vogel方法②Vogel方程2．1油井流入动态a.计算qomax

已知地层压力Pr和一个工作点（Pwf,q）：③利用Vogel方程绘制IPR曲线的步骤c.根据给定的流压及计算的相应产量绘制IPR曲线。b.给定不同流压Pwf，计算相应的产量qo：2．1油井流入动态2.1.2油气两相渗流时的流入动态

（1）Vogel方法已知B井油藏平均压力13MPa，流压11MPa时的产量为30m3/d，试用Vogel方程绘制该井的IPR曲线。（1）计算qomax（2）预测不同流压下产量流压（MPa）121110975产量（m3/d）15.630.043.355.676.893.62．1油井流入动态2.1.2油气两相渗流时的流入动态

（1）Vogel方法2．1油井流入动态油藏压力未知，已知两个工作点a.油藏平均压力的确定b.计算qomaxd.根据给定的流压及计算的相应产量绘制IPR曲线c.给定不同流压，计算相应的产量2．1油井流入动态a.按Vogel方程计算的IPR曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量。一般，误差低于5％。虽然，随着采出程度的增加，到开采末期误差上升到20％，但其绝对值却很小。④Vogel曲线与数值模拟IPR曲线的对比b.如果用测试点的资料按直线外推时，最大误差可达70～80％，只是在开采末期约30%。c.采出程度N对油井流入动态影响大，而kh/μ、B0、k、S0等参数对其影响不大。1-用测试点按直线外推2-数值模拟计算3-用Vogel方程计算2．1油井流入动态2.1.2油气两相渗流时的流入动态

（2）费特柯维奇方法溶解气驱油藏：假设与压力成直线关系，则：2．1油井流入动态2.1.2油气两相渗流时的流入动态

式中：则：令：当时：所以：（2）费特柯维奇方法2．1油井流入动态2.1.2油气两相渗流时的流入动态

（3）非完善井Vogel方程的修正油水井的非完善性：打开性质不完善；如射孔完成打开程度不完善；如未全部钻穿油层打开程度和打开性质双重不完善油层受到损害酸化、压裂等措施改变油井的完善性，从而增加或降低井底附近的压力降，影响油井流入动态关系。2．1油井流入动态2.1.2油气两相渗流时的流入动态

（3）非完善井Vogel方程的修正完善井和非完善井周围的压力分布示意图（条件：同产量）2．1油井流入动态2.1.2油气两相渗流时的流入动态

完善井:非完善井:令：非完善井附加压力降:则：（3）非完善井Vogel方程的修正表皮系数2．1油井流入动态2.1.2油气两相渗流时的流入动态

（3）非完善井Vogel方程的修正油井的流动效率（FE）：油井的理想生产压差与实际生产压差之比。油层受污染的或不完善井，完善井,增产措施后的超完善井，2．1油井流入动态2.1.2油气两相渗流时的流入动态

利用流动效率计算非完善直井流入动态的方法①Standing方法(FE=0.5~1.5)FE≠1时的无因次IPR曲线(standingIPR曲线)2．1油井流入动态Standing方法计算不完善井IPR曲线的步骤：a.根据已知Pr和Pwf计算在FE=1时最大产量c.根据计算结果绘制IPR曲线b.预测不同流压下的产量2．1油井流入动态②Harrison方法(FE=1~2.5)计算高流动效率井的IPR曲线和预测低流压下的产量2．1油井流入动态Harrison方法计算超完善井IPR曲线的步骤：a.计算FE=1时的qomax(FE=1)

先求Pwf/Pr，然后查图1-7中对应的FE曲线上的相应值qo/qomax(FE=1)。则c.计算不同流压下的产量d.根据计算结果绘制IPR曲线b.求FE对应的最大产量，即Pwf=0时的产量2．1油井流入动态2.1.3组合型流入动态（Pr>Pb>Pwf时的流入动态）

（1）基本公式当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般表达式为:2．1油井流入动态2.1.3组合型流入动态（Pr>Pb>Pwf时的流入动态）

2．1油井流入动态2.1.3组合型流入动态（Pr>Pb>Pwf时的流入动态）

（2）实用计算方法(组合型IPR方法)组合型IPR曲线2．1油井流入动态2.1.3组合型流入动态（Pr>Pb>Pwf时的流入动态）

（2）实用计算方法(组合型IPR方法)流压等于饱和压力时的产量为：①当时，由于油藏中全部为单相液体流动。流入动态公式为：②当后，油藏中出现两相流动。2．1油井流入动态2.1.3组合型流入动态（Pr>Pb>Pwf时的流入动态）

（2）实用计算方法(组合型IPR方法)直线段采油指数2．1油井流入动态2.1.4油气水三相流入动态Petrobras提出了计算三相流动IPR曲线的方法。综合IPR曲线的实质:

按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时，是按产量加权平均；当预测产量或流压时是按流压加权平均。油气水三相IPR曲线2．1油井流入动态2.1.4油气水三相流入动态(一)采液指数计算已知一个测试点(pwftest、qtest)和饱和压力Pb及油藏压力Pr①当时：②当时：③当时：其中：思考题：推导时的采液指数计算式。2．1油井流入动态2.1.4油气水三相流入动态①，则：②，则按流压加权平均进行推导：

(二)某一产量下的流压计算所以：2．1油井流入动态2.1.4油气水三相流入动态(二)某一产量下的流压计算因为：所以：③若，则综合IPR曲线的斜率可近似为常数思考题：试推导2．1油井流入动态2.1.4油气水三相流入动态1）p1>pbAndp2>pbpr=-q1*(p1-p2)/(q1-q2)+p12）p1>pbAndp2<pbAndp2>=pwfga2=1-0.2*p2/pb-0.8*(p2/pb)^2pr=((1-fw)*(-pb+pb*a2/1.8)-fw*p2+p1*q2/q1)/(q2/q1-1)程序编制：已知两个测试点(p1,q1)和(p2,q2)、qb、fw，计算IPR，分六种情况讨论。2．1油井流入动态2.1.4油气水三相流入动态3）p1<pbAndp1>=pwfgAndp2<pbAndp2>=pwfga1=1-0.2*p1/pb-0.8*(p1/pb)^2a2=1-0.2*p2/pb-0.8*(p2/pb)^2pr=(q2*(1-fw)*(-pb+pb*a1/1.8)-q2*fw*p1-q1*(1-fw)*(-pb+pb*a2/1.8)+q1*fw*p2)/(q1-q2)4）p1>=pbAndp2<pwfgpr=((p2-4/9*fw*pb)*q1/(8*fw-9)+q2*p1-4*q1*pb/9)/(q2-q1)2．1油井流入动态2.1.4油气水三相流入动态5）p1<pbAndp1>=pwfgAndp2<pwfg迭代计算：pr=pb*1.1:cdc=0Jl=(8*fw-9)*q2/(p2+(8*fw-9)*(pr-4/9*pb)-4/9*fw*pb)qb=Jl*(pr-pb)qomax=qb+Jl*pb/1.8pr_11=p1+fw*q1/Jl-0.125*(1-fw)*pb*(-1+Sqr(81-80*(q1-qb)/(qomax-qb)))pr_11=pr_11/fw2．1油井流入动态2.1.4油气水三相流入动态6）p1<pwfgAndp2<pwfgpr=(p1*q2-p2*q1-4/9*pb*(8*fw-9)*(q2-q1)-4/9*fw*pb*(q2-q1))/(8*fw-9)/(q1-q2)‘计算采液指数1）p1>=pb，PI=q1/(pr-p1)2）p1<pbAndp1>=pwfga=1-0.2*p1/pb-0.8*(p1/pb)^2PI=q1/(pr-fw*p1-(1-fw)*pb*(1-a/1.8))2．1油井流入动态2.1.4油气水三相流入动态3）p1<pwfgPI=q1*(8*fw-9)/(p1-4*pb*fw+pr*(8*fw-9)+4*pb)‘计算其他参数

qb=PI*(pr-pb)qomax=qb+PI*pb/1.8qtmax=qomax+fw*(pr-qomax/PI)*PI/(9-8*fw)2．1油井流入动态2.1.5斜井和水平井流入动态(1)常规的水平井IPR计算方法①Cheng方程

,

只需一组测试点，原始油藏压力，即可求得IPR曲线。2．1油井流入动态2.1.5斜井和水平井流入动态(1)常规的水平井IPR计算方法②Bendakhlia方程

,

Bendakhlia等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关系，得到了不同条件下的IPR曲线。结果表明：早期IPR曲线近乎于直线，随采出程度增加，曲度增加；接近衰竭时曲度减小。2．1油井流入动态2.1.5斜井和水平井流入动态(1)常规的水平井IPR计算方法②Bendakhlia方程,

2．1油井流入动态2.1.5斜井和水平井流入动态(2)水平井稳态产能方程,

①Borisov公式②Giger公式假设供给边界为椭圆，井内无阻力的条件下得到的。若水平井长度L<<reh，而椭圆形供给边界的长短轴又近似相等，即reh≈a≈b时。2．1油井流入动态2.1.5斜井和水平井流入动态(3)水平井拟稳态产能方程

,

假定油藏在所有方向上都是有封闭边界的，而且水平井的位置在矩形的有界泄油面积中是任意的。这些方法的差别在于所用数学解法和边界条件不同，但三种方法计算的产量差别很小。①Mutalik、Godbole和Joshi产能方程2．1油井流入动态2.1.5斜井和水平井流入动态(3)水平井拟稳态产能方程,

②Mutalik等人的修正公式③Economides、Brand和Frick水平井方程2．1油井流入动态2.1.5斜井和水平井流入动态(4)部分射开的水平井流入动态(5)定向井的IPR曲线缺点是方程没有归一化2．1油井流入动态新型油井流入动态注水保持压力开发的油田，当井底流压低于饱和压力以后，由于井底附近油层中渗流条件发生变化，指示曲线向压力轴偏转，并出现最大产量点2．1油井流入动态新型油井流入动态IPR曲线拐点因素分析：根据达西渗流定律，在某一瞬时油井的产量为：引入相对渗透率概念，则有：2．1油井流入动态新型油井流入动态IPR曲线拐点因素分析：（1）储层渗透率的影响在油田开发过程中,地层要发生弹性、弹－塑性和塑性变形。变形引起储层绝对渗透率和孔隙度随压力降低而减小。（2）流体相渗透率的影响压力降低到低于饱和压力以后,气体从原油中分离出来,气相饱和度Sg变大,液相饱和度SL变小。相对渗透率曲线变化规律是：气相饱和度刚开始增加时,液相相对渗透率急剧下降,但气相饱和度大于某一值(或液相饱和度小于某一值)时,液相相对渗透率下降逐渐变缓。2．1油井流入动态新型油井流入动态IPR曲线拐点因素分析：（3）湍流附加阻力影响在低压下,地层渗流中出现湍流,即存在流体的湍流附加粘度。流体内摩擦力急剧增大,地层中渗阻力急剧成倍增加,导致生产压差增大、产量却减小的可能性。使得IPR曲线出现最大产量点。2．1油井流入动态新型油井流入动态流入动态方程的建立井底附近油层中油、气、水三相流动时,油层中油的相对流动能力为：采出1t地面油时,井底条件下的油、气、水体积流量可以分别用下述各式进行计算2．1油井流入动态新型油井流入动态流入动态方程的建立如果井底压力大于饱和压力,则R=0。井底三相流动中,液相(油和水)的相对流动能力为：①当井底压力大于饱和压力,且含水率为零时,油的相对流动能力为1;②当井底压力大于饱和压力,且为油水两相流时,液相的相对流动能力亦为1;2．1油井流入动态新型油井流入动态流入动态方程的建立③当井底压力低于饱和压力,且含水率为零时,油的相对流动能力为：随油井含水率和井底气油比上升,油的相对流动能力下降;随气油比上升,液相流入能力下降;随含水率上升,液相流入能力增加。2．1油井流入动态新型油井流入动态流入动态方程的建立油的流动方程为：油水两相的流动方程为：2．1油井流入动态新型油井流入动态第一点是直线弯曲的始点,该点处的流动压力等于饱和压力,流动压力低于该点以后,采油指数降低,产量增长速度减慢;第二个特征点为最大产量点,该点对应的压力可称为油井最低允许流动压力,流动压力低于该点以后,产量开始降低主要原因：流动压力下降到一定程度以后,井底附近可能出现油气两相流动,使油相的流动能力急剧下降。生产压差对产量的贡献已经小于采油指数下降对产量所产生的影响。2．1油井流入动态新型油井流入动态油井最低允许流动压力为：2．1油井流入动态新型油井流入动态其他参数如下:Bo=1.121;α=3.97m3/(m3

·MPa);T=49K;Do=0.848t/m3;B(pb)=1.3219m3/t;pb=10.4MPa;fw=63.0%。例：某井稳定试井数据见下表2．1油井流入动态新型油井流入动态（1）计算溶解气油比（2）计算平均采液指数（3）油井流入动态方程2．1油井流入动态新型油井流入动态（4）最低允许流动压力值和最大理论产量2．1油井流入动态2.1.6多层油藏油井流入动态（1）多油层油井流入动态—迭加型IPR随着流压的降低，由于参加工作的小层数增多，产量将大幅度增加，采油指数也随之而增大。2．1油井流入动态2.1.6多层油藏油井流入动态（2）含水油井流入动态含水油井流入动态与含水变化

()

含水油井流入动态曲线

()转渗现象2．1油井流入动态小结(1)上述介绍的方法阐明了油井流入动态的物理意义，也是目前现场最常用的计算方法。(2)油井流入动态研究主要有三种途径：基于Vogel、Fetkovich、Petrobras方法的完善。建立不同类型油藏和井底条件的渗流模型。利用单井流入动态的油藏数值模拟技术。(3)油井流入动态是采油工程各项技术措施设计、分析与评价的依据。(4)作业：习题1.1、习题1.3、习题1.5、习题1.6、习题1.9。2.2井筒气液多相流井筒多相流理论：研究各种举升方式油井生产规律基本理论研究特点：流动复杂性、无严格数学解研究途径：基本流动方程实验资料相关因次分析近似关系2.2井筒气液多相流2.2.1气液混合物在垂直管中的流动特性(一)气液两相流动与单相液流的比较2.2井筒气液多相流2.2.1气液混合物在垂直管中的流动特性流动型态（流动结构、流型）：

流动过程中油、气的分布状态。影响流型的因素各相介质的体积比例介质的流速各相的物理及化学性质(密度、粘度界面张力等)流道的几何形状壁面特性管道的安装方式(二)气液混合物在垂直管中的流动结构变化流动型态的划分方法：两类第一类划分方法：根据两相介质分布的外形划分，包括泡状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流

垂直气液两相流流型

水平气液两相流流型第二类划分方法：按流动的数学模型或流体的分散程度划分，包括分散流、间歇流、分离流

两种分类方法比较

第一类划分方法较为直观第二类划分方法便于进行数学处理第一类划分方法泡状流弹状流或团状流层状流波状流段塞流或冲击流环状流雾状流第二类划分方法分散流间歇流分离流分离流间歇流分离流分散流两类划分结果的对应关系,

(二)气液混合物在垂直管中的流动结构变化②泡流当井筒压力稍低于饱和压力时，溶解气开始从油中分离出来，气体都以小气泡分散在液相中。滑脱现象：混合流体流动过程中，由于流体间的密度差异，引起的小密度流体流速大于大密度流体流速的现象。如：油气滑脱、气液滑脱、油水滑脱等。特点：气体是分散相，液体是连续相；气体主要影响混合物密度，对摩擦阻力影响不大；滑脱现象比较严重。①纯液流:当井筒压力大于饱和压力时，天然气溶解在原油中，产液呈单相液流。,

特点：气体呈分散相，液体呈连续相；一段气一段液交替出现；气体膨胀能得到较好的利用；滑脱损失变小；摩擦损失变大。③段塞流

当混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡将合并成大气泡，直到能够占据整个油管断面时，井筒内将形成一段液一段气的结构。,

特点：气液两相都是连续相；气体举油作用主要是靠摩擦携带；滑脱损失变小；摩擦损失变大。④环流

油管中心是连续的气流而管壁为油环的流动结构。,

特点：气体是连续相，液体是分散相；气体以很高的速度携带液滴喷出井口；气、液之间的相对运动速度很小；气相是整个流动的控制因素。⑤雾流

气体的体积流量增加到足够大时，油管中内流动的气流芯子将变得很粗，沿管壁流动的油环变得很薄，绝大部分油以小油滴分散在气流中。,

(二)气液混合物在垂直管中的流动结构变化总结：油井生产中可能出现的流型自下而上依次为：纯油(液)流、泡流、段塞流、环流和雾流。实际上，在同一口井内，一般不会出现完整的流型变化。,

2.2井筒气液多相流2.2.2气液混合物在水平管中的流动特性雾状流（1）泡状流：此时气体量很少，气体以小气泡的形式在管道中与液体一起流动。（2）团状流：随着气体量增多，气泡合并成较大的气团。气团在管道中与液体一起流动。（3）层状流：气体量再增多，气团联成一片。气相和液相分成具有光滑界面的气体层和液体层。（4）波状流：气体量进一步增多，流速提高，在气液界面上引起波浪。（5）冲击流：又称段塞流。气体的流速更大时，波浪加剧，并且波浪的顶部不时高达管壁的上部。（6）环状流：气体量和流速进一步提高，随着气体流速的提高，液体断面将进一步变薄，并且沿管壁搭接成环形断面。（7）雾状流：当气体流速很大时，液体被气流吹散，以液滴或雾的形式随着高速气流向前流动。,

2.2井筒气液多相流2.2.3井筒气液两相流能量平衡方程及压力分布计算

两个流动断面间的能量平衡关系：(一)能量平衡方程推导,

(一)能量平衡方程推导倾斜多相管流断面1和断面2的流体的能量平衡关系为：,

(一)能量平衡方程推导适合于各种管流的通用压力梯度方程：则：令：,

(二)多相垂直管流压力分布计算步骤（1）按深度增量迭代的步骤①已知任一点(井口或井底)的压力作为起点，任选一个合适的压力降作为计算的压力间隔p。②估计一个对应的深度增量h

。③计算该管段的平均温度及平均压力，并确定流体性质参数。④判断流型，并计算该段的压力梯度dp/dh。⑧以计算段下端压力为起点，重复②～⑦步，计算下一段的深度和压力，直到各段的累加深度等于管长为止。⑥重复②～⑤的计算，直至。⑤计算对应于的该段管长(深度差)h。⑦计算该段下端对应的深度及压力。,

(二)多相垂直管流压力分布计算步骤（2）按压力增量迭代的步骤（略）思考题：根据上述步骤整理出计算压力分布的程序流程框图。说明：a.计算压力分布过程中，温度和压力是相关的；b.流体物性参数计算至关重要，但目前方法精度差；c.不同的多相流计算方法差别较大，因此在实际应用中有必要根据油井的实际情况筛选精度相对高的方法。,

2.2井筒气液多相流2.2.4常用多相流计算方法2.4.1.1Orkiszewski方法

综合了Griffith&Wallis和Duns&Ros等方法

处理过渡性流型时，采用Ros方法(内插法)

针对每种流动型态提出存容比及摩擦损失的计算方法

提出了四种流型，即泡流、段塞流、过渡流及环雾流

把Griffith段塞流相关式改进后推广到了高流速区

1967年提出，适用于垂直管流计算,

2.2.4常用多相流计算方法2.4.1.1Orkiszewski方法出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：（一）压力降公式及流动型态划分界限由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和：所以压降计算式为：,

2.4.1.1Orkiszewski方法Orkiszewski方法流型划分界限不同流动型态下和的计算方法不同。流动型态界

限泡

流BtgLqq<段

塞

流SgBtgLvLqq<>,过

渡

流SgMLvL>>雾

流MgLv>,

气相存容比(含气率)Hg

：管段中气相体积与管段容积之比值。液相存容比(持液率)HL

：管段中液相体积与管段容积之比值。平均密度:（二）平均密度及摩擦损失梯度的计算(1)泡流滑脱速度：气相流速与液相流速之差。则：泡流摩擦损失梯度按液相进行计算：

图1-21摩擦阻力系数曲线(2)段塞流平均密度:段塞流的摩擦梯度：段塞流计算中，关键是滑脱速度vs的计算。目前，vs的计算方法有两种：查图迭代法和经验公式法。,

泡流雷诺数：雷诺数：滑脱速度的计算——迭代法滑脱速度：滑脱速度的计算——经验公式计算法,

(3)过渡流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。,

雾流混合物平均密度计算公式与泡流相同：由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得摩擦梯度:(4)雾流所以：,

2.2井筒气液多相流,

2.2井筒气液多相流

当油气集输管线穿越丘陵及多山地带或铺设在海底并向上倾斜延伸到海岸时，集输管线与水平方向存在一定的倾角，因而将水平管中气液两相流动的压差计算方法用于这种倾斜流动时常是不成功的。这是因为，倾斜管中气液两相混合物在上升时具有重位压差，且重位压差常大于其摩阻压差；同时，由于下坡流动时气液混合物的密度及持液率通常要比上坡时小得多，从而使气液两相倾斜管流上坡时的压力损失在下坡地段是难以象单相流动那样完全恢复。,

2.2井筒气液多相流2.2.4常用多相流计算方法2.4.1.2Beggs-Brill方法

水和空气、聚丙烯管实验基础上总结的方法

建立流型分布图，将七种流型归为三类，增加了过渡流

计算时先按水平管流计算，然后采用倾斜校正系数校正成相应的倾斜管流

倾斜度-90°~+90°，分上坡和下坡流动

1973年提出，适用于水平、垂直和任意倾斜管流计算2.2井筒气液多相流2.2.4常用多相流计算方法2.4.1.2Beggs-Brill方法分离流分层流波状流环状流间歇流团状流段塞流分散流泡流雾流2.4.1.2Beggs-Brill方法（一）基本方程

单位质量气液混合物稳定流动的机械能量守恒方程为：(1)位差压力梯度：消耗于混合物静水压头的压力梯度。假设条件:气液混合物既未对外作功，也未受外界功。(2)摩擦压力梯度：克服管壁流动阻力消耗的压力梯度。(3)加速度压力梯度：由于动能变化而消耗的压力梯度。,

2.4.1.2Beggs-Brill方法忽略液体压缩性、考虑到气体质量流速变化远远小于气体密度变化，则：(4)总压力梯度（Beggs-Brill方法的基本方程）,

2.4.1.2Beggs-Brill方法（二）Beggs&Brill方法的流型分布图及流型判别式Ⅰ-分离流；Ⅱ-间歇流；Ⅲ-分散流；Ⅳ-过渡流,

（二）Beggs&Brill方法的流型分布图及流型判别式Beggs-Brill法流型判别条件,

（三）持液率及混合物密度确定实验结果表明，倾斜校正系数与倾斜角、无滑脱持液率、弗洛德数及液体速度数有关。不同EL下的倾斜校正系数曲线大约都在与水平方向成±50°的地方发生翻转。重力和粘滞力对液相的综合影响(1)持液率Beggs&Brill方法计算倾斜管流时首先按水平管计算,然后进行倾斜校正。2.2井筒气液多相流2.2.4常用多相流计算方法2.4.1.2Beggs-Brill方法（三）持液率及混合物密度确定气液流量配比一定的情况下，气液滑脱程度越大，则持液率数值越大在不同的气液两相倾斜流动状况下，管道倾角的变化将引起流体的重新分布，从而使得重力和粘滞阻力对液相产生不同的影响，最终引起滑脱速度和持液率的变化2.2井筒气液多相流上坡流动：当管子角度沿正方向增加时，作用在液体上的重力使液体的流速减小，因而增加了滑脱和持液率。随着角度的进一步增加，液体在管子里全部搭接起来，减少了两相之间的滑脱，因而也减小了持液率。角度对持液率影响2.2井筒气液多相流角度对持液率影响下坡流动：当角度沿着负方向增加时，使液体的流速增加，因而减小了滑脱和持液率。随着角度沿负方向的进一步增加，有更多的液体与管壁接触，粘性的拖拽使得液体的流速减小，持液率增加。,

（三）持液率及混合物密度确定根据实验结果回归的倾斜校正系数的相关式为:对于垂直管:系数C与无滑脱持液率、弗洛德数和液相速度数有关。系数d、e、f、g,

（三）持液率及混合物密度确定对于过渡流型，先分别用分离流和间歇流计算，之后采用内插法确定其持液率。利用持液率计算流动条件下混合物实际密度：,

（四）阻力系数气液两相流阻力系数与无滑脱气液两相流阻力系数的比值与持液率和无滑脱持液率(入口体积含液率)之间的关系：当1<y<1.2时其中：两相流动的雷诺数：也可用Moody图上的光滑管曲线来确定:气液两相流阻力系数:,

2.2