第1章油井流入动态及多相流动计算(4.7教室)

第一章油井流入动态与井筒多相流动计算★油井流入动态★井筒气液两相流基本概念★计算气液两相垂直管流的Orkiszewski方法★计算井筒多相管流的Beggs-Brill方法油井生产过程井口到分离器的流动－近似水平管流油气从油藏流到井底－地层中的渗流从井底流到井口－多相管流四个基本流动过程：通过油嘴的流动－嘴流第一节油井流入动态（IPR曲线）wfPoQ油气井流入动态：在一定的油层压力下，流体（油，气，水）产量与相应的井底流压的关系，反映了油藏向该井供油气的能力。达西定律线性流径向流2lnooeoowkhPqrBroookAPqLIPR曲线表示产量与井底流压关系的曲线(InflowPerformanceRelationshipCurve),称为流入动态曲线，简称IPR曲线。复合型：PwfPbPe，单--两相渗流直线型：PwfPb,单相渗流，牛顿流体（刚性水驱）曲线型：a.非牛顿流体单相驱，b.PwfPb,两相流，溶解气驱，粘弹流体。IPR曲线的影响因素：油藏驱动类型；完井状况；油藏及流体物性。反映了油藏向井的供油能力；反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响；通过油井流入动态研究为油藏工程提供检验资料；为采油工程的下一步工作提供依据；检查钻井、固井、完井和各项工艺措施等技术水平的优劣。IPR曲线的主要作用：sXBhakJooo21ln2)(wfroPPqJ一、单相流体流入动态asrrBPPhkqweoowfroo21ln)(2供给边缘压力不变圆形地层中心一口井的产量公式为：asrrBPPhkqweoowfroo43ln)(2圆形封闭油藏，拟稳态条件下的油井产量公式为：1.符合线性渗流规律时的流入动态（一）单相液体流入动态对于非圆形封闭泄油面积的油井产量公式，可根据泄油面积和油井位置进行校正。Xrrwe泄油面积形状与油井的位置系数图2()1ln2orwfoookhPPqaBXs2()3ln4orwfoookhPPqaBXs油井产量公式变为：sXBhakJooo43ln2sXBhakJooo21ln2单相流动时，油层物性及流体性质基本不随压力变化:)(wfroPPJq)(wfroPPqJJm1采油(液)指数：单位生产压差下的油井产油(液)量，M3/(d.Mpa)。►产油量与生产压差之比；►每增加单位生产压差时，油井产量的增加值；►油井IPR曲线斜率的负倒数。其它定义形式：采油指数J的确定：◆根据试井资料绘制asrrBPPhkqweoowfroo21ln)(2)(wfroPPJq◆根据油藏参数计算：改变油井工作制度，当油井稳定生产后，测得3～5个稳定工作制度下的产量及其流压，便可绘制IPR曲线。再根据IPR曲线求取。注意问题：对于非直线型IPR曲线，由于斜率不是定值，按不同定义求得的采油指数不同。在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。例：A井位于正方形泄油面积的中心，,3,2.1,1.0,10424sBrmAow根据系统试井，计算（流动系数）ooRohkPJ,,采油指数pwf11.210.29.79.1q16.135.946.157.9。MPapwf时的产量及8.8例题1）绘制IPR曲线解：pwf11.210.29.79.1q16.135.946.157.9平均地层压力MpamppqqJwfwfooo321122091120602）计算采油指数3）查表得mPasmmpasmsrrBJhkrrCweooowex2394107.0))/((104107.02)43(ln11421.040000571.0,571.04）直线外推至q=0,求.12MpaPR2.符合非线性渗流规律时的流入动态条件：油井产量很高时，在井底附近不再符合线性渗流，呈现高速非线性渗流。2DqCqPPwfrhakSxBCooo2)43(lnworhBD222134103396.1201.1710906.1k55.061008.1kg胶结地层的紊流速度系数：非胶结地层紊流速度系数：如果试井资料在单相渗流呈现非线性渗流，可绘制与的关系曲线。()RwfppqqDqCqPPwfr由此可以看出，与呈线性关系。绘制的直线的斜率为D，其截距则为C。qPPwfr/q（二）单相气体流入动态1.符合线性渗流规律时的流入动态drdprhkqggr2条件：定压边界、圆形气层中心有一口气井稳定生产时，距井轴r处的流量为：根据气体连续方程和状态方程，将半径r处的流量折算为标准状态下的气井产量qg。ppwscscscggwfdpZprrTpZhTKq2ln引用假(拟)压力的概念:02pppdpZ)(lnwfwscscscggrrTpZhTKqwfwfoppppppdpZpdpZpdpZp222wfreerr因为所以wfe用数值积分法或其它方法求得拟压力后，再求得气井产量。在工程中常近似地用平均压力求μ和Z2)(wfepppdpprrZTpZhTKqewfppwegscscscgg2ln即积分得wegscwfescscggrrZTpppZhTKqln)(22wegscscscgrrZTpZhTKDln令则)(22wfegppDq第八章1)绘制与的关系曲线；)(22wfeppgq应用：2)求出斜率D；3)预测稳定线性渗流条件下的气井流入动态。22()/ewfgDppq)(22wfegppDq2.符合非线性渗流规律时的流入动态(1)二项式方程矿场上将与的关系曲线称为二项式特征曲线。gqp2gq222ggwfeBqAqppggwfeBqAqpp2222ggBqAqpBpBAAqg2422利用上式求得一组不同下的，即可绘出气井的流入动态曲线。wfpgq将==1.03×105Pa时的产量称为气井的绝对无阻流量wfpapAOFqBppBAAwfe2)(4222用途：衡量气井生产能力的大小；进行气井间生产能力的比较；是气井配产的依据。nwfegppCq)(22C——与气层及流体性质有关n——渗流指数，0.5≤＜1。)lg(lglg22wfegppnCqnaeAOFppCq)(22利用试井资料求出C，n值，求得其绝对无阻流量与液体相比，气体具有更大的压缩性，应注意问题：气体的产量均指标准状态下的产量气井的压力均采用绝对压力而非表压力（2）指数式方程二、油气两相渗流时的流入动态o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制IPR曲线。平面径向流，直井油气两相渗流时油井产量公式为：(一)垂直井油气两相渗流时的流入动态KKKroo1.Vogel方法①假设条件：a.圆形封闭油藏，油井位于中心；b.均质油层，含水饱和度恒定；c.忽略重力影响；d.忽略岩石和水的压缩性；e.油、气组成及平衡不变；f.油、气两相的压力相同；g.拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。•1968年，沃格尔对不同流体性质、油气比、相对渗透率、井距、压裂井、污染井等各种情况下的21个溶解气驱油藏进行了计算。低高排除高粘度原油及严重污染的油井后，绘制了一条参考曲线，这一曲线被称为沃格尔曲线。Vogel曲线及方程]8.02.01[2maxrtestwfrtestwftestooPPPPqqmaxoqa.计算max28.02.01orwfrwfoqPPPPqc.根据给定的流压及计算的相应产量绘制IPR曲线。b.给定不同流压，计算相应的产量：已知地层压力和一个测试点：利用Vogel方程绘制IPR曲线的步骤AACBBPr242121qqA12212.0wfwfPPqqB2122218.0wfwfPPqqC◆油藏压力未知，已知两个工作点a.确定油藏平均压力maxoqb.计算d.根据给定的流压及计算的相应产量绘制IPR曲线c.给定不同流压，计算相应的产量dpBksrrkhqrwfPPooroweo43ln22.费特柯维奇方法溶解气驱油藏：假设与压力成直线关系，则：ooroBkprwfPPweocpdpsrrkhq43ln222243ln2wfrwePPcsrrkhrpoororKPc)(1rwfrPooroweoPPPKsrrkhqr243ln222rpooroweoPKsrrkhJr2143ln2令：当时：0wfP243ln2maxrPooroweoPKsrrkhqr所以：2max1rwfooPPqq)(22wfroPPJ3.非完善井Vogel方程的修正油水井的非完善性：◆打开性质不完善；如射孔完成◆打开程度不完善；如未全部钻穿油层◆打开程度和打开性质双重不完善◆油层受到损害◆酸化、压裂等措施完善井和非完善井周围的压力分布示意图非完善井:wssseooowfeorrkrrkBPPhqln1ln1)(2shkBqPoooosk2非完善井附加压力降:1ln2oooosskoswqBkrPkhkr则：wssorrkksln1FEwfrwfrpppp=wfrskwfrpppppwfwfskPPPskP为“正”称“正”表皮，油井不完善；skP为“负”称“负”表皮，油井超完善。油井的流动效率（FE）：油井的理想生产压差与实际生产压差之比。0s1FE油层受污染的或不完善井，0s1FE完善井,0s1FE增产措施后的超完善井，shkBqPoooosk2利用流动效率计算非完善直井流入动态无因次IPR曲线①Standing方法(FE=0.5~1.5)2)1max(8.02.01rwfrwfFEooPPPPqqStanding方法不完善井IPR曲线的步骤：2)1max(8.02.01rwfrwfoFEoPPPPqqFEPPPPwfrrwf)(a.根据已知Pr和Pwf计算在FE=1时最大产量21max8.02.01rwfrwfFEooPPPPqqb.预测不同流压下的产量c.根据计算结果绘制IPR曲线)5.2~1(FE②Harrison方法Harrison无因次IPR曲线(FE1)1maxmax1maxmaxFEoFEoFEoFEoqqqqb.求FE对应的最大产量，即Pwf=0时的产量Harrison方法计算超完善井IPR曲线步骤：a.计算FE=1时的qomax(FE=1)先求Pwf/Pr，然后查图1-7中对应的FE曲线上的相应值qo/qomax(FE=1)。则)1max()1max(/FEoooFEoqqqq1max1maxFEooFEooqqqqc.计算不同流压下的产量d.根据计算结果绘制IPR曲线2PCPBAq(二)斜井和水平井的IPR曲线Cheng对溶解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似Vogel方程的不同井斜角井的IPR回归方程：P′=Pwf/P