油井流入动态与多相流PPT学习教案

1、会计学1油井流入动态与多相流油井流入动态与多相流油井流入动态曲线（IPR曲线）： 表示产量与井底流压关系的曲线，简称IPR曲线。油井流入动态： 油井产量与井底流动压力的关系。它反映了油藏向井的供油能力，反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响，是采油工程与油藏工程的衔接点。作用：通过油井流入动态研究为油藏工程提供检验资料；为采油工程的下一步工作提供依据；检查钻井、固井、完井和各项工艺措施等技术水平的优劣。第1页/共102页第2页/共102页IPRIPR发展历程发展历程 IPR IPR 最初只是经验地描述了油井产量与给定平均地层压力、最初只是经验地描述了油井产量与给定平

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

3、的两位。油井工作情况的研究者中的两位。 第3页/共102页IPRIPR发展历程发展历程 （2 2）19681968年，年，VogelVogel选用选用21 21 个油田的实例数据个油田的实例数据( (油藏岩石油藏岩石和流体性质有较大的变化范围和流体性质有较大的变化范围) ) 进行数值模拟得到一系列进行数值模拟得到一系列IPR IPR 关系数据。分析这些数据时关系数据。分析这些数据时,Vogel ,Vogel 首先注意到这些首先注意到这些实例的生产实例的生产压力关系曲线非常相似。他将每一个点的压力关系曲线非常相似。他将每一个点的压压力力除以油藏平均压力、将每个点的除以油藏平均压力、将每个点的产量

4、产量除以油井最大产量除以油井最大产量进行进行无量纲化无量纲化, , 发现这些无量纲化的发现这些无量纲化的IPR IPR 数据点最后落在数据点最后落在一个狭小的范围内一个狭小的范围内, , 经回归得到了后来称为经回归得到了后来称为Vogel Vogel 方程的方程的IPR IPR 曲线。曲线。 第4页/共102页IPRIPR发展历程发展历程 （3 3）19731973年，年，FetkovichFetkovich曾经建议用油井曾经建议用油井等时试井等时试井数据来数据来评价其生产能力评价其生产能力, , 他在他在气井气井产能经验方程基础上产能经验方程基础上, , 根据对根据对6 6 个油田、个油田、

5、40 40 口不同的油井生产数据分析结果口不同的油井生产数据分析结果, , 提出了提出了后来称为后来称为Fetkovich Fetkovich 方程的方程的IPRIPR关系式。关系式。 （4 4）19761976年，年，JonesJones、Blount Blount 和和Glaze Glaze 通过研究用通过研究用多流多流量短时测试量短时测试预测油井流入动态预测油井流入动态, , 考虑到考虑到非达西流动非达西流动的影响的影响, , 根据根据Forchheimer Forchheimer 方程得到一种二项式方程得到一种二项式IPR IPR 方程。方程。 第5页/共102页IPRIPR发展历程发

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

7、起的流动效率的变化。 第6页/共102页IPRIPR发展历程发展历程 第7页/共102页IPRIPR研究方法研究方法 利用井底压力计测量测试井的地面产量利用井底压力计测量测试井的地面产量q q 和井底压力和井底压力p pwfwf, , 并建立其之间关系的井底流入动态关系式并建立其之间关系的井底流入动态关系式, , 主要有适用于主要有适用于未饱和油藏的直线未饱和油藏的直线IPR IPR 方程和描述饱和油藏中油井流入动方程和描述饱和油藏中油井流入动态的态的Vogel (1968) Vogel (1968) 方程方程; ; 利用利用试井方法试井方法确定油藏流入动态确定油藏流入动态; ; 利用利用油藏

8、数值计算方法油藏数值计算方法来研究油藏来研究油藏. . 上述方法有很多不同的上述方法有很多不同的假设条件假设条件和和适用范围条件适用范围条件限制限制, , 最根本的假设是将油气藏视为最根本的假设是将油气藏视为线性动态系统线性动态系统, , 且需建立描且需建立描述油气藏动态的数学模型述油气藏动态的数学模型, ,用算法辨识来进行其参数辨识用算法辨识来进行其参数辨识, , 从而达到建立动态模型的目的。从而达到建立动态模型的目的。 第8页/共102页利用利用BPBP神经网络技术神经网络技术进行油井流入动态分析方法研究进行油井流入动态分析方法研究 将油井流入动态视为将油井流入动态视为非线性动态系统非线性

9、动态系统, , 用神经网络用神经网络预测油井产量随井底流压的变化情况预测油井产量随井底流压的变化情况, , 建立油井流入动建立油井流入动态的神经网络模型态的神经网络模型, , 从而进行油井的流入动态分析从而进行油井的流入动态分析; ; 该方法不需建立描述油井流入动态的复杂数学模型该方法不需建立描述油井流入动态的复杂数学模型, , 其变量的内在关系归结为相应的网络参数其变量的内在关系归结为相应的网络参数, , 只需对其动只需对其动态系统的输入输出进行网络训练就可建立神经网络预测态系统的输入输出进行网络训练就可建立神经网络预测模型。模型。参考文献：陈军斌等，利用参考文献：陈军斌等，利用BP BP

