一种基于油井示功图的能量计算动液面方法与流程

[0001]
本发明属于油田数字化及采油工艺技术领域，具体涉及一种有杆抽油示功图获取系统及计算动液面方法。


背景技术：

[0002]
追求产量和综合效率最大化是油田开发的终极目标，而地层的供液能力将是制约实现这一目标的根本因素，所以了解掌握地层供液能力并加以科学地利用，就成为油田开发管理活动中不可回避的课题。而油井动液面是反应地层供液能力的重要指标，是进行机采井设计的重要依据。因此，及时准确地获得动液面将是非常重要和有意义的
[0003]
长久以来，抽油机井环空动液面都是使用回声仪测试的，测试方法是人工定期进行操作，使用的设备比较危险、笨重，而且也很难长时间连续测试。目前所采用的动液面测试仪器由于受油井套管环空内死油环、井下狗腿等因素的影响，造成测试分析结果与实际液面情况不符，目前也有使用其他测量方法(诸如电动气枪、电动氮气瓶)，但因工艺结构复杂，成本较高，使用周期短，因此推广起来比较困难。所以有必要寻求新的动液面获取方法，以实现油田开发降本增效之目的。


技术实现要素：

[0004]
本发明提供了一种有杆抽油示功图获取系统及计算动液面方法，目的在于通过采集功图数据，结合已有分析结果，基于油井示功图的能量法计算动液面方法，以达到指导生产，降本增效的目的。
[0005]
为了实现上述目的，本发明采用的技术方案是：
[0006]
一种基于油井示功图的能量计算动液面方法，包括如下步骤，
[0007]
步骤一：通过油井有杆抽油功图系统获取数据，得到油井示功图；
[0008]
步骤二：建立数学模型，并通过建立的数学模型及步骤一得到的油井示功图中参数，计算得到油井泵功图；
[0009]
步骤三：利用步骤二得到的油井泵功图中参数数据，求取泵有效冲程；
[0010]
步骤四：计算泵功图面积；
[0011]
步骤五：通过步骤四计算得到的泵功图面积，计算得到油井动液面参数值；
[0012]
步骤六：将得到的油井动液面参数值，用于油田生产的指导。
[0013]
所述的步骤一中建立数学模型并通过建立的数学模型计算得到油井泵功图的方法如下：
[0014]
第一步：建立有杆泵采油系统的故障诊断简化模型；
[0015]
设抽油杆柱的位移s向下为正，向上为负；抽油杆柱的轴力f拉为正，压为负；把整个抽油杆柱分为m个单元，对每个单元做如下的等效处理：去掉单元的弯曲刚度，将单元的质量集中到单元一端的结点上，形成多个集中质量m
i
，各个集中质量之间用拉压刚度为k
i
的弹簧来代替，其影响等效成阻尼系数为c
i
的阻尼器，设抽油杆和油管之间的摩擦系数为μ
i
；
[0016]
拉压刚度为k
i
的弹簧两端的作用力都为f
i-1
；集中质量m
i
两端的位移都为s
i
，于是有：
[0017]
m
i
＝q
ri
l
i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0018]
式中：
[0019]
q
ri
为第i个单元的单长杆重，kg/m；
[0020]
l
i
为第i个单元的长度，m；
[0021][0022]
式中，
[0023]
e
i
，第i个单元材料的弹性模量pa；
[0024]
a
i
，第i个单元的横截面面积，m2；
[0025]
第二步：推导抽油杆柱各截面示功图的递推格式；
[0026]
先考虑拉压刚度为k
i
的弹簧，由胡克定律知：
[0027]
f
i-1
＝k
i
(s
i-s
i-1
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0028]
解出：
[0029][0030]
再考虑集中质量m
i
，由牛顿第二定律知：
[0031][0032]
将式(4)代入式(5)并进行整理，得
[0033][0034]
将式(4)和式(6)统一记作：
[0035][0036]
此式即为求解抽油杆柱各截面示功图的递推格式；
[0037]
第三步：快速传递矩阵算法
[0038]
式(7)中的s
i-1
＝s
i-1
(θ)、f
i-1
＝f
i-1
(θ)、s
i
＝s
i
(θ)和f
i
＝f
i
(θ)都是曲柄转角θ的函数，把它们展开为fourier级数：
[0039]
[0040][0041][0042][0043]
其中：
[0044]
i指的是第i个抽油杆单元，i＝{0、1
……
m}，m是抽油杆最大个数；
[0045]
j指的是数据对个数，每个单元的位移与载荷是由j个数据对构成的， j＝{0、1
……
n}，n是最大个数；
[0046][0047][0048][0049][0050]
式中：
[0051]
k——2π周期内等分的份数，份数不少于144；
[0052]
p——从0到k各点序号；
[0053]
把上面(8)、(9)、(10)和(11)四个fourier级数代入式(7)，并设曲柄匀速转动，(n为冲次，min-1
)，整理可得：
[0054][0055]
[0056][0057]
令
[0058][0059]
有
[0060][0061][0062][0063]
简写为
[0064][0065][0066][0067]
其中：
[0068]
ω—曲柄角速度；
[0069]
m—抽油杆柱第m个单元；
[0070]
悬点位置的位移函数s0＝s0(θ)＝-s
pr
(θ)与轴力函数由实测光杆示功图求得，为已知；把它们展开为fourier级数，获得系数和然
后通过式(19)、(20)和(21)获得泵处的系数和最后通过fourier级数式(8)～(11) 的计算就可以求出泵处的位移函数s
pump
(θ)＝-s
m
(θ)＝-s
m
和轴力函数 f
pump
(θ)＝f
m
(θ)＝f
m
，即得到油井泵示功图。
[0071]
所述的步骤二中的油井示功图中参数包括油井悬点载荷和位移。
[0072]
所述的步骤三求取泵有效冲程具体过程如下：
[0073]
第一步：将步骤二得到的油井泵功图与数据库中的常见故障泵功图进行比较；
[0074]
第二步：通过第一步比较得出的结果，寻找出故障；
[0075]
第三步：结合生产参数确定泵有效冲程。
[0076]
所述的步骤四中计算泵功图面积是采用下述积分式得到的
[0077][0078]
其中：
[0079]
s是泵功图面积；
[0080]
x是拉杆与泵连接处位移；
[0081]
f
u
是上冲程泵筒所受载荷；
[0082]
f
d
是下冲程泵筒所受载荷。
[0083]
所述的步骤五计算油井动液面参数值的方法如下：
[0084]
液体据升到地面做的功w
[0085][0086]
其中，ρ
液
—油井产液密度；
[0087]
d—泵筒内径；
[0088]
h—有效扬程，m；
[0089]
s
e
—为泵有效冲程；
[0090]
g—为重力加速度；
[0091]
其中，有效扬程由下式求得
[0092][0093]
式中，h——动液面深度，m，
[0094]
p
c
——油管压力，mpa；
[0095]
p
t
——套管压力，mpa。
[0096]
液体密度由下式求得
[0097]
ρ
液
＝(1-f
w
)ρ
o
+f
w
ρ
w
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(25)
[0098]
式中，f
w
——含水率，
[0099]
ρ
o
——油的密度，t/m3，
[0100]
ρ
w
——水的密度，t/m3，
[0101]
由式s＝w
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26)
[0102]
可得
[0103]
所述步骤一中的油井有杆抽油功图系统包括数据采集点、数据处理点和服务器；所述的数据采集点与数据处理点之间电信号连接；所述的数据处理点和服务器之间电信号连接。
[0104]
所述数据采集点包括抽油机、电参数传感器、位移传感器、载荷传感器、rtu和电机监测装置；所述的电参数传感器设置在抽油机的电控箱上，所述的位移传感器设置在抽油机的游梁中轴上方；所述的载荷传感器设置在抽油机的悬绳器与光杆卡子上；所述的电参数传感器、位移传感器和载荷传感器分别与rtu电信号连接，rtu与数据处理点电信号连接；所述的电机监测装置分别与抽油机和rtu电连接。
[0105]
所述数据处理点包括数据接收器、数据发送器和数据处理器，所述的数据接收器与数据采集点电信号连接，数据处理器分别与数据接收器和数据发送器电信号连接，所述的数据发送器与服务器电信号连接。
[0106]
有益效果：
[0107]
本发明通过已采集功图数据，结合已有分析结果，基于油井示功图的能量法计算动液面方法，达到了指导生产、降本增效的目的。
[0108]
上述说明仅是本发明技术方案的概述，为了能够更清楚的了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例，详细说明如后。
附图说明
[0109]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0110]
图1为本发明有杆泵采油系统的故障诊断简化模型示意图；
[0111]
图2为本发明抽油机传感器安装示意图；
[0112]
图3为本发明有杆抽油功图系统结构图；
[0113]
图4为本发明的泵有效冲程识别流程图；
[0114]
图5为本发明采集到的光杆示功图和据此分析得到的泵示功图；
[0115]
图6为本发明抽油机井抽油井筒示意图；
[0116]
图7为本发明工况分析与功图计产系统。
[0117]
图中：1-电参数传感器；2-位移传感器；3-载荷传感器；4-数据采集点；5-数据处理点；6-服务器；7-抽油机。
[0118]
上述说明仅是本发明技术方案的概述，为了能够更清楚的了解本发明的技术手段，并可依照说明书的内容予以实施，以下以本发明的较佳实施例，详细说明如后。
具体实施方式
[0119]
下面将结合实施例，对本发明的技术方案进行清楚、完整地描述，显然，所描述的
实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0120]
实施例一：
[0121]
一种基于油井示功图的能量计算动液面方法，包括如下步骤，
[0122]
步骤一：通过油井有杆抽油功图系统获取数据，得到油井示功图；
[0123]
步骤二：建立数学模型，并通过建立的数学模型及步骤一得到的油井示功图中参数，计算得到油井泵功图；
[0124]
步骤三：利用步骤二得到的油井泵功图中参数数据，求取泵有效冲程；
[0125]
步骤四：计算泵功图面积；
[0126]
步骤五：通过步骤四计算得到的泵功图面积，计算得到油井动液面参数值；
[0127]
步骤六：将得到的油井动液面参数值，用于油田生产的指导。
[0128]
在实际使用时，本发明通过六个步骤，将已采集功图数据，结合已有分析结果，基于油井示功图的能量法计算动液面方法，达到了指导生产、降本增效的目的。
[0129]
本发明的方法简单易行，且不增加现有生产成本。
[0130]
实施例二：
[0131]
参照图1、图4-图7所示的一种基于油井示功图的能量计算动液面方法，在实施例一的基础上，所述的步骤一中建立数学模型并通过建立的数学模型计算得到油井泵功图的方法如下：
[0132]
第一步：建立有杆泵采油系统的故障诊断简化模型；
[0133]
设抽油杆柱的位移s向下为正，向上为负；抽油杆柱的轴力f拉为正，压为负；把整个抽油杆柱分为m个单元，对每个单元做如下的等效处理：去掉单元的弯曲刚度，将单元的质量集中到单元一端的结点上，形成多个集中质量m
i
，各个集中质量之间用拉压刚度为k
i
的弹簧来代替，其影响等效成阻尼系数为c
i
的阻尼器，设抽油杆和油管之间的摩擦系数为μ
i
；(如图 1)
[0134]
拉压刚度为k
i
的弹簧两端的作用力都为f
i-1
；集中质量m
i
两端的位移都为s
i
，于是有：
[0135]
m
i
＝q
ri
l
i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0136]
式中：
[0137]
q
ri
为第i个单元的单长杆重，kg/m；
[0138]
l
i
为第i个单元的长度，m；
[0139][0140]
式中，
[0141]
e
i
，第i个单元材料的弹性模量pa；
[0142]
a
i
，第i个单元的横截面面积，m2；
[0143]
第二步：推导抽油杆柱各截面示功图的递推格式；
[0144]
先考虑拉压刚度为k
i
的弹簧，由胡克定律知：
[0145]
f
i-1
＝k
i
(s
i-s
i-1
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0146]
解出：
[0147][0148]
再考虑集中质量m
i
，由牛顿第二定律知：
[0149][0150]
将式(4)代入式(5)并进行整理，得
[0151][0152]
将式(4)和式(6)统一记作：
[0153][0154]
此式即为求解抽油杆柱各截面示功图的递推格式；
[0155]
第三步：快速传递矩阵算法
[0156]
式(7)中的s
i-1
＝s
i-1
(θ)、f
i-1
＝f
i-1
(θ)、s
i
＝s
i
(θ)和f
i
＝f
i
(θ)都是曲柄转角θ的函数，把它们展开为fourier级数：
[0157][0158][0159][0160][0161]
其中：
[0162]
i指的是第i个抽油杆单元，i＝{0、1
……
m}，m是抽油杆最大个数；
[0163]
j指的是数据对个数，每个单元的位移与载荷是由j个数据对构成的， j＝{0、1
……
n}，n是最大个数；
[0164][0165][0166][0167][0168]
式中：
[0169]
k——2π周期内等分的份数，份数不少于144；
[0170]
p——从0到k各点序号；
[0171]
把上面(8)、(9)、(10)和(11)四个fourier级数代入式(7)，并设曲柄匀速转动，(n为冲次，min-1
)，整理可得：
[0172][0173][0174][0175]
令
[0176][0177]
有
[0178][0179][0180][0181]
简写为
[0182][0183][0184][0185]
其中：
[0186]
ω—曲柄角速度；
[0187]
m—抽油杆柱第m个单元；
[0188]
悬点位置的位移函数s0＝s0(θ)＝-s
pr
(θ)与轴力函数由实测光杆示功图求得，为已知；把它们展开为fourier级数，获得系数和然后通过式(19)、(20)和(21)获得泵处的系数和最后通过fourier级数式(8)～(11) 的计算就可以求出泵处的位移函数s
pump
(θ)＝-s
m
(θ)＝-s
m
和轴力函数 f
pump
(θ)＝f
m
(θ)＝f
m
，即得到油井泵示功图。
[0189]
进一步的，所述的步骤二中的油井示功图中参数包括油井悬点载荷和位移。
[0190]
进一步的，所述的步骤三求取泵有效冲程具体过程如下(如图4)：
[0191]
第一步：将步骤二得到的油井泵功图与数据库中的常见故障泵功图进行比较；
[0192]
第二步：通过第一步比较得出的结果，寻找出故障；
[0193]
第三步：结合生产参数确定泵有效冲程。
[0194]
进一步的，所述的步骤四中计算泵功图面积是采用下述积分式得到的 (如图5)
[0195]
[0196]
其中：
[0197]
s是泵功图面积；
[0198]
x是拉杆与泵连接处位移；
[0199]
f
u
是上冲程泵筒所受载荷；
[0200]
f
d
是下冲程泵筒所受载荷。
[0201]
进一步的，所述的步骤五计算油井动液面参数值的方法如下：
[0202]
液体据升到地面做的功w(如图6)
[0203][0204]
其中，ρ
液
—油井产液密度；
[0205]
d—泵筒内径；
[0206]
h—有效扬程，m；
[0207]
s
e
—为泵有效冲程；
[0208]
g—为重力加速度；
[0209]
其中，有效扬程由下式求得
[0210][0211]
式中，h——动液面深度，m，
[0212]
p
c
——油管压力，mpa；
[0213]
p
t
——套管压力，mpa。
[0214]
液体密度由下式求得
[0215]
ρ
液
＝(1-f
w
)ρ
o
+f
w
ρ
w
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(25)
[0216]
式中，f
w
——含水率，
[0217]
ρ
o
——油的密度，t/m3，
[0218]
ρ
w
——水的密度，t/m3，
[0219]
由式s＝w
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(26)
[0220]
可得
[0221]
其中的s
e
可由已有工况分析与功图计产系统获得(如图7)。
[0222]
采用本发明的技术方案，在不增加现有生产成本的基础上，简单易行的得到了动液面参数值，达到了指导生产，降本增效的目的。
[0223]
实施例三：
[0224]
参照图2和图3所示的一种基于油井示功图的能量计算动液面方法，所述步骤一中的油井有杆抽油功图系统包括数据采集点4、数据处理点5和服务器6；所述的数据采集点4与数据处理点5之间电信号连接；所述的数据处理点5和服务器6之间电信号连接。
[0225]
进一步的，所述数据采集点4包括抽油机7、电参数传感器1、位移传感器2、载荷传感器4、rtu和电机监测装置；所述的电参数传感器1设置在抽油机7的电控箱上，所述的位移传感器2设置在抽油机7的游梁中轴上方；所述的载荷传感器4设置在抽油机7的悬绳器与光
杆卡子上；所述的电参数传感器1、位移传感器2和载荷传感器4分别与rtu电信号连接， rtu与数据处理点5电信号连接；所述的电机监测装置分别与抽油机7和 rtu电连接。
[0226]
进一步的，所述数据处理点5包括数据接收器、数据发送器和数据处理器，所述的数据接收器与数据采集点4电信号连接，数据处理器分别与数据接收器和数据发送器电信号连接，所述的数据发送器与服务器6电信号连接。
[0227]
在具体应用时，本发明通过安装在井口的载荷传感器4、位移传感器2，每10分钟获取一张油井示功图，然后建立数学模型计算得到油井泵功图，由泵功图求取泵有效冲程及其面积，又由于泵功图的面积就相当于把本次冲程泵筒内的液量由动液面处举升到地面所做的功，最后实现实时获取泵功图面积，进而可得到油井动液面，最后将得到的油井动液面参数值，用于油田生产的指导。
[0228]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
[0229]
在不冲突的情况下，本领域的技术人员可以根据实际情况将上述各示例中相关的技术特征相互组合，以达到相应的技术效果，具体对于各种组合情况在此不一一赘述。
[0230]
需要说明，本发明实施例中所有方向性指示(诸如上、下、左、右、前、后
……
)仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
[0231]
以上所述，只是本发明的较佳实施例而已，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖性特点相一致的最宽的范围。依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。