CN105257279A - 一种抽油机井动液面的测量方法 - Google Patents
一种抽油机井动液面的测量方法 Download PDFInfo
- Publication number
- CN105257279A CN105257279A CN201510702552.1A CN201510702552A CN105257279A CN 105257279 A CN105257279 A CN 105257279A CN 201510702552 A CN201510702552 A CN 201510702552A CN 105257279 A CN105257279 A CN 105257279A
- Authority
- CN
- China
- Prior art keywords
- parameter
- pump
- load
- rod string
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000005086 pumping Methods 0.000 title abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims abstract description 59
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000000704 physical effect Effects 0.000 claims abstract description 11
- 238000004519 manufacturing process Methods 0.000 claims abstract description 9
- 238000009826 distribution Methods 0.000 claims description 28
- 239000007788 liquid Substances 0.000 claims description 27
- 230000005484 gravity Effects 0.000 claims description 10
- 238000013210 evaluation model Methods 0.000 claims description 4
- 239000003129 oil well Substances 0.000 abstract description 21
- 238000010586 diagram Methods 0.000 abstract description 4
- 239000003921 oil Substances 0.000 description 36
- 239000007789 gas Substances 0.000 description 14
- 230000003068 static effect Effects 0.000 description 9
- 239000010779 crude oil Substances 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- OVSKIKFHRZPJSS-UHFFFAOYSA-N 2,4-D Chemical compound OC(=O)COC1=CC=C(Cl)C=C1Cl OVSKIKFHRZPJSS-UHFFFAOYSA-N 0.000 description 2
- 229910000975 Carbon steel Inorganic materials 0.000 description 2
- 241001442234 Cosa Species 0.000 description 2
- 244000089409 Erythrina poeppigiana Species 0.000 description 2
- 235000009776 Rathbunia alamosensis Nutrition 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000010962 carbon steel Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Landscapes
- Reciprocating Pumps (AREA)
Abstract
本发明公开了一种抽油机井动液面的测量方法,根据阻力系数计算模型,获取阻力系数;根据所述阻力系数、抽油杆柱位移运动模型和抽油杆柱载荷计算模型,获取所述抽油杆柱在不同深度位置处的位移参数及其载荷参数;从而获取泵的位移参数及其载荷参数;根据所述泵的位移参数及其载荷参数,获取所述泵的泵功图;根据所述泵功图和沉没压力计算模型以及流体物性计算模型,获取所述抽油机井的动液面参数。本发明提供的抽油机井动液面的测量方法,在确保计算精度的情况下,能够有效提高其适时性,而且还能应用于控制套压生产的油井,使得计算出的动液面的误差降低。
Description
技术领域
本发明涉及石油技术领域,尤其涉及一种抽油机井动液面的测量方法。
背景技术
抽油机井动液面是了解油井的供液情况、诊断油井故障,进行采油工艺适应性评价和优化的重要参数。目前油井动液面的主要测试方法为声波反射法,此方法存在的主要问题有:⑴依靠工人手工操作;⑵井口发声装置存在一定的不安全因素;⑶每月测试一次,不能在线跟踪动液面的变化;⑷对于动液面深度大于1000米的油井,容易受到一些其它工况的影响,测试误差比较大。
近年来,由于示功图测试手段的发展,使得根据示功图计算动液面技术从理论向应用跨越,动液面的变化趋势跟踪技术逐步建立并不断完善,但是,在有通过测绘实测静载示功图计算油井动液面时,可以有效消除惯性载荷、振动载荷、摩擦载荷等因素的影响,计算结果精度相对较高,但同时,上下静载的获取比较复杂和繁琐,且校核困难,一旦油井工况发生改变,也需要对其静载示功图进行修正,适时性较差,而且采用该方法计算油井动液面时,没有考虑套压的影响,对于控制套压生产的油井,计算的动液面将会出现误差较高的问题。
发明内容
本发明提供了一种抽油机井动液面的测量方法,在确保计算精度的情况下,能够有效提高其适时性,而且还能应用于控制套压生产的油井,使得计算出的动液面的误差降低。
本发明实施例提供了一种抽油机井动液面的测量方法,包括:
根据阻力系数计算模型,获取与所述抽油机井对应的阻力系数;
根据所述阻力系数和抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数;以及
根据所述阻力系数和抽油杆柱载荷计算模型,获取所述抽油杆柱的载荷参数;
根据所述抽油杆柱的位移参数,获取所述抽油杆柱下部的泵的位移参数;以及根据所述抽油杆柱的载荷参数,获取所述泵的载荷参数;
根据所述泵的位移参数和所述泵的载荷参数,获取所述泵的泵功图;
根据所述泵功图,获取所述泵的吸入口沉没压力参数;
根据所述抽油机井的井筒的温度分布和流体物性计算模型,获取所述井筒的流体物性参数;
根据所述吸入口沉没压力参数和所述流体物性参数,获取所述抽油井的油管与套管之间的环形空间内的压力分布参数;
根据所述压力分布参数,获取所述抽油机井的动液面参数。
可选的,所述阻力系数计算模型,具体为:
其中,
式中μ表示液体粘度;ρr表示抽油杆的密度;Ar表示抽油杆的截面积;Dt表示油管直径;Dr表示抽油杆直径;L表示抽油杆长度;C表示阻力系数。
可选的,所述根据所述阻力系数和抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数,具体为:
将所述阻力系数代入所述抽油杆柱位移运动模型中,获取所述抽油杆柱的位移参数,其中,所述抽油杆柱位移运动模型具体为:
式中:
U(x,t)表示t时刻抽油杆柱x断面处的位移;
a表示声波在抽油杆柱中的传播速度。
可选的,所述根据所述阻力系数和抽油杆柱载荷计算模型,获取所述抽油杆柱的载荷参数,具体包括:
将所述阻力系数代入所述抽油杆柱载荷计算模型中,获取所述抽油杆柱的载荷参数,其中,所述抽油杆柱载荷计算模型包括抽油杆重力计算模型、液柱载荷计算模型、流体通过凡尔孔的阻力计算模型、抽油杆柱惯性载荷计算模型,液柱惯性载荷计算模型、泵筒与柱塞摩擦载荷计算模型、抽油杆与液体的摩擦载荷计算模型、管液摩擦载荷计算模型和抽油杆柱上托力计算模型。
可选的,所述根据所述泵功图,获取所述泵的吸入口沉没压力参数,具体包括:
选取所述泵功图中的上下载荷点;
根据所述上下载荷点和泵吸入口沉没压力计算模型,获取所述吸入口沉没压力参数。
可选的,所述泵吸入口沉没压力计算模型具体为:
其中,ΔP表示油层产出液经过固定阀和游动阀孔时产生的压力降;Wl表示油井液柱载荷。
通过一个实施例或多个实施例,本发明具有以下有益效果或者优点:
由于本申请实施例中的抽油机井动液面的测量方法,根据阻力系数计算模型、抽油杆柱位移运动模型和抽油杆柱载荷计算模型,获取所述抽油杆柱的位移参数和所述抽油杆柱的载荷参数;根据所述抽油杆柱的位移参数,获取所述抽油杆柱下部的泵的位移参数;以及根据所述抽油杆柱的载荷参数,获取所述泵的载荷参数;根据所述泵的位移参数和所述泵的载荷参数,获取所述泵的泵功图;根据所述泵功图,获取所述泵的吸入口沉没压力参数;根据所述抽油机井的井筒的温度分布和流体物性计算模型,获取所述井筒的流体物性参数;根据所述吸入口沉没压力参数和所述流体物性参数,获取所述抽油井的油管与套管之间的环形空间内的压力分布参数;根据所述压力分布参数,获取所述抽油机井的动液面参数,如此,可以通过阻力系数计算模型、抽油杆柱位移运动模型和抽油杆柱载荷计算模型来消除惯性、振动及摩擦载荷的影响,并求解得到泵功图的静载荷,使得获取的泵功图的静载荷的精确度更高,而使得再根据泵功图的静载荷进行动液面计算获得的动液面参数的精确度也得以提高,而且通过上述计算模型能够实时获取到动液面面参数,如此,能够有效提高其适时性,而且本申请提供的测量方法是根据所述环形空间内的压力分布参数来获取当动液面参数,如此,使得本申请提供的测量方法还能应用于控制套压生产的油井,使得计算出的动液面的误差降低。
附图说明
图1为本发明实施例中泵的工作原理图的;
图2为本发明实施例中井筒流体分布图。
图中有关附图标记如下:
10——泵简,11——柱塞,12——游动阀,13——固定阀,14——抽油杆,20——气柱段,21——含气油柱段,22——混合液段。
具体实施方式
本发明提供了一种抽油机井动液面的测量方法,在确保计算精度的情况下,能够有效提高其适时性,而且还能应用于控制套压生产的油井,使得计算出的动液面的误差降低。
为了更好的理解上述技术方案,下面将结合说明书附图以及具体的实施方式对上述技术方案进行详细的说明。
在叙述本发明一种抽油机井动液面的测量方法之前,首先需要了解泵的工作原理,下面具体以有杆泵为例,具体的,参见图1,泵包括泵简10、柱塞11、游动阀12、固定阀13和抽油杆14,在泵工作过程中,泵筒10内压力P(t)随柱塞运动方向的改变,由吸入压力Pi升至排出压力PP或由PP降至Pi,柱塞11完成卸载或加载:当固定阀13开启后,液体经固定阀13的孔吸入泵腔,此时P(t)=Pi,柱塞11加载完成,泵载保持不变;当游动阀12开启后,液体经游动阀12的孔排出泵腔,此时P(t)=PP,柱塞11卸载完成,泵载保持不变。
如果忽略柱塞与液体的惯性力,则作用于柱塞上的平衡方程应是:
上冲程:
Fp(t)=Pp·(fp-fr)-P(t)·fp+Wp+f公式1
下冲程:
Fp(t)=Pp·fr-P(t)·fp+Wp-f公式2
在公式1和公式2中,Fp(t)表示泵的载荷,单位:N;Pp表示游动阀12上部的压力,单位:Pa;P(t)表示泵筒10内压力,单位:Pa;Wp表示柱塞11重量,单位:N;f表示柱塞11与泵筒10间的摩擦阻力,单位:N;fp表示柱塞11的截面积,单位:m2;fr表示抽油杆14的截面积,单位:m2。
其中,泵沉没度对应的泵吸入口沉没压力参数与上冲程时泵的吸入压力之间的存在确定关系,因此可由示功图求出所述泵吸入口沉没压力参数,然后利用井筒多相流模型计算出所述抽油井的油管与套管之间的环形空间内的压力分布参数,所述泵吸入口沉没压力参数与所述压力分布参数的交点即为动液面的深度,如此,本申请的测量方法是用于获取所述泵吸入口沉没压力参数与所述压力分布参数,以此来获取动液面面参数,所述动液面面参数包括动液面面的深度。
其中,本发明提供的一种抽油机井动液面的测量方法,包括以下步骤:
步骤100:根据阻力系数计算模型,获取与所述抽油机井对应的阻力系数。
在具体实施过程中,泵的实际冲程可以通过泵功图获得,然而准确计算泵功图之前必须先知道阻力系数。在这种情况下,可以通过吉布斯阻力系数计算模型来获取所述阻力系数,其中,所述吉布斯阻力系数计算模型具体为:
在公式3中,μ表示液体粘度,单位:Pa·s;ρr表示抽油杆的密度,单位:Kg/m3;Ar表示抽油杆的截面积,单位:m2;Dt表示油管直径,单位:m;Dr表示抽油杆直径,单位:m;L表示抽油杆长度,单位:m;a表示声波在抽油杆柱中的传播速度(声波在碳素钢中传播速度大致为4900m/s);ω表示曲柄角速度;m表示油管内径与抽油杆直径之比。
具体的,公式3只是单级抽油杆柱计算模型,只适用于单级抽油杆杆柱,多级抽油杆杆柱可将求解分为多级计算,每一级的尾端数据作为下一级计算的边界条件进行多级求解即可得到。
步骤101:根据所述阻力系数和抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数。
在具体实施过程中,抽油杆柱位移运动模型可用带阻尼的波动方程来描述,具体如下:
在公式4中,U(x,t)表示t时刻抽油杆柱x断面处的位移;a表示声波在抽油杆柱中的传播速度(声波在碳素钢中传播速度大致为4900m/s);c表示阻尼系数.如此,在通过步骤100获取到c的值之后,将所述c的值带入所述抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数具体可以用U(x,t)表示。
其中,公式4的边界条件为悬点动负荷函数D(t)和光杆位移函数U(t),两个方程都用截尾傅里叶级数的方式由公式给出给出:
因为方程4中没有考虑抽油杆柱重力,所以实际动载荷函数D(t)应为悬点载荷减去抽油杆重力,并用动负荷函数作为模型的边界条件。D(t)及U(t)的傅里叶系数σ0、σn、τn及ν0、νn、δn求解方程如下所示:
其中,上述公式中,ω表示曲柄角速度;T表示抽汲周期。
具体的,在实际工作中D(t)及U(t)是以示功图数据点的形式给出的,可以通过傅里叶系数可用近似的面积积分来确定,所述傅里叶级数n最大取值可以为6,积分求解采用分割方法将积分图像分割成144块等宽举行进行近似求解。用分离变量的方法对方程进行化简求解,得到t时刻抽油杆柱x断面处的位移变化关系如下式所示:
t时刻抽油杆柱x断面处的动负荷函数的变化关系如式所示:
在t时刻,抽油杆柱x断面上的总载荷=F(x,t)+抽油杆柱x断面以下的杆柱重量计算模型如下所示:
On(x)=(Knchβnx+δnshβnx)sinanx+(μnshβnx+νnchβnx)cosanx公式48
Pn(x)=(Knshβnx+δnchβnx)cosanx+(μnchβnx+νnshβnx)sinanx公式49
步骤102:根据所述阻力系数和抽油杆柱载荷计算模型,获取所述抽油杆柱的载荷参数。
在具体实施过程中,将所述阻力系数代入所述抽油杆柱载荷计算模型中,获取所述抽油杆柱的载荷参数,其中,所述抽油杆柱载荷计算模型包括抽油杆重力计算模型、液柱载荷计算模型、流体通过凡尔孔的阻力计算模型、抽油杆柱惯性载荷计算模型,液柱惯性载荷计算模型、泵筒与柱塞摩擦载荷计算模型、抽油杆与液体的摩擦载荷计算模型、管液摩擦载荷计算模型和抽油杆柱上托力计算模型。
具体的,所述抽油杆重力模型用于消除抽油杆柱重力的影响,在分段计算的悬点载荷中直接减去对应长度抽油杆柱的重力即可,计算模型如下:
Fr=qrgL公式11
公式11中,Fr表示抽油杆柱在空气中的重力,单位:N;L表示抽油杆柱长度,单位:m;qr表示抽油杆柱的每米质量,单位:kg/m。
具体的,所述液柱载荷计算模型具体为:
Fl=(Ap-Am)×Pout-Ap×Pin公式12
公式12中,Fl表示作用在柱塞上的液柱载荷,单位:N;Am表示最下一级抽油杆截面积,单位:m2;Ap表示抽油泵活塞截面积,m2;Pout表示泵排出口处压力,单位:Pa;Pin表示泵吸入口处压力,单位:Pa。
具体的,所述流体通过凡尔孔的阻力计算模型具体为
公式13中,Fv表示流体通过凡尔孔的阻力,单位:N;ρl表示流体密度,kg/m3;f0表示凡尔孔过流面积,单位:m2;Sp表示活塞有效冲程,单位:m;N表示冲数,单位:rpm;μ表示由实验确定的凡尔流量系数,由下式计算:
公式14和公式15中,dv0表示凡尔孔直径,单位:m;μl表示流体粘度,单位:Pa.S。
具体的,所述抽油杆柱惯性载荷计算模型具体为:
公式16
公式17
公式16和公式17中,Fri表示抽油杆柱惯性载荷,单位:N;S表示冲程,单位:m;r表示抽油机曲柄半径,单位:m;l表示抽油机连杆长度,单位:m。
具体的,所述液柱惯性载荷计算模型具体为:
公式18中,Fli表示液柱惯性载荷,单位:N;At表示油管截面积,单位:m2。
具体的,所述泵筒与柱塞摩擦载荷计算模型具体为:
公式19中,Fp表示泵筒与柱塞的摩擦载荷,单位:N;dp表示泵柱塞直径,单位:m;de表示柱塞与衬套的间隙,单位:m;
具体的,所述抽油杆与液体的摩擦载荷计算模型具体为:
公式20中,Frl表示杆液摩擦载荷,单位:N;L表示抽油杆长度,单位:m;m表示油管内径与抽油杆直径之比,dt表示油管内径,单位:m;dr表示抽油杆内径,单位:m。
具体的,所述管液摩擦载荷计算模型具体为
公式21中,Ftl表示管液摩擦载荷,单位:N。
具体的,所述抽油杆柱上托力计算模型具体为
Frxj=Pj×(Arj-Arj+1)(j=1,2,...m)公式22
公式22中,Frxj表示第j级抽油杆下端面处所受的流体上托力,单位:N;Pj表示第j级抽油杆下端面处所受的流体压力,单位:Pa;Arj表示第j级抽油杆截面积,单位:m2;Arj+1表示第j+1级抽油杆截面积,单位:m2;m表示抽油杆级数。
步骤103:根据所述抽油杆柱的位移参数,获取所述抽油杆柱下部的泵的位移参数。
在具体实施过程中,由于所述抽油杆柱的位移参数与所述泵的位移参数相对应,如此,在通过步骤101获取到所述抽油杆柱的位移参数之后,根据所述抽油杆柱的位移参数即可获取所述泵的位移参数,其中,所述抽油杆柱的位移参数与所述泵的位移参数相同。
步骤104:根据所述抽油杆柱的载荷参数,获取所述泵的载荷参数。
在具体实施过程中,由于所述抽油杆柱的载荷参数与所述泵的载荷参数相对应,如此,在通过步骤102获取到所述抽油杆柱的载荷参数之后,根据所述抽油杆柱的载荷参数即可获取所述泵的载荷参数,其中,所述抽油杆柱的载荷参数与所述泵的载荷参数相同。
步骤105:根据所述泵的位移参数和所述泵的载荷参数,获取所述泵的泵功图。
在具体实施过程中,在获取到所述泵的位移参数和所述泵的载荷参数之后,直接根据所述泵的位移参数和所述泵的载荷参数,即可获取所述泵的泵功图。
具体来讲,通过建立数学模型进行计算,消除了抽油杆动载,杆柱伸缩等的影响,得到形状规整、能准确反映泵实际工作状况的泵功图。由于利用泵功图数据进行计算可以有效消除抽油杆动载、杆柱伸缩和杆柱摩擦阻力等的影响,如此,能够有效降低计算误差,使得根据所述泵功图而计算得到的动液面参数也更精确。
步骤106:根据所述泵功图,获取所述泵的吸入口沉没压力参数。
在具体实施过程中,选取所述泵功图中的上下载荷点;根据所述上下载荷点和泵吸入口沉没压力计算模型,获取所述吸入口沉没压力参数。
具体来讲,选取所述泵功图中的上下载荷点,再根据所述泵吸入口沉没压力计算模型计算得到所述吸入口沉没压力参数,油层产出液经过固定阀和游动阀孔时均会产生一个压力降ΔP,其中,所述泵吸入口沉没压力计算模型具体为:
公式23中:ρl表示原油的密度,单位:kg/m3;ξ表示阀流量系数(与阀孔直径、液体粘度、液体流速三者之间存在明确的函数关系,根据所述函数关系即可计算出)fo表示阀孔截面积,单位:m2。
公式24中,NRe表示雷诺数;do表示凡尔孔径,单位:m;Vf表示液流速度,单位:m/s;ν表示液体运动粘度,单位:m2/s。
其中,游动阀上部的压力Pp计算模型如下所示:
Pp=Ph+ρlgL公式25
上冲程泵载荷计算模型根据公式3至公式22可知,具体为:
P(t)=Pn-ΔP公式26
Fpu=Pp·(fp-fr)-(Pn-ΔP)·fp+Wp+f公式27
游动阀下冲程打开后至关闭前泵载荷可由如下公式计算:
P(t)=Pn+ΔP公式28
Fpd=Pp·fr-(Pp+ΔP)·fp+Wp-f公式29
根据公式28和公式29,可得到所述吸入口沉没压力参数,具体如下:
忽略掉柱塞与泵筒间的摩擦阻力f,fpu-fpd即为油井液柱载荷Wl,可由示功图数据得到,所述吸入口沉没压力参数具体为:
根据油井示功图得到油井液柱载荷Wl,即可利用公式31计算出所述吸入口沉没压力参数。
步骤107:根据所述抽油机井的井筒的温度分布和流体物性计算模型,获取所述井筒的流体物性参数。
在具体实施过程中,由于地下温度压力会对井筒流体物性参数产生一定影响,所以在计算过程中必须考虑温度压力对流体物性参数的影响以及井筒温度梯度分布即为所述温度分布参数。
其中,原油的密度计算模型如下所示:
公式32中,ρo表示在平均压力及平均温度条件下的原油密度,单位:kg/m3;γo表示地面条件下的原油相对密度;γg表示地面条件下的气体相对密度;Rs表示在平均压力及平均温度条件下的溶解气油比,单位:m3/m3;Bo表示在平均压力及平均温度条件下的原油体积系数。
具体,所述原油API度的计算模型如下:
具体的,原油的体积系数具体为:
Bo=0.972+0.000147F1.175公式34
其中:
其中,溶解气油比:
公式36和公式37中,T表示平均温度,单位:℃;P表示平均压力,单位:Pa。
具体的,油水混合液体的密度具体为:
ρl=ρo(1-fw)+ρwfw公式38
公式38中,fw表示体积含水率,单位:%。
具体的,天然气的密度具体为:
公式39中,ρg表示在给定温度和压力条件下天然气的密度,单位:kg/m3。
具体的,所述井筒温度场计算模型如下:
由于地层温度随着深度的变化不断改变,而温度的变化会导致井筒流体物性参数的改变,所以在设计中必须计算不同位置处地层温度,根据温度确定地层流体物性参数,根据经验公式可以计算沿井筒的温度分布:
公式40、公式41、公式42和公式43中,Ql表示油井产液量,单位:t/d;fw表示质量含水率,单位:%;t0表示地表恒温层温度,单位:℃;tr表示油层温度,单位:℃;H表示油层中部深度,单位:m;L表示计算点深度,单位:m。
步骤108:根据所述吸入口沉没压力参数和所述流体物性参数,获取所述抽油井的油管与套管之间的环形空间内的压力分布参数。
在具体实施过程中,参见图2,油管与套管之间的环形空间内的流体由于重力分异,一般形成三段:气柱段20、含气油柱段21及油、气、水进行混合的混合液段22,所述吸入口沉没压力参数为气柱段20和含气油柱段21两段压力之和:
Pn=PD+ΔPo公式44
公式44中,Pn表示吸入口沉没压力,单位:Pa;PD表示动液面处的压力,单位:Pa;ΔPo表示动液面到泵处的油柱产生的压力,单位:Pa。
其中,环空气柱压力PD计算方式如下:
由于抽油机井环空过流断面相对较大,而气流量一般来说不会很大,因此可以忽略不计由于摩阻和动能引起的压力损失。考虑到气柱的密度随压力、温度的变化而变化,动液面处的压力由下式求出:
公式45中,x表示从井口算起的深度,单位:m;g表示重力加速度,单位:m/s2;T0、P0表示标准状况下的温度和压力,单位:K和MPa;ρg0表示标准状况下的气体密度,单位:kg/m3;Tav表示平均温度下的气体压缩因子;Zav表示平均温度和压力下的气体压缩因子;Pso表示井口套压(取其绝对值),单位:MPa;PD(x)表示x处的气柱压力(取其绝对值),单位:MPa。
进一步的,动液面到泵处的油柱产生的压力ΔPo计算方式如下:
具体的,动液面到泵处的油柱产生的压力ΔPo,可由多相流模型计算油套环空中油柱压力梯度分布得到。
步骤109:根据所述压力分布参数,获取所述抽油机井的动液面参数。
在具体实施过程中,可以通过绘制抽油机井的示功图;首先确定液柱载荷;然后计算泵功图,确定液载;接着计算所述入口沉没压力参数;再计算环空液柱的压力分布参数;紧接着计算环形空间内的压力分布参数,最后确定出环空液柱压力分布参数与环形空间内的压力分布参数的交点即为动液面参数,其中,本申请文件中的参数均可以使用曲线表示。
在实际应用过程中,通过本申请的测量方法获取的动液面参数与实际动液面深度的对比表具体如下表1所示。
表1
其中,通过建立抽油机井动液面预测模型,将地面功图转化为泵功图,再对泵功图进行定量处理,进而结合油井的基础数据和原油物性参数计算得到油井动液面。经过现场试验,预测精度3%左右,具有精度高、可靠性好的优点,而且能够实时进行预测,提高其预测的实时性。
在实际应用过程中,以某实际生产油井为例通过本申请的测量方法进行了油井动液面深度的测量,首先采集所述该油井的相关技术资料,包括油井正常生产的基础生产资料,抽油机光杆冲程、冲次以及抽油泵泵径、下泵深度、原油密度、含水率等资料。对于多级抽油杆柱组合***还应包括各级抽油杆柱的直径、杆长以及钢级信息,示功图数据为载荷位移数据,具体获取油井数据如表2和表3所示。
油层中深/m | 1921.95 | 泵径/mm | 44 |
气相对密度 | 0.7 | 泵深/m | 1696.3 |
原油相对密度 | 0.9 | 冲程/m | 6 |
地层温度/℃ | 60 | 冲次 | 1.8 |
含水率 | 0.633 | 产液量/t | 8.04 |
原油粘度/mPa.s | 640 | 产油量/t | 2.95 |
生产汽油比 | 10 | 油压/MPa | 0.86 |
油管内径/mm | 62.0 | 套压/MPa | 1.83 |
套管内径/mm | 139.7 | 回压/MPa | 0.83 |
表2
表3
具体的,根据表2和表3中的数据,通过本申请的测量方法对所述油井进行测量,获得所述油井的泵功图数据如表4所示。
表4
进一步的,从预测的泵功图即从表4中选取载荷最大值和载荷最小值,然后结合基础数据进行动液面预测,其中,表4中泵功图载荷最大值为24.1KN和最小值为-1.12KN,载荷差为25.22KN,计算得到泵吸入口压力为0.13MPa,再结合所述油井的示功图测量出所述油井的动液面参数为1668m,实际动液面深度为1680m,误差为-0.714%。
通过一个实施例或多个实施例,本发明具有以下有益效果或者优点:
由于本申请实施例中的抽油机井动液面的测量方法,根据阻力系数计算模型、抽油杆柱位移运动模型和抽油杆柱载荷计算模型,获取所述抽油杆柱的位移参数和所述抽油杆柱的载荷参数;根据所述抽油杆柱的位移参数,获取所述抽油杆柱下部的泵的位移参数;以及根据所述抽油杆柱的载荷参数,获取所述泵的载荷参数;根据所述泵的位移参数和所述泵的载荷参数,获取所述泵的泵功图;根据所述泵功图,获取所述泵的吸入口沉没压力参数;根据所述抽油机井的井筒的温度分布和流体物性计算模型,获取所述井筒的流体物性参数;根据所述吸入口沉没压力参数和所述流体物性参数,获取所述抽油井的油管与套管之间的环形空间内的压力分布参数;根据所述压力分布参数,获取所述抽油机井的动液面参数,如此,可以通过阻力系数计算模型、抽油杆柱位移运动模型和抽油杆柱载荷计算模型来消除惯性、振动及摩擦载荷的影响,并求解得到泵功图的静载荷,使得获取的泵功图的静载荷的精确度更高,而使得再根据泵功图的静载荷进行动液面计算获得的动液面参数的精确度也得以提高,而且通过上述计算模型能够实时获取到动液面面参数,如此,能够有效提高其适时性,而且本申请提供的测量方法是根据所述环形空间内的压力分布参数来获取当动液面参数,如此,使得本申请提供的测量方法还能应用于控制套压生产的油井,使得计算出的动液面的误差降低。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将使显而易见的,本文所定义的一般原理可以在不脱离发明的精神或范围的情况下,在其他实施例中实现。因此,本发明将不会被限制与本文所示的这些实施例,而是要符合与本文所公开的原理和新颖性特点相一致的最宽的范围。
Claims (6)
1.一种抽油机井动液面的测量方法,其特征在于,包括:
根据阻力系数计算模型,获取与所述抽油机井对应的阻力系数;
根据所述阻力系数和抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数;以及
根据所述阻力系数和抽油杆柱载荷计算模型,获取所述抽油杆柱的载荷参数;
根据所述抽油杆柱的位移参数,获取所述抽油杆柱下部的泵的位移参数;以及
根据所述抽油杆柱的载荷参数,获取所述泵的载荷参数;
根据所述泵的位移参数和所述泵的载荷参数,获取所述泵的泵功图;
根据所述泵功图,获取所述泵的吸入口沉没压力参数;
根据所述抽油机井的井筒的温度分布和流体物性计算模型,获取所述井筒的流体物性参数;
根据所述吸入口沉没压力参数和所述流体物性参数,获取所述抽油井的油管与套管之间的环形空间内的压力分布参数;
根据所述压力分布参数,获取所述抽油机井的动液面参数。
2.如权利要求1所述的方法,其特征在于,所述阻力系数计算模型,具体为:
其中,
式中μ表示液体粘度;ρr表示抽油杆的密度;Ar表示抽油杆的截面积;Dt表示油管直径;Dr表示抽油杆直径;L表示抽油杆长度;C表示阻力系数。
3.如权利要求2所述的方法,其特征在于,所述根据所述阻力系数和抽油杆柱位移运动模型,获取所述抽油杆柱的位移参数,具体为:
将所述阻力系数代入所述抽油杆柱位移运动模型中,获取所述抽油杆柱的位移参数,其中,所述抽油杆柱位移运动模型具体为:
式中:
U(x,t)表示t时刻抽油杆柱x断面处的位移;
a表示声波在抽油杆柱中的传播速度。
4.如权利要求3所述的方法,其特征在于,所述根据所述阻力系数和抽油杆柱载荷计算模型,获取所述抽油杆柱的载荷参数,具体包括:
将所述阻力系数代入所述抽油杆柱载荷计算模型中,获取所述抽油杆柱的载荷参数,其中,所述抽油杆柱载荷计算模型包括抽油杆重力计算模型、液柱载荷计算模型、流体通过凡尔孔的阻力计算模型、抽油杆柱惯性载荷计算模型,液柱惯性载荷计算模型、泵筒与柱塞摩擦载荷计算模型、抽油杆与液体的摩擦载荷计算模型、管液摩擦载荷计算模型和抽油杆柱上托力计算模型。
5.如权利要求4所述的方法,其特征在于,所述根据所述泵功图,获取所述泵的吸入口沉没压力参数,具体包括:
选取所述泵功图中的上下载荷点;
根据所述上下载荷点和泵吸入口沉没压力计算模型,获取所述吸入口沉没压力参数。
6.如权利要求5所述的方法,其特征在于,所述泵吸入口沉没压力计算模型具体为:
其中,ΔP表示油层产出液经过固定阀和游动阀孔时产生的压力降;Wl表示油井液柱载荷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510702552.1A CN105257279A (zh) | 2015-10-26 | 2015-10-26 | 一种抽油机井动液面的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510702552.1A CN105257279A (zh) | 2015-10-26 | 2015-10-26 | 一种抽油机井动液面的测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105257279A true CN105257279A (zh) | 2016-01-20 |
Family
ID=55097169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510702552.1A Pending CN105257279A (zh) | 2015-10-26 | 2015-10-26 | 一种抽油机井动液面的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105257279A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106351645A (zh) * | 2016-10-26 | 2017-01-25 | 北京安控科技股份有限公司 | 一种抽油机井动液面连续测量方法及装置 |
CN106837305A (zh) * | 2016-12-28 | 2017-06-13 | 中国石油天然气股份有限公司 | 确定抽油井井下液面深度的方法和装置 |
CN107237615A (zh) * | 2017-08-09 | 2017-10-10 | 大连西贝通达科技有限公司 | 一种基于光杆载荷的游梁式抽油机运行速度优化控制方法 |
CN106777709B (zh) * | 2016-12-21 | 2018-11-06 | 中国石油大学(华东) | 基于地面实测电参数的潜油电泵井动液面计算的方法 |
CN109779612A (zh) * | 2017-11-14 | 2019-05-21 | 中国石油化工股份有限公司 | 一种碳杆举升***井下泵功图测试装置及方法 |
CN109899056A (zh) * | 2017-12-07 | 2019-06-18 | 中国石油天然气股份有限公司 | 抽油机的参数确定方法、装置及计算机可读存储介质 |
CN110008494A (zh) * | 2018-12-13 | 2019-07-12 | 胜利油田森诺胜利工程有限公司 | 一种油田机械采油***能耗在线优化方法 |
CN110529079A (zh) * | 2019-08-14 | 2019-12-03 | 中国石油化工股份有限公司 | 一种采油产液量控制方法及装置 |
CN111475989A (zh) * | 2020-03-30 | 2020-07-31 | 中国石油天然气股份有限公司 | 一种基于油井示功图获取原油含水量的方法 |
CN111946331A (zh) * | 2020-08-20 | 2020-11-17 | 中联煤层气有限责任公司 | 一种井底流压的测试方法及获取粘滞阻力的方法 |
CN111963151A (zh) * | 2020-09-01 | 2020-11-20 | 中国石油天然气股份有限公司 | 一种通过抽油机悬点静载荷确定地层压力的方法 |
CN112392461A (zh) * | 2020-12-09 | 2021-02-23 | 中国石油天然气股份有限公司 | 一种快速计算油井井筒中混合液体含水率的方法 |
CN113338915A (zh) * | 2021-07-13 | 2021-09-03 | 西南石油大学 | 一种判断气井是否积液及预测积液高度的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102268989A (zh) * | 2010-06-07 | 2011-12-07 | 任丘市华北石油浩普科技有限公司 | 获取抽油机井动液面在线实时数据的方法以及*** |
CN102877832A (zh) * | 2012-09-14 | 2013-01-16 | 中国石油天然气股份有限公司 | 一种油井动液面的测量方法及*** |
WO2015143626A1 (zh) * | 2014-03-25 | 2015-10-01 | 中国石油天然气股份有限公司 | 一种基于示功图的二流量试井分析的方法及*** |
-
2015
- 2015-10-26 CN CN201510702552.1A patent/CN105257279A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102268989A (zh) * | 2010-06-07 | 2011-12-07 | 任丘市华北石油浩普科技有限公司 | 获取抽油机井动液面在线实时数据的方法以及*** |
CN102877832A (zh) * | 2012-09-14 | 2013-01-16 | 中国石油天然气股份有限公司 | 一种油井动液面的测量方法及*** |
WO2015143626A1 (zh) * | 2014-03-25 | 2015-10-01 | 中国石油天然气股份有限公司 | 一种基于示功图的二流量试井分析的方法及*** |
Non-Patent Citations (3)
Title |
---|
李娟: "稠油抽油机井功图量化分析技术研究", 《中国优秀硕士学位论文数据库》 * |
程曦月等: "有杆抽油***的数学建模及诊断", 《数学的实践与认识》 * |
***等: "《海洋石油工程》", 31 July 2015 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106351645A (zh) * | 2016-10-26 | 2017-01-25 | 北京安控科技股份有限公司 | 一种抽油机井动液面连续测量方法及装置 |
CN106351645B (zh) * | 2016-10-26 | 2020-01-03 | 北京安控科技股份有限公司 | 一种抽油机井动液面连续测量方法及装置 |
CN106777709B (zh) * | 2016-12-21 | 2018-11-06 | 中国石油大学(华东) | 基于地面实测电参数的潜油电泵井动液面计算的方法 |
CN106837305B (zh) * | 2016-12-28 | 2020-06-09 | 中国石油天然气股份有限公司 | 确定抽油井井下液面深度的方法和装置 |
CN106837305A (zh) * | 2016-12-28 | 2017-06-13 | 中国石油天然气股份有限公司 | 确定抽油井井下液面深度的方法和装置 |
CN107237615A (zh) * | 2017-08-09 | 2017-10-10 | 大连西贝通达科技有限公司 | 一种基于光杆载荷的游梁式抽油机运行速度优化控制方法 |
CN109779612A (zh) * | 2017-11-14 | 2019-05-21 | 中国石油化工股份有限公司 | 一种碳杆举升***井下泵功图测试装置及方法 |
CN109779612B (zh) * | 2017-11-14 | 2023-12-01 | 中国石油化工股份有限公司 | 一种碳杆举升***井下泵功图测试装置及方法 |
CN109899056A (zh) * | 2017-12-07 | 2019-06-18 | 中国石油天然气股份有限公司 | 抽油机的参数确定方法、装置及计算机可读存储介质 |
CN109899056B (zh) * | 2017-12-07 | 2021-11-30 | 中国石油天然气股份有限公司 | 抽油机的参数确定方法、装置及计算机可读存储介质 |
CN110008494A (zh) * | 2018-12-13 | 2019-07-12 | 胜利油田森诺胜利工程有限公司 | 一种油田机械采油***能耗在线优化方法 |
CN110529079A (zh) * | 2019-08-14 | 2019-12-03 | 中国石油化工股份有限公司 | 一种采油产液量控制方法及装置 |
CN110529079B (zh) * | 2019-08-14 | 2021-06-18 | 中国石油化工股份有限公司 | 一种采油产液量控制方法及装置 |
CN111475989A (zh) * | 2020-03-30 | 2020-07-31 | 中国石油天然气股份有限公司 | 一种基于油井示功图获取原油含水量的方法 |
CN111946331A (zh) * | 2020-08-20 | 2020-11-17 | 中联煤层气有限责任公司 | 一种井底流压的测试方法及获取粘滞阻力的方法 |
CN111946331B (zh) * | 2020-08-20 | 2024-05-24 | 中联煤层气有限责任公司 | 一种井底流压的测试方法及获取粘滞阻力的方法 |
CN111963151A (zh) * | 2020-09-01 | 2020-11-20 | 中国石油天然气股份有限公司 | 一种通过抽油机悬点静载荷确定地层压力的方法 |
CN111963151B (zh) * | 2020-09-01 | 2024-05-28 | 中国石油天然气股份有限公司 | 一种通过抽油机悬点静载荷确定地层压力的方法 |
CN112392461A (zh) * | 2020-12-09 | 2021-02-23 | 中国石油天然气股份有限公司 | 一种快速计算油井井筒中混合液体含水率的方法 |
CN112392461B (zh) * | 2020-12-09 | 2023-10-31 | 中国石油天然气股份有限公司 | 一种快速计算油井井筒中混合液体含水率的方法 |
CN113338915A (zh) * | 2021-07-13 | 2021-09-03 | 西南石油大学 | 一种判断气井是否积液及预测积液高度的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105257279A (zh) | 一种抽油机井动液面的测量方法 | |
Szilas | Production and transport of oil and gas | |
CN106437682B (zh) | 一种预测油井示功图的方法 | |
CN105089638A (zh) | 一种利用抽油机泵功图在线计算油井动液面方法 | |
CN109827734A (zh) | 一种评估内外流作用下深海立管涡激振动的方法 | |
CN108694254A (zh) | 一种变产变压生产气井产量经验递减曲线分析方法 | |
Liu et al. | Dynamic behavior of coalbed methane flow along the annulus of single-phase production | |
CN106351645A (zh) | 一种抽油机井动液面连续测量方法及装置 | |
CN109708707A (zh) | 一种气体流量测量装置及测量方法 | |
CN105178943A (zh) | 一种实时校正井筒压力的方法 | |
CN109162697A (zh) | 一种天然气井油套环空流动动态描述方法 | |
CN104776971A (zh) | 一种气流携液携砂可视化实验装置 | |
Baxendell | Producing wells on casing flow-an analysis of flowing pressure gradients | |
CN102587887B (zh) | 一种气井井况预测方法 | |
Tan et al. | Review of variable speed drive technology in beam pumping units for energy-saving | |
CN109815543A (zh) | 计算气田动态储量的方法 | |
CN105224716A (zh) | 一种基于泵功图载荷变化的油井结蜡预警方法 | |
NO343700B1 (no) | Fremgangsmåte og apparat for å måle volumstrømningsraten for en væske | |
CN109915128A (zh) | 地层承压能力动态测试方法及固井方法 | |
CN106593408A (zh) | 油井流速的获取方法与装置 | |
CN108021769A (zh) | 一种判断气井是否正常携液的方法 | |
CN101871344B (zh) | 一种称重式气井井筒液面位置的确定方法 | |
RU2372480C1 (ru) | Способ определения дебита нефтяной скважины | |
CN107939367B (zh) | 一种压裂水马力确定方法 | |
US6863128B2 (en) | Method of predicting friction pressure drop of proppant-laden slurries using surface pressure data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160120 |
|
RJ01 | Rejection of invention patent application after publication |