10、网络技术进行油井流入动态分析方网络技术进行油井流入动态分析方法研究，西安石油学院学报法研究，西安石油学院学报( (自然科学版自然科学版) )，2002 2002 年年11 11 月，第月，第1717卷卷第第6 6期期 IPRIPR研究方法研究方法 第9页/共102页IPRIPR研究方法研究方法 利用利用模拟退火算法模拟退火算法进行油井流入动态研究进行油井流入动态研究 VogelVogel曲线曲线仅限于产水少或不产水的井，而且提出通用仅限于产水少或不产水的井，而且提出通用方程时有很多假设条件；方程时有很多假设条件；StandingStanding方法方法由于要求知道油层由于要求知道油层的体积系数

11、、原油粘度和相对渗透率，难于应用；的体积系数、原油粘度和相对渗透率，难于应用；陈元千陈元千推导的曲线通式虽然考虑了采出程度和油井不完善性的影推导的曲线通式虽然考虑了采出程度和油井不完善性的影响，但也仅适用于低含水率的油藏；响，但也仅适用于低含水率的油藏； 近年来近年来, , 越来越多的稳定试井资料证实越来越多的稳定试井资料证实, , 注水保持压注水保持压力开发的油田力开发的油田, ,当当井底流压低于饱和压力井底流压低于饱和压力以后以后, , 由于井底由于井底附近油层渗流条件发生了变化附近油层渗流条件发生了变化, , 指示曲线向压力轴偏转指示曲线向压力轴偏转, , 产量出现最大点产量出现最大点,

12、 , 此时就不能用达西公式和此时就不能用达西公式和Vogel Vogel 方程来方程来计算油井的产量。因此计算油井的产量。因此, , 需找到一种非数值的方法进行油需找到一种非数值的方法进行油井的流入动态研究。井的流入动态研究。第10页/共102页IPRIPR研究方法研究方法 利用利用模拟退火算法模拟退火算法进行油井流入动态研究进行油井流入动态研究 模拟退火算法源于对固体退火过程的模拟模拟退火算法源于对固体退火过程的模拟, , 固体退固体退火是先将固体加热至熔化火是先将固体加热至熔化, , 再徐徐冷却使之凝固成规整再徐徐冷却使之凝固成规整晶体的热力学过程。晶体的热力学过程。 利用模拟退火方法进行

13、油井流入动态研究时利用模拟退火方法进行油井流入动态研究时不需要不需要把单相流和多相流分开处理把单相流和多相流分开处理, , 适用范围较广适用范围较广; ; 该方法该方法无无任何前提假设条件任何前提假设条件, , 比比VogelVogel方法优越。方法优越。 参考文献：陈军斌等，利用模拟退火算法进行油井流入动态研究，参考文献：陈军斌等，利用模拟退火算法进行油井流入动态研究，数学的实践与认识，数学的实践与认识，20032003年年7 7月，第月，第3333卷第卷第7 7期期第11页/共102页IPRIPR的应用范围（作用）的应用范围（作用） 确定油气井合理产能确定油气井合理产能预测产量预测产量制定

14、油气井最优工作制度制定油气井最优工作制度分析、评价油气井动态分析、评价油气井动态确定井底合理流压确定井底合理流压确定启动压力确定启动压力确定地层压力确定地层压力确定地层污染情况确定地层污染情况确定油井的流动效率和沃格参数及饱和压力确定油井的流动效率和沃格参数及饱和压力第12页/共102页一、 单相液体流入动态(基于达西定律)asrrBPPhkqweoowfroo21ln)(2供给边缘压力不变圆形地层中心一口井的产量公式为：（1-1）asrrBPPhkqweoowfroo43ln)(2圆形封闭油藏，拟稳态条件下的油井产量公式为：（1-2）第13页/共102页图1-2 泄油面积形状与油井的位置系数

15、对于非圆形封闭泄油面积的油井产量公式，可根据泄油面积和油井位置进行校正。Xrrwe第14页/共102页sXBhakJooo43ln2sXBhakJooo21ln2)(wfroPPqJ采油采油( (液液) )指数指数： 单位生产压差下的油井产油(液)量，反映油层性质、厚度、流体物性、完井条件及泄油面积等与产量之间关系的综合指标。单相流动时，油层物性及流体性质基本不随压力变化，产量公式可表示为:)(wfroPPJq第15页/共102页采油指数J的获得： 油藏参数计算 试井资料：测得35个稳定工作制度下的产量及其流压，便可绘制该井的实测IPR曲线 对于单相液体流动的直线型IPR曲线，采油指数可定义为

16、产油量与生产压差之比，或者单位生产压差下的油井产油量；也可定义为每增加单位生产压差时，油井产量的增加值，或油井IPR曲线斜率的负倒数。注意事项： 对于多相流动的非直线型IPR曲线，由于其斜率不是定值，按上述几种定义所求得的采油指数则不同。所以，对于具有非直线型IPR曲线的油井，在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。第16页/共102页2DqCqPPwfrhakSxBCooo2)43(lnworhBD222134103396. 1201. 1710906. 1k55. 061008. 1kg当油井产量很高时井底附近将出现非达西渗流非

17、达西渗流：胶结地层的紊流速度系数：非胶结地层紊流速度系数：第17页/共102页二、油气两相渗流时的流入动态o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制通常结合生产资料来绘制IPRIPR曲线。曲线。平面径向流，直井油气两相渗流时油井产量公式为：(一)垂直井油气两相渗流时的流入动态drdpBhrkqoooo2kkkorodpBKrrkhqewfPPooroweoln2第18页/共102页1.1.Vogel Vogel 方法方法 假设条件：a.圆形封闭油藏，油井位于中心；b.均质油层，含水饱和度恒定；c.忽略重力影响；d.忽略岩石和水的压缩性； e.油、气

18、组成及平衡不变；f.油、气两相的压力相同；g.拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。第19页/共102页2max8 . 02 . 01rwfrwfooPPPPqqVogel方程图1-3 Vogel 曲线第20页/共102页8 . 02 . 01 2maxrtestwfrtestwftestooPPPPqqmaxoq a.计算max28 . 02 . 01orwfrwfoqPPPPqc.根据给定的流压及计算的相应产量绘制IPR曲线。b.给定不同流压，计算相应的产量： 已知地层压力和一个工作点：已知地层压力和一个工作点：利用Vogel方程绘制IPR曲线的步骤第21页/共102页A

19、ACBBPr242121qqA12212 . 0wfwfPPqqB2122218 . 0wfwfPPqqC 油藏压力未知，已知两个工作点油藏压力未知，已知两个工作点a. 油藏平均压力的确定maxoq b.计算d.根据给定的流压及计算的相应产量绘制IPR曲线c. 给定不同流压，计算相应的产量第22页/共102页Vogel曲线与数值模拟IPR曲线的对比a.按Vogel方程计算的IPR曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量。一般，误差低于5。虽然，随着采出程度的增加，到开采末期误差上升到20，但其绝对值却很小。b.如果用测试点的资料按直线外推时，最大误差可达 7080，只是在开采

20、末期约30%。C.采出程度N对油井流入动态影响大，而kh/、B0、k、S0等参数对其影响不大。图1-4 不同方法计算的油井IPR曲线1-用测试点按直线外推；2-计算机计算的；3-用Vogel方程计算的第23页/共102页dpBksrrkhqrwfPPooroweo43ln22.费特柯维奇方法溶解气驱油藏：假设 与压力 成直线关系，则：ooroBkprwfPPweocpdpsrrkhq43ln222243ln2wfrwePPcsrrkh第24页/共102页rpoororKPc)(1式中：rwfrPooroweoPPPKsrrkhqr243ln222则：rpooroweoPKsrrkhJr2143

21、ln2令：当 时：0wfP243ln2maxrPooroweoPKsrrkhqr所以：2max1rwfooPPqq)(22wfroPPJ第25页/共102页3.非完善井Vogel方程的修正油水井的非完善性： 打开性质不完善；如射孔完成 打开程度不完善；如未全部钻穿油层 打开程度和打开性质双重不完善 油层受到损害 酸化、压裂等措施改变油井的完善性，从而增加或降低井底附近的压力降，影响油井流入动态关系。第26页/共102页 图1-5 完善井和非完善井周围 的压力分布示意图第27页/共102页完善井:weoowfeoorrlnB)PP(hk2q非完善井:wssseoowfeorrkrrkBPPhql

22、n1ln1)(20令：wssorrkksln1shkBqPoooosk2非完善井附加压力降:wssooooowfwfskrrkkhkBqPPPln12则：第28页/共102页油井的流动效率（FE）：油井的理想生产压差与实际生产压差之比。0s1FE油层受污染的或不完善井，0s1FE完善井,0s1FE增产措施后的超完善井，shkBqPoooosk2wfrskwfrwfrwfrpppppppppFE第29页/共102页利用流动效率计算非完善直井流入动态的方法图1-6 Standing 无因次IPR曲线 Standing方法(FE=0.51.5)，扩大了 Vogel的使用范围，可以适用于那些污阻井或经

23、过增产措施的井第30页/共102页2)1max(8 . 02 . 01rwfrwfoFEoPPPPqqFEPPPPwfrrwf)(a.根据已知Pr和Pwf计算在FE=1时最大产量21max8 . 02 . 01rwfrwfFEooPPPPqqb.预测不同流压下的产量c.根据计算结果绘制IPR曲线第31页/共102页 Harrison方法 提供了FE=1 2.5的无因次IPR曲线,扩大了Standing曲线的范围,它可用来计算高流动效率井的IPR曲线和预测低流压下的产量。图1-7 Harrison无因次IPR曲线(FE1) )1(maxFEqqoo图 2-7 Harrison 无 因 次 IPR

24、 曲 线 (FE1)第32页/共102页1maxmax1maxmaxFEoFEoFEoFEoqqqqb.求FE对应的最大产量，即Pwf=0时的产量a.计算FE=1时的qomax(FE=1) 先求Pwf/Pr，然后查图1-7中对应的FE曲线上的相应值qo/qomax(FE=1)。则)1max()1max(/FEoooFEoqqqq1max1maxFEooFEooqqqqc.计算不同流压下的产量d.根据计算结果绘制IPR曲线第33页/共102页(1)(1)常规的水平井常规的水平井IPRIPR计算方法计算方法ChengCheng方程方程 , 4.斜井和水平井的IPR曲线2PCPBAq Cheng对溶

25、解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似Vogel方程的不同井斜角井的IPR回归方程：P=Pwf/Pr； q=qo/qomax ；A、B、C为取决于井斜角的系数。第34页/共102页 BendakhliaBendakhlia方程方程 , BendakhliaBendakhlia等用两种三维三相等用两种三维三相黑油模拟器黑油模拟器研究了多种情况下溶研究了多种情况下溶解气驱油藏中解气驱油藏中水平井水平井的流入动态关系，得到了不同条件下的的流入动态关系，得到了不同条件下的IPRIPR曲线。结果表明：早期曲线。结果表明：早期IPRIPR曲线近乎于直线，随采出程度增加，曲线近乎于

26、直线，随采出程度增加，曲度增加；接近衰竭时曲度减小。曲度增加；接近衰竭时曲度减小。nrwfrwfooPPvPPvqq2max11第35页/共102页(2)(2)水平井稳态产能方程水平井稳态产能方程 , BorisovBorisov公式公式 wehoohrhLhLrBpKhq2ln4ln)/(2GigerGiger公式公式)2(ln2/211ln)/(22wehehoohrhrLrLhLBpKLq假设供给边界为椭圆，井内无阻力的条件下得到的。假设供给边界为椭圆，井内无阻力的条件下得到的。若水平井长度若水平井长度LrLreheh，而椭圆形供给边界的长短轴，而椭圆形供给边界的长短轴又近似相等，即又近

27、似相等，即rehab时。时。 第36页/共102页(3)(3)水平井拟稳态产能方程水平井拟稳态产能方程 , 假定油藏在所有方向上都是有封闭边界的，而且水平井的位置在矩形的有界泄油面积中是任意的。这些方法的差别在于所用数学解法和边界条件不同，但三种方法计算的产量差别很小。 Mutalik、Godbole和和Joshi产能方程产能方程 DqCSSSArrBKhppqJhcAmfweoowfrhh ,ln)/(2ArewfrLS4lnMutalik等人的修正公式等人的修正公式 DqSLhKKCSSrrBphKKqzxhcAfweooyxh,738. 0ln2DqSLxpBppKxqodoowfre2

28、22.877Economides、Brand和和Frick水平井方程水平井方程第37页/共102页(4)(4)部分射开的水平井流入动态部分射开的水平井流入动态 mWDoohSpBKhJ*)/(2(5)(5)定向井的定向井的IPRIPR曲线曲线 2PCPBAqrwfppp maxooqqq 缺点是方程没有归一化缺点是方程没有归一化 第38页/共102页5.5.新型油井流入动态新型油井流入动态注水保持压力开注水保持压力开发的油田，发的油田，当井当井底流压低于饱和底流压低于饱和压力以后压力以后，由于，由于井底附近油层中井底附近油层中渗流条件发生变渗流条件发生变化，指示曲线向化，指示曲线向压力轴偏转，

29、并压力轴偏转，并出现出现最大产量点最大产量点第39页/共102页IPRIPR曲线拐点因素分析曲线拐点因素分析：根据达西渗流定律，在某一瞬时油井的产量为：根据达西渗流定律，在某一瞬时油井的产量为：引入相对渗透率概念，则有：引入相对渗透率概念，则有：第40页/共102页IPRIPR曲线拐点因素分析：曲线拐点因素分析：（1 1）储层渗透率的影响）储层渗透率的影响：在油田开发过程中：在油田开发过程中, ,地层要发生弹性、地层要发生弹性、弹塑性和塑性变形。变形引起储层绝对渗透率和孔隙度随压弹塑性和塑性变形。变形引起储层绝对渗透率和孔隙度随压力降低而减小。力降低而减小。 （2 2）流体相渗透率的影响）流体

30、相渗透率的影响： 压力降低到低于饱和压力以后压力降低到低于饱和压力以后, , 气体从原油中分离出来气体从原油中分离出来, , 气相饱和度气相饱和度S Sg g 变大变大, , 液相饱和度液相饱和度S SL L 变小。相对渗透率曲线变化规律是变小。相对渗透率曲线变化规律是：气相饱和度刚开始增加时气相饱和度刚开始增加时, ,液相相对渗透率急剧下降液相相对渗透率急剧下降, ,但气相饱和度大于某一值但气相饱和度大于某一值( (或液相饱或液相饱和度小于某一值和度小于某一值) ) 时时, ,液相相对渗透率下降逐渐变缓。液相相对渗透率下降逐渐变缓。 （3 3）湍流附加阻力影响）湍流附加阻力影响：在低压下：在

31、低压下, ,地层渗流中出现湍流地层渗流中出现湍流, ,即即存在流体的湍流附加粘度。流体内摩擦力急剧增大存在流体的湍流附加粘度。流体内摩擦力急剧增大, ,地层中渗地层中渗阻力急剧成倍增加阻力急剧成倍增加, ,导致生产压差增大、产量却减小的可能性。导致生产压差增大、产量却减小的可能性。使得使得IPR IPR 曲线出现最大产量点。曲线出现最大产量点。 第41页/共102页流入动态方程的建立流入动态方程的建立 井底附近油层中油、气、水三相流动时井底附近油层中油、气、水三相流动时, , 油层中油的相对流油层中油的相对流动能力为：动能力为：采出采出1 t 1 t 地面油时地面油时, ,井底条件下的油、气、

32、水体积流量井底条件下的油、气、水体积流量可以分别用下述各式进行计算可以分别用下述各式进行计算如果井底压力大于饱和压力如果井底压力大于饱和压力, ,则则R R = 0 = 0 。第42页/共102页当井底压力低于饱和压力当井底压力低于饱和压力, ,且含水率为零时且含水率为零时, ,油的相对油的相对流动能力为：流动能力为：井底三相流动中井底三相流动中, ,液相液相( (油和水油和水) ) 的相对流动能力为：的相对流动能力为： 当井底压力大于饱和压力当井底压力大于饱和压力, ,且含水率为零时且含水率为零时, ,油的相对流油的相对流动能力为动能力为1 ;1 ; 当井底压力大于饱和压力当井底压力大于饱和

33、压力, ,且为油水两相流时且为油水两相流时, ,液相的相液相的相对流动能力亦为对流动能力亦为1 ;1 ;第43页/共102页油的流动方程为：油的流动方程为：油水两相的流动方程为：油水两相的流动方程为：随油井含水率和井底气油比随油井含水率和井底气油比上升上升, ,油的相对流动能力下油的相对流动能力下降降; ; 随气油比上升随气油比上升, ,液相流入能液相流入能力下降力下降; ;随含水率上升随含水率上升, ,液相液相流入能力增加。流入能力增加。第44页/共102页新型油井流入动态新型油井流入动态第一点是直线弯曲的始点第一点是直线弯曲的始点, ,该点处的流动压力等于饱和压力该点处的流动压力等于饱和压

34、力, ,流动压力低于该点以后流动压力低于该点以后, ,采油指数降低采油指数降低, ,产量增长速度减慢产量增长速度减慢; ;第二个特征点为最大产量点第二个特征点为最大产量点, ,该点对应的压力可称为油井最低允许流动压力该点对应的压力可称为油井最低允许流动压力, ,流动压力低于该点以后流动压力低于该点以后, ,产量开始降低产量开始降低主要原因：流动压力下降到一定程度以后主要原因：流动压力下降到一定程度以后, ,井底附近井底附近可能出现油气两相流动可能出现油气两相流动, ,使油相的流动能力急剧下降。使油相的流动能力急剧下降。生产压差对产量的贡献已经小于采油指数下降对产量生产压差对产量的贡献已经小于采

35、油指数下降对产量所产生的影响。所产生的影响。第45页/共102页新型油井流入动态新型油井流入动态油井最低允许流动压力为：油井最低允许流动压力为：第46页/共102页其它水平井产能计算模型：Borisov模型： srhLhLrBhkJwehooohh2/ln/4ln/543. 0Giger模型：srhLhrLrLBhkJwehehooohh2/ln/2/2/11ln/543. 02Joshi模型：shrhLhLLaaBhkJwooohh2/2/ln/2/2/ln/543. 022222Renard & Dupuy模型： srhLhXBhkJwooohh12/ln/cosh/543. 0第47页/

36、共102页dpBKrrkhqrwfPPooroweo43ln2（1）基本公式 当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般表达式为:三、 时的流入动态wfbrPPP第48页/共102页图1-11 组合型IPR曲线（2）实用计算方法(组合型IPR方法)第49页/共102页)(8 . 02 . 01 2bwfbwfcboPPPPqqq) 1(8 . 18 . 1brbbcPPqJPq)(brbPPJq流压等于饱和压力时的产量为：)(8 . 0)(2 . 01 8 . 12bwfbwfbbroPPPPPPPqJ 当 时，由于油藏中全部为单相

37、液体流动。bwfPP 当 后，油藏中出现两相流动。bwfPP)(wfroPPJq流入动态公式为：直线段采油指数第50页/共102页四、油气水三相IPR 曲线Petrobras提出了计算三相流动IPR曲线的方法。综合IPR曲线的实质: 按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算采液指数时，是按产量加权平均；当预测产量或流压时是按流压加权平均。 图1-12 油气水三相IPR 曲线第51页/共102页已知一个测试点( 、 )和饱和压力 及油藏压力wftestpttestqbPrP 当 时：bwftestPPwftestrttestlPPqJ思考题：推导 时的采液指数计算式。

38、bwftestPP 当 时：)(8 . 1)1 (wftestrwbbrwttestlPPfAPPPfqJbwftestwfGPPP 当 时：wfGwftestPP0bwrwbwftestwttestlpfpfppfqJ49849828 . 02 . 01bwftestbwftestPPPPA其中：bwwfGpfP94直线段采油指数第52页/共102页 ，则：btqq 0ltrwfJqPP ，则按流压加权平均进行推导：maxotbqqqwfwaterwwfoilwwfPfPfP)1 (bobtbwfoilqqqqPPmax80811125. 0ltrwfwaterJqPp (二)某一产量下的流

39、压计算)(80811)1 (125. 0maxbobtbwltrwwfqqqqPfJqPfP所以：第53页/共102页lwqqtwfJfdqdPot98max因为：lwotlorwwfJfqqJqPfP)98)()(maxmax所以：若 ，则综合IPR曲线的斜率可近似为常数。maxmaxttoqqq思考题：试推导lwqqtwfJfdqdPot98max第54页/共102页五、多层油藏油井流入动态（1）多油层油井流入动态迭加型IPR图1-13 多层油藏油井流入动态第55页/共102页（2）含水油井流入动态图1-14 含水油井流入动态与含水变化 ( )soswPP 图1-15 含水油井流入动态曲线

40、 ( )swsoPP 第56页/共102页小 结(1) 上述介绍的方法阐明了油井流入动态的物理意义，也是目前现场最常用的计算方法。(2) 油井流入动态研究主要有三种途径： 基于Vogel、Fetkovich、Petrobras方法的完善。 建立不同类型油藏和井底条件的渗流模型。 利用单井流入动态的油藏数值模拟技术。(3) 油井流入动态是采油工程各项技术措施设计、分析与评价的依据。第57页/共102页第二节 井筒气液两相流基本概念井筒多相流理论井筒多相流理论：研究各种举升方式油井生产规律基本理论研究特点研究特点：流动复杂性、无严格数学解研究途径研究途径：基本流动方程 实验资料相关因次分析 近似关

41、系第58页/共102页(一)气液两相流动与单相液流的比较比较项目单相液流气液两相流能量来源井底流压井底流压气体膨胀能能量损失重力损失摩擦损失重力损失摩擦损失动能损失流动型态基本不变流型变化能量关系简单复杂第59页/共102页流动型态（流动结构、流型）：流动型态（流动结构、流型）： 流动过程中油、气的分布状态。流动过程中油、气的分布状态。影响流型的因素影响流型的因素各相介质的体积比例各相介质的体积比例介质的流速介质的流速各相的物理及化学性质各相的物理及化学性质(密度、粘度界面张力等密度、粘度界面张力等)流道的几何形状流道的几何形状壁面特性壁面特性管道的安装方式管道的安装方式( (二二) )气液混

42、合物在垂直管中的流动结构变化气液混合物在垂直管中的流动结构变化第60页/共102页流动型态的划分方法：两类流动型态的划分方法：两类第一类划分方法：第一类划分方法：根据两相介质分布的外形划分，根据两相介质分布的外形划分，包括泡包括泡状流、弹状流或团状流、（层状流、波状流）、段塞流或状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流冲击流、环状流、雾状流 垂直气液两相流流型垂直气液两相流流型 水平气液两相流流型水平气液两相流流型第61页/共102页第二类划分方法：第二类划分方法：按流动的数学模型或流体的分散程度划分，按流动的数学模型或流体的分散程度划分，包括包括分散流、间歇流

43、、分离流分散流、间歇流、分离流 两种分类方法比较两种分类方法比较 第一类划分方法较为直观第一类划分方法较为直观第二类划分方法便于进行数学处理第二类划分方法便于进行数学处理第一类划分方法第一类划分方法泡状流泡状流弹状流或团状流弹状流或团状流层状流层状流波状流波状流段塞流或冲击流段塞流或冲击流环状流环状流雾状流雾状流第二类划分方法第二类划分方法分散流分散流间歇流间歇流分离流分离流分离流分离流间歇流间歇流分离流分离流分散流分散流两类划分结果两类划分结果的对应关系的对应关系第62页/共102页, ( (二二) )气液混合物在垂直管中的流动结构变化气液混合物在垂直管中的流动结构变化泡流泡流 当当井筒压力

44、稍低于饱和压力井筒压力稍低于饱和压力时，溶解气开始从油中分时，溶解气开始从油中分离出来，气体都以小气泡分散在液相中。离出来，气体都以小气泡分散在液相中。滑脱现象滑脱现象： 混合流体流动过程中，由于流体间的密度差异，引混合流体流动过程中，由于流体间的密度差异，引起的小密度流体流速大于大密度流体流速的现象。如：起的小密度流体流速大于大密度流体流速的现象。如：油气滑脱、气液滑脱、油水滑脱等。油气滑脱、气液滑脱、油水滑脱等。特点特点：气体是分散相，液体是连续相；气体主要影响混合物密度，气体是分散相，液体是连续相；气体主要影响混合物密度，对摩擦阻力影响不大；滑脱现象比较严重。对摩擦阻力影响不大；滑脱现象

45、比较严重。 纯液流纯液流: :当当井筒压力大于饱和压力井筒压力大于饱和压力时，天然气溶解在原油中，产液时，天然气溶解在原油中，产液呈单相液流。呈单相液流。第63页/共102页段塞流 当混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡将合并成大气泡，直到能够占据整个油管断面时，井筒内将形成一段液一段气的结构。特点：气体呈分散相，液体呈连续相； 一段气一段液交替出现； 气体膨胀能得到较好的利用； 滑脱损失变小； 摩擦损失变大。第64页/共102页环流 油管中心是连续的气流而管壁为油环的流动结构。特点：气液两相都是连续相； 气体举油作用主要是靠摩擦携带； 滑脱损失变小； 摩擦损失变大。第65页

46、/共102页 雾流 气体的体积流量增加到足够大时，油管中内流动的气流芯子将变得很粗，沿管壁流动的油环变得很薄，绝大部分油以小油滴分散在气流中。特点：气体是连续相，液体是分散相； 气体以很高的速度携带液滴喷出井口； 气、液之间的相对运动速度很小； 气相是整个流动的控制因素。第66页/共102页总结： 油井生产中可能出现的流型自下而上依次为：纯油(液)流、泡流、段塞流、环流和雾流。 实际上，在同一口井内，一般不会出现完整的流型变化。图1-17 油气沿井筒喷出时的流型变化示意图纯油流；泡流；段塞流；环流；雾流第67页/共102页mmmlmff实际计算：直接求存在滑脱混合物密度或包括滑脱在内的摩擦阻力

47、系数。(三)滑脱损失概念因滑脱而产生的附加压力损失称为滑脱损失。单位管长上滑脱损失为：图1-18 气液两相流流动断面简图滑脱损失的实质： 液相的流动断面增大引起混合物密度的增加。第68页/共102页量从断面流出的流体能的能量在断面和之间耗失体额外所做的功在断面和之间对流进入断面的流体能量二、井筒气液两相流能量平衡方程及压力分布计算步骤 两个流动断面间的能量平衡关系：(一)能量平衡方程推导222,VPU222,2mghmV121111,2,mghmVVPU图2-19 倾斜管流能量平衡关系示意图第69页/共102页222,VPU222,2mghmV121111,2,mghmVVPU图2-19 倾斜

48、管流能量平衡关系示意图2222221121112sin2sinVPmvmgZUqVPmvmgZU0)(sindqPVddZmgmvdvdU倾斜多相管流断面1和断面2的流体的能量平衡关系为：0sin1wdIdZgvdvdP第70页/共102页2sin2mmmmmmmvdfdZdvvgdhdP适合于各种管流的通用压力梯度方程：加速度摩擦举高)()()(dZdPdZdPdZdPdZdP则：令：2)()(sin)(2vdfdZdIdZdPdZdvvdZdPgdZdPw摩擦加速度举高第71页/共102页以计算段下端压力为起点，重复步，计算下一段的深度和压力，直到各段的累加深度等于管长为止。(2)多相垂直

49、管流压力分布计算步骤重复的计算，直至 。估计计算hh1)按深度增量迭代的步骤已知任一点(井口或井底)的压力作为起点，任选一个合适的压力降作为计算的压力间隔p。估计一个对应的深度增量h 。计算该管段的平均温度及平均压力，并确定流体性质参数。判断流型，并计算该段的压力梯度dp/dh。计算对应于p的该段管长(深度差)h。计算该段下端对应的深度及压力。第72页/共102页2)按压力增量迭代的步骤（略）思考题：根据上述步骤整理出计算压力分布的程序流程框图。说明：a. 计算压力分布过程中，温度和压力是相关的；b. 流体物性参数计算至关重要，但目前方法精度差；c. 不同的多相流计算方法差别较大，因此在实际应

50、用中有必要根据油井的实际情况筛选精度相对高的方法。第73页/共102页第三节 Orkiszewski方法 综合了Griffith & Wallis 和 Duns & Ros 等方法 处理过渡性流型时，采用Ros方法(内插法) 针对每种流动型态提出存容比及摩擦损失的计算方法 提出了四种流型，即泡流、段塞流、过渡流及环雾流 把Griffith段塞流相关式改进后推广到了高流速区 1967年提出，适用于垂直管流计算第74页/共102页mmmmfdvvdhgdhdPdpPAqWdvvpgtmmm2出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：一、压力降公式及流动型态划分界限

51、由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和：kPkpgtfmhPAqWg12所以压降计算式为：第75页/共102页表1-3 Orkiszewski方法流型划分界限流 动 型 态界 限泡 流BtgLqq段 塞 流SgBtgLvLqq,过 渡 流SgMLvL雾 流MgLv不同流动型态下 和 的计算方法不同。mf第76页/共102页二、平均密度及摩擦损失梯度的计算气相存容比气相存容比( (含气率含气率) )Hg ：管段中气相体积与管段容 积之比值。液相存容比液相存容比( (持液率持液率) )HL ：管段中液相体积与管段容 积之比值。(1)泡流1gLHHggLgggLLmHH

52、HH)1 (平均密度:第77页/共102页)1 (1gpgtgpggsLgsgsHAqqHAqHvHvv22LHLtvDf)1 (gpLLHHAqv滑脱速度：气相流速与液相流速之差。4)1 (1 212psgpstpstgAvqAvqAvqH则：泡流摩擦损失梯度按液相进行计算：第78页/共102页 图1-21 摩擦阻力系数曲线第79页/共102页(2)段塞流lpstpsltmAvqAvW平均密度:段塞流的摩擦梯度：)(22pstpsltlfAvqAvqDvf段塞流计算中，关键是滑脱速度vs的计算。目前， vs的计算方法有两种：查图迭代法和经验公式法。第80页/共102页llsbDvN泡流雷诺数

53、：图1-22 C1Nb曲线lltDvNRe雷诺数：图1-23 C2NRe曲线滑脱速度的计算迭代法gDCCvs21滑脱速度：第81页/共102页滑脱速度的计算经验公式计算法详见教材p3839公式（1-54）（1-58）。MiSMsgSLSMgMmLLLvLLvL(3)过渡流 过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。MiSMggSLSMgMtLLLvLLvL第82页/共102页ggLgggLLmHHHH)1 (雾流混合物平均密度计算公式与泡流相同：由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。雾流摩擦系数可根据气体雷诺数和

54、液膜相对粗糙度查图得。Dvfsggf22摩擦梯度:(4)雾流gLggqqqH所以：第83页/共102页以井口油压或井底流压为起点，选择合适的压力间隔P，假设h计算平均 P 和 T，并求得在此 P 和 T 下的流体性质参数和流动参数，以及相应的流动型态界限 LB、Lg和 LM确定流动型态雾流计算气相存容比、平均密度及摩擦梯度过渡流分别按段塞流和雾流计算平均密度及摩擦梯度，并进行内插段塞流计算滑脱速度、液体分布系数、平均密度和摩擦梯度泡流计算气相存容比、平均密度和摩擦梯度计算并比较h，重复上述计算使h的计算值与假设值相等或在允许的误差范围内重复上述步骤，直到h的等于或大于油层深度为止图1-24 O

55、rkiszewski方法计算流程框图第84页/共102页第四节 Beggs & Brill 方法 水和空气、聚丙烯管实验基础上总结的方法 建立流型分布图，将七种流型归为三类，增加了过渡流 计算时先按水平管流计算，然后采用倾斜校正系数校正成相应的倾斜管流 倾斜度 -90 +90，分上坡和下坡流动 1973年提出，适用于水平、垂直和任意倾斜管流计算第85页/共102页Beggs & Brill 两相水平管流型分离流分层流波状流环状流间歇流团状流段塞流分散流泡 流雾 流第86页/共102页加速度摩擦位差dZdPdZdPdZdPdZdPDAGDvdZdp2/22摩擦一、基本方程 单位质量气液混合物稳定

56、流动的机械能量守恒方程为：(1)位差压力梯度：消耗于混合物静水压头的压力梯度。(2)摩擦压力梯度：克服管壁流动阻力消耗的压力梯度。假设条件:气液混合物既未对外作功，也未受外界功。sin1singHHgdZdpLgLL位差第87页/共102页dZdvvdZdp加速度dZdpPvvdZdpsg加速度AQvgsg/PvvHHDAGgHHdZdPsgLgLLLgLL/)1 (12sin)1 (3)加速度压力梯度：由于动能变化而消耗的压力梯度。忽略液体压缩性、考虑到气体质量流速变化远远小于气体密度变化，则：(4)总压力梯度（Beggs-Brill方法的基本方程）第88页/共102页 图1-26 Beggs-Brill流型分布图分离流； 间歇流； 分散流； 过渡流二、 Beggs & Brill方法的流