CN106777709B - 基于地面实测电参数的潜油电泵井动液面计算的方法 - Google Patents
基于地面实测电参数的潜油电泵井动液面计算的方法 Download PDFInfo
- Publication number
- CN106777709B CN106777709B CN201611194201.5A CN201611194201A CN106777709B CN 106777709 B CN106777709 B CN 106777709B CN 201611194201 A CN201611194201 A CN 201611194201A CN 106777709 B CN106777709 B CN 106777709B
- Authority
- CN
- China
- Prior art keywords
- pump
- temperature
- pressure
- well
- section
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000012530 fluid Substances 0.000 claims abstract description 108
- 238000004364 calculation method Methods 0.000 claims abstract description 58
- 239000007788 liquid Substances 0.000 claims abstract description 45
- 238000012546 transfer Methods 0.000 claims abstract description 30
- 238000010438 heat treatment Methods 0.000 claims abstract description 12
- 230000005855 radiation Effects 0.000 claims abstract description 5
- 238000004519 manufacturing process Methods 0.000 claims description 17
- 238000005086 pumping Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 230000001012 protector Effects 0.000 claims description 5
- 239000004568 cement Substances 0.000 claims description 4
- 238000010792 warming Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000003068 static effect Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 5
- 238000005457 optimization Methods 0.000 abstract description 3
- 239000003129 oil well Substances 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005499 meniscus Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/04—Measuring depth or liquid level
- E21B47/047—Liquid level
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geochemistry & Mineralogy (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
一种基于地面实测电参数的潜油电泵井动液面计算的方法，考虑油套环空中液面以上气体热辐射、液面以下液体热传导，且电机、电缆发热等影响，构建了“四段法”潜油电泵井井筒流体温度计算模型；研究分析了实时监测的地面电参数与潜油电机输入功率、泵扬程、泵吸入口压力的内在关系以及潜油电泵的工作特性，推导了基于地面实测电参数的潜油电泵井动液面计算模型。本发明用于潜油电泵井动液面预测，为潜油电泵井智能化管理和实时优化提供理论和技术支持。
Description
技术领域
本发明涉及一种基于地面实测电参数的潜油电泵井动液面计算的方法，属于计算机油田信息化建设的技术领域。
背景技术
潜油电泵举升工艺是油井保持或提高产量的主要途径，其动液面是判别油井供液能力和分析油井生产状况的重要资料。目前，潜油电泵井的地面电参数和井口油压、回压、套压、温度等资料已经实现了实时监测，但动液面无法实时监测，仍然采用声波法，由人工定期监测，不能满足油井生产智能分析和实时优化的需要，影响了油井生产系统效率和效益。因此，推导建立基于地面实测电参数的潜油电泵井动液面计算方法，实现了潜油电泵井动液面的自动实时监测，有利于油井的实时工况分析和智能化、信息化管理水平的提高，为潜油电泵井实时优化生产、提高运行效率与效益提供理论和技术支持。
发明内容
针对现有技术的不足，本发明提出一种基于地面实测电参数的潜油电泵井动液面计算的方法。本发明的目的是计算潜油电泵井的动液面，以解决潜油电泵井动液面无法实时监测的问题，实现油井智能化管理与实时优化运行，提高油井生产效率和效益。
本发明的技术方案如下：
基于地面实测电参数的潜油电泵井动液面计算的方法，包括如下步骤，
步骤一：利用“四段法”建立潜油电泵井井筒流体温度计算模型；
步骤二：建立基于地面电参数的潜油电泵所需输入功率计算模型：
步骤三：建立潜油电泵扬程计算模型：
步骤四：建立泵吸入口压力模型：
步骤五：潜油电泵井生产时，井口到管鞋的油套环空中流体分为两段，即气柱段和油柱段,由气柱段压力计算模型和液面以下油柱段压力计算模型联合求解，即确定动液面深度。
根据本发明优选的，所述步骤一：利用“四段法”建立潜油电泵井井筒流体温度计算模型：
(a)综合考虑油套环空液面以上气体热辐射、液面以下液体热传导以及电机、电缆加热等影响，提出合理的模型假设条件，确定井筒内产液温度沿井深的变化规律：
井筒内流体温度分布规律由方程(I)描述，
采用常数变易法求解，即：
式(I)、(II)中，T为井筒内某点处的产液温度，℃；Tc为井口产出液温度，℃；Ts为地表恒温层温度，℃；kl为传热系数，W/(m·K)；G为流体质量流量，kg/s；cp为流体比定压热容，J/(kg·K)；x为以井口为原点，所述计算点的深度，m；m为地温梯度，℃/m；qv为内热源强度，W/m；g为重力加速度，m/s2；
(b)根据潜油电泵井井筒结构和动液面位置，将井筒分为四段：即井口至动液面段、动液面至泵出口段、泵机组段和电机至井底段，各段的温度计算模型为：
(1)井口至动液面段温度计算模型为
(2)动液面至泵出口段温度计算模型为
(3)泵机组段温度计算模型为
Tpint＝Tpout-ΔTdj-ΔTdl(V)
(4)电机至井底段温度计算模型为
式(III)-(VI)中，kl1为井口至动液面段内的井筒与地层之间的传热系数，W/(m·K)；Td为动液面处的井筒流体温度，℃；Hd为动液面深度，m；kl2为动液面至泵出口段内的井筒与地层之间的传热系数，W/(m·K)；Tpint为泵吸入口处流体的温度，℃；ΔTdj为电机增温，℃；ΔTdl为泵内小扁电缆的发热增温，℃；Hp为下泵深度，m；kl3为电机至井底段内的井筒与地层之间的传热系数，W/(m·K)。
根据本发明优选的，所述步骤二：建立基于地面电参数的潜油电泵所需输入功率计算模型：
式(VII)中，Np为潜油电泵输入功率，kW；U为地面电压，V；I为地面电流，A；cosφ为功率因数，小数；Rc为小扁电缆电阻，Ω/m；Lp为下放潜油电泵所需的电缆长度，m；Ls为电泵内小扁电缆的长度m；ηm为潜油电机工作时的泵效，小数；P分离器和P保护器分别为分离器和保护器所消耗的功率(一般地，每级分离器的机械损耗功率为1～1.5kW，每级保护器的机械损耗功率为1.0kW)，kW。
根据本发明优选的，所述步骤三：建立潜油电泵扬程计算模型：
式(VIII)中，H为潜油电泵有效举升高度，m；ηo为潜油电泵实际泵效，小数；Q为泵的实际排量，m3/d；γl为井液平均相对密度，小数；
根据潜油电泵工作特性曲线，其实际泵效由式(IX)表示：
η0＝kη[anQn+…+a2Q2+a1Q+a0] (IX)
式(IX)中，kη为泵效的修正系数，小数；an，···，a2，a1，a0为潜油电泵特性曲线泵效与排量曲线回归系数，常数。
根据本发明优选的，所述步骤四：建立泵吸入口压力模型：
以泵为节点，计算出泵吸入口的压力，即沉没压力，
pin＝pout-ρlgH (X)
式(X)中，pin为泵吸入口处的压力，MPa；pout为泵排出口处的压力(由Beggs-Brill方法确定)，MPa；ρl为泵内流体的混合密度，kg/m3。
根据本发明优选的，所述步骤五：潜油电泵井生产时，井口到管鞋的油套环空中流体分为两段，即气柱段和油柱段，以Δh为计算单元步长，将井筒深度划分为个网格，计算每个网格的压力，以提高模型求解精度，由气柱段压力计算模型和液面以下油柱段压力计算模型联合求解，即确定动液面深度：
式(XI)中，pg(i)为第i网格深度处的气柱压力，MPa；pg(i-1)为第i-1网格深度处的气柱压力，MPa；ρg0为标准状况下的气体密度，kg/m3；T0为标准状况下的温度，K；p0为标准状况下的压力，MPa；Δh为网格划分单元格的精度，m；i为网格划分的网格数，整数；Tav为第i-1和第i深度该段Δh内的平均温度，K；Zav为第i-1和第i深度该段Δh内的平均温度和平均压力下的气体压缩因子，小数；po(x)为x深度处的油柱压力，MPa；ρo(x)为微元段dx内平均压力和平均温度下的油的密度，kg/m3；pc为井口套管压力，MPa；
根据上述模型，计算流程如下：
①以计算单元步长Δh将井筒网格划分，初设一个动液面假设值Hd，计算井筒流体温度和压力分布(以井口为起点)，得到泵排出口压力pout；
②利用公式(VII)，通过地面实测电参数，计算潜油电泵所需的轴功率Np；
③利用公式(VIII)，根据潜油电泵的工作特性，计算有效举升高度H；
④利用公式(X)，以泵为节点，得到泵吸入口压力pin；
⑤利用公式(XI)所述的气柱段压力计算模型和液面以下油柱段压力计算模型联合求解动液面Hd’；
⑥比较Hd’与Hd，若|Hd’-Hd|＜ε，输出结果Hd’；若不满足上述条件，将令Hds＝0.618Hd’+0.382Hd，将Hds作为新的动液面假设值，重复步骤①～⑤计算新的动液面计算值Hd’，直至满足|Hd’-Hds|＜ε要求，输出新的动液面计算值Hd’。
本发明的技术优势在于：
相对于人工定期的回声仪测量，本发明的优点在于：
(1)本发明只需根据地面实测电参数和生产数据，实现潜油电泵井动液面的自动实时监测，且具有较高的精度；
(2)提高油井实时工况分析和智能化、信息化管理水平；
(3)降低劳动强度，减少用工量、提升生产效率、节约生产成本。
附图说明
图1是本发明所述基于电参数的电潜泵井动液面计算方法的流程图；
图2是17口油井动液面计算结果对比图。
具体实施方式
下面结合实施例和说明书附图对本发明做详细的说明，但不限于此。
实施例1、
步骤一：利用“四段法”建立潜油电泵井井筒流体温度计算模型
电潜泵井生产过程中，油套环空中液面以上气体热辐射、液面以下液体热传导，且电机、电缆设备发热等，影响井筒流体温度及其流动规律。基于电潜泵井实际物理模型，提出如下假设：
(1)传热过程为井筒内的稳态传热过程和地层内的非稳态传热过程；
(2)井内流动为一维稳定流动；
(3)井筒周围地层温度，受地温梯度影响，并且随着深度线性变化；
(4)电缆为均匀发热的线热源，潜油电机为点热源；
(5)考虑实际生产时套管产气量很小，认为环空流体为静止状态，忽略环空流体对流传热影响；
(6)以井口为坐标原点，沿井筒向下为坐标轴的正方向。
(a)综合考虑油套环空液面以上气体热辐射、液面以下液体热传导以及电机、电缆加热等影响，提出合理的模型假设条件，确定井筒内产液温度沿井深的变化规律：
井筒内流体温度分布规律由方程(I)描述，
采用常数变易法求解，即：
式(I)、(II)中，T为井筒内某点处的产液温度，℃；Tc为井口产出液温度，℃；Ts为地表恒温层温度，℃；kl为传热系数，W/(m·K)；G为流体质量流量，kg/s；cp为流体比定压热容，J/(kg·K)；x为以井口为原点，所述计算点的深度，m；m为地温梯度，℃/m；qv为内热源强度，W/m；g为重力加速度，m/s2；
(b)根据潜油电泵井井筒结构和动液面位置，将井筒分为四段：即井口至动液面段、动液面至泵出口段、泵机组段和电机至井底段，各段的温度计算模型为：
(1)井口至动液面段温度计算模型为
其中，
kl1为井口至动液面段内的井筒与地层之间的传热系数，W/(m·K)；hl为产液与油管内表面的表面传热系数，W/(m·K)；dti为油管内径，m；λtub为油管的导热系数，W/(m·K)；dto为油管外径，m；λr为油管外表面与套管内表面间辐射传热过程的当量导热系数(由于油管外表温度和套管内表面温度未知，为计算传热系数kl，需采用迭代法求解)，W/(m·K)；λcas为套管的导热系数，W/(m·K)；dco为套管外径，m；dci为套管内径，m；λcem为水泥环的导热系数，W/(m·K)；dh为水泥环外缘直径，m；Rle为单位长度地层热阻，(m·K)/W；σb为Stefan-Boltzmann常数，取值为5.763×10-8W/(m2·K4)；Tto为油管外表面温度，K；Tci为油管外表面温度，K；εc，εt分别为套管内表面和油管外表面的发射率，小数；I为泵机组工作时的工作电流，A；R0为大扁电缆单位长度的电阻值，Ω/m。
动液面处的井筒流体温度为：
Td为动液面处的井筒流体温度，℃；Hd为动液面深度，m；
(2)动液面至泵出口段温度计算模型为
其中，
式中，kl2为动液面至泵出口段内的井筒与地层之间的传热系数，W/(m·K)；λl为流体导热系数，W/(m·K)。
泵出口处的井筒流体温度为：
式中，Tpout为泵排出口处流体的温度，℃；Hp为下泵深度，m。
(3)泵机组段温度计算模型为
Tpint＝Tpout-ΔTdj-ΔTdl (V)
其中，
式中，Tpint为泵吸入口处流体的温度，℃；ΔTdj为电机增温，℃；ΔTdl为泵内小扁电缆的发热增温，℃；Nm为电机功率，kW；ηm为电机效率，小数；Rc为小扁电缆单位长度的电阻值，Ω/m；Ls为泵内小扁电缆长度，m。
(4)电机至井底段温度计算模型为
其中，
式中，kl3为电机至井底段内的井筒与地层之间的传热系数，W/(m·K)。
式(III)-(VI)中，kl1为井口至动液面段内的井筒与地层之间的传热系数，W/(m·K)；Td为动液面处的井筒流体温度，℃；Hd为动液面深度，m；kl2为动液面至泵出口段内的井筒与地层之间的传热系数，W/(m·K)；Tpint为泵吸入口处流体的温度，℃；ΔTdj为电机增温，℃；ΔTdl为泵内小扁电缆的发热增温，℃；Hp为下泵深度，m；kl3为电机至井底段内的井筒与地层之间的传热系数，W/(m·K)；
步骤二：建立基于地面电参数的潜油电泵所需输入功率计算模型：
式(VII)中，Np为潜油电泵输入功率，kW；U为地面电压，V；I为地面电流，A；cosφ为功率因数，小数；Rc为小扁电缆电阻，Ω/m；Lp为下放潜油电泵所需的电缆长度，m；Ls为电泵内小扁电缆的长度m；ηm为潜油电机工作时的泵效，小数；P分离器和P保护器分别为分离器和保护器所消耗的功率(一般地，每级分离器的机械损耗功率为1～1.5kW，每级保护器的机械损耗功率为1.0kW)，kW；
步骤三：建立潜油电泵扬程计算模型：
式(VIII)中，H为潜油电泵有效举升高度，m；ηo为潜油电泵实际泵效，小数；Q为泵的实际排量，m3/d；γl为井液平均相对密度，小数；
根据潜油电泵工作特性曲线，其实际泵效由式(IX)表示：
η0＝kη[anQn+…+a2Q2+a1Q+a0] (IX)
式(IX)中，kη为泵效的修正系数，小数；an，···，a2，a1，a0为潜油电泵特性曲线泵效与排量曲线回归系数，常数。
步骤四：建立泵吸入口压力模型：
以泵为节点，计算出泵吸入口的压力，即沉没压力，
pin＝pout-ρlgH (X)
式(X)中，pin为泵吸入口处的压力，MPa；pout为泵排出口处的压力(由Beggs-Brill方法确定)，MPa；ρl为泵内流体的混合密度，kg/m3。
步骤五：潜油电泵井生产时，井口到管鞋的油套环空中流体分为两段，即气柱段和油柱段，以Δh为计算单元步长，将井筒深度划分为个网格，计算每个网格的压力，以提高模型求解精度。由气柱段压力计算模型和液面以下油柱段压力计算模型联合求解，即确定动液面深度：
式(XI)中，pg(i)为第i网格深度处的气柱压力，MPa；pg(i-1)为第i-1网格深度处的气柱压力，MPa；ρg0为标准状况下的气体密度，kg/m3；T0为标准状况下的温度，K；p0为标准状况下的压力，MPa；Δh为网格划分单元格的精度，m；i为网格划分的网格数，整数；Tav为第i-1和第i深度该段Δh内的平均温度，K；Zav为第i-1和第i深度该段Δh内的平均温度和平均压力下的气体压缩因子，小数；po(x)为x深度处的油柱压力，MPa；ρo(x)为微元段dx内平均压力和平均温度下的油的密度，kg/m3；pc为井口套管压力，MPa；
根据上述模型与研究，计算流程如下：
①以计算单元步长Δh将井筒网格划分，假设动液面Hd，计算井筒流体温度和压力分布(以井口为起点)，得到泵排出口压力pout；
②利用公式(VII)，通过地面实测电参数，计算潜油电泵所需的轴功率Np；
③利用公式(VIII)，根据潜油电泵的工作特性，计算有效举升高度H；
④利用公式(X)，以泵为节点，得到泵吸入口压力pin；
⑤利用公式(XI)所述的气柱段压力计算模型和液面以下油柱段压力计算模型联合求解动液面Hd’；
⑥比较Hd’与Hd，若|Hd’-Hd|＜ε，输出结果Hd’；若不满足上述条件，将令Hds＝0.618Hd’+0.382Hd，将Hds作为新的动液面假设值，重复
步骤①～⑤计算新的动液面计算值Hd’，直至满足|Hd’-Hds|＜ε要求，
输出新的动液面计算值Hd’。
经现场17口油井资料，计算与分析可知(见表1，图2)：①基于“四段法”的井筒流体温度计算的最大相对误差为17.70％，最小的相对误差为3.39％，其平均相对误差为8.75％；②动液面计算的最大相对误差为14.61％，最小相对误差为1.65％，其平均相对误差为6.98％。结果表明：动液面计算模型具有较好的精度。
表1井筒油层流体及动液面计算结果
本发明所述的基于地面实测电参数的潜油电泵井动液面的计算方法，能够根据地面电参数及日常生产数据，实现实时、准确的计算单井动液面深度，解决动液面实时监测的困难。
本发明应用于潜油电泵井日常工况参数检测领域，能够有效地配合对潜油电泵井的工况分析，为潜油电泵井智能化管理和实时优化调整提供理论和技术支持。
Claims (1)
1.基于地面实测电参数的潜油电泵井动液面计算的方法，其特征在于，该方法包括如下步骤，
步骤一：利用“四段法”建立潜油电泵井井筒流体温度计算模型；
步骤二：建立基于地面电参数的潜油电泵所需输入功率计算模型：
步骤三：建立潜油电泵扬程计算模型：
步骤四：建立泵吸入口压力模型：
步骤五：潜油电泵井生产时，井口到管鞋的油套环空中流体分为两段，即气柱段和油柱段,由气柱段压力计算模型和液面以下油柱段压力计算模型联合求解，即确定动液面深度；
所述步骤一：利用“四段法”建立潜油电泵井井筒流体温度计算模型：
基于电潜泵井实际物理模型，令：
(1)传热过程为井筒内的稳态传热过程和地层内的非稳态传热过程；
(2)井内流动为一维稳定流动；
(3)井筒周围地层温度，受地温梯度影响，并且随着深度线性变化；
(4)电缆为均匀发热的线热源，潜油电机为点热源；
(5)考虑实际生产时套管产气量很小，认为环空流体为静止状态，忽略环空流体对流传热影响；
(6)以井口为坐标原点，沿井筒向下为坐标轴的正方向；
(a)确定井筒内产液温度沿井深的变化规律：
井筒内流体温度分布规律由方程(I)描述，
采用常数变易法求解，即：
式(I)、(II)中，T为井筒内某点处的产液温度，℃；Tc为井口产出液温度，℃；Ts为地表恒温层温度，℃；kl为传热系数，W/(m·K)；G为流体质量流量，kg/s；cp为流体比定压热容，J/(kg·K)；x为以井口为原点，所述计算点的深度，m；m为地温梯度，℃/m；qv为内热源强度，W/m；g为重力加速度，m/s2；(b)根据潜油电泵井井筒结构和动液面位置，将井筒分为四段：即井口至动液面段、动液面至泵出口段、泵机组段和电机至井底段，各段的温度计算模型为：
(1)井口至动液面段温度计算模型为
其中，
kl1为井口至动液面段内的井筒与地层之间的传热系数，W/(m·K)；hl为产液与油管内表面的表面传热系数，W/(m·K)；dti为油管内径，m；λtub为油管的导热系数，W/(m·K)；dto为油管外径，m；λr为油管外表面与套管内表面间辐射传热过程的当量导热系数, W/(m·K),由于油管外表温度和套管内表面温度未知，为计算传热系数kl，需采用迭代法求解；λcas为套管的导热系数，W/(m·K)；dco为套管外径，m；dci为套管内径，m；λcem为水泥环的导热系数，W/(m·K)；dh为水泥环外缘直径，m；Rle为单位长度地层热阻，(m·K)/W；σb为Stefan-Boltzmann常数，取值为5.763×10-8W/(m2·K4)；Tto为油管外表面温度，K；Tci为油管外表面温度，K；εc，εt分别为套管内表面和油管外表面的发射率，小数；I为泵机组工作时的工作电流，A；R0为大扁电缆单位长度的电阻值，Ω/m；
(2)动液面至泵出口段温度计算模型为
其中，
式中，kl2为动液面至泵出口段内的井筒与地层之间的传热系数，W/(m·K)；λl为流体导热系数，W/(m·K)；Td为动液面处的井筒流体温度，℃；Hd为动液面深度，m；
泵出口处的井筒流体温度为：
式中，Tpout为泵排出口处流体的温度，℃；Hp为下泵深度，m；
(3)泵机组段温度计算模型为
Tpint＝Tpout-ΔTdj-ΔTdl (V)
式中，Tpint为泵吸入口处流体的温度，℃；ΔTdj为电机增温，℃；ΔTdl为泵内小扁电缆的发热增温，℃；
(4)电机至井底段温度计算模型为
其中，
式中，kl3为电机至井底段内的井筒与地层之间的传热系数，W/(m·K)；Ls为泵内小扁电缆长度，m；
所述步骤二：建立基于地面电参数的潜油电泵所需输入功率计算模型：
式(VII)中，Np为潜油电泵输入功率，kW；U为地面电压，V；I为地面电流，A；cosφ为功率因数，小数；Rc为小扁电缆电阻，Ω/m；Lp为下放潜油电泵所需的电缆长度，m；Ls为电泵内小扁电缆的长度m；ηm为潜油电机工作时的泵效，小数；P分离器和P保护器分别为分离器和保护器所消耗的功率，kW；
所述步骤三：建立潜油电泵扬程计算模型：
式(VIII)中，H为潜油电泵有效举升高度，m；ηo为潜油电泵实际泵效，小数；Q为泵的实际排量，m3/d；γl为井液平均相对密度，小数；
根据潜油电泵工作特性曲线，其实际泵效由式(IX)表示：
η0＝kη[anQn+…+a2Q2+a1Q+a0] (IX)
式(IX)中，kη为泵效的修正系数，小数；an，···，a2，a1，a0为潜油电泵特性曲线泵效与排量曲线回归系数，常数；
所述步骤四：建立泵吸入口压力模型：
以泵为节点，计算出泵吸入口的压力，即沉没压力，
pin＝pout-ρlgH (X)
式(X)中，pin为泵吸入口处的压力，MPa；pout为泵排出口处的压力，MPa；ρl为泵内流体的混合密度，kg/m3；
所述步骤五：潜油电泵井生产时，井口到管鞋的油套环空中流体分为两段，即气柱段和油柱段，以Δh为计算单元步长，将井筒深度划分为个网格，计算每个网格的压力，以提高模型求解精度；由气柱段压力计算模型和液面以下油柱段压力计算模型联合求解，即确定动液面深度：
式(XI)中，pg(i)为第i网格深度处的气柱压力，MPa；pg(i-1)为第i-1网格深度处的气柱压力，MPa；ρg0为标准状况下的气体密度，kg/m3；T0为标准状况下的温度，K；p0为标准状况下的压力，MPa；Δh为网格划分单元格的精度，m；i为网格划分的网格数，整数；Tav为第i-1和第i深度该段Δh内的平均温度，K；Zav为第i-1和第i深度该段Δh内的平均温度和平均压力下的气体压缩因子，小数；po(x)为x深度处的油柱压力，MPa；ρo(x)为微元段dx内平均压力和平均温度下的油的密度，kg/m3；pc为井口套管压力，MPa；
计算流程如下：
①以计算单元步长Δh将井筒网格划分，假设动液面Hd，计算井筒流体温度和压力分布，得到泵排出口压力pout；
②利用公式(VII)，通过地面实测电参数，计算潜油电泵所需的轴功率Np；
③利用公式(VIII)，根据潜油电泵的工作特性，计算有效举升高度H；
④利用公式(X)，以泵为节点，得到泵吸入口压力pin；
⑤利用公式(XI)所述的气柱段压力计算模型和液面以下油柱段压力计算模型联合求解动液面Hd’；
⑥比较Hd’与Hd，若|Hd'-Hd|<ε，输出结果Hd’；若不满足上述条件，将令Hds＝0.618Hd’+0.382Hd，将Hds作为新的动液面假设值，重复步骤①～⑤计算新的动液面计算值Hd’，直至满足|Hd'-Hds|<ε要求，输出新的动液面计算值Hd’。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611194201.5A CN106777709B (zh) | 2016-12-21 | 2016-12-21 | 基于地面实测电参数的潜油电泵井动液面计算的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611194201.5A CN106777709B (zh) | 2016-12-21 | 2016-12-21 | 基于地面实测电参数的潜油电泵井动液面计算的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777709A CN106777709A (zh) | 2017-05-31 |
| CN106777709B true CN106777709B (zh) | 2018-11-06 |
Family
ID=58899100
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611194201.5A Expired - Fee Related CN106777709B (zh) | 2016-12-21 | 2016-12-21 | 基于地面实测电参数的潜油电泵井动液面计算的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777709B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107387030B (zh) * | 2017-07-26 | 2020-04-10 | 中国海洋石油总公司 | 一种利用本井气的接力举升采油设计方法 |
| CN107729631B (zh) * | 2017-09-28 | 2021-06-01 | 中国石油天然气股份有限公司 | 一种利用液面测试数据进行井底压力计算的方法及装置 |
| CN107630688B (zh) * | 2017-11-09 | 2020-06-02 | 中国石油大港油田勘探开发研究院 | 一种高凝油油藏地层掺稀冷采新技术评价方法 |
| CN107939378B (zh) * | 2017-11-30 | 2021-04-20 | 中国石油大学(华东) | 一种地面驱动螺杆泵井动液面实时获取方法 |
| CN108266176B (zh) * | 2018-01-08 | 2021-06-01 | 中国海洋石油集团有限公司 | 一种基于井筒模型的天然气井口流量计算方法 |
| CN108798612B (zh) * | 2018-05-16 | 2020-11-03 | 辽宁瑞邦石油技术发展有限公司 | 一种无杆油井举升系统的智能控制方法 |
| CN109740661B (zh) * | 2018-12-28 | 2022-10-28 | 沈阳工业大学 | 一种基于模糊评估的油井动液面自适应预测方法 |
| CN110284875A (zh) * | 2019-06-13 | 2019-09-27 | 中国石油化工股份有限公司 | 一种基于电流参数的电潜泵井工况诊断方法和装置 |
| CN114412439A (zh) * | 2021-12-31 | 2022-04-29 | 沈阳中科奥维科技股份有限公司 | 一种螺杆泵井的产量预测方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105257279A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种抽油机井动液面的测量方法 |
| CN105298471A (zh) * | 2014-12-04 | 2016-02-03 | 大庆高新区中环电力控制系统有限公司 | 油井动液面连续远程监控系统 |
-
2016
- 2016-12-21 CN CN201611194201.5A patent/CN106777709B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105298471A (zh) * | 2014-12-04 | 2016-02-03 | 大庆高新区中环电力控制系统有限公司 | 油井动液面连续远程监控系统 |
| CN105257279A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油天然气股份有限公司 | 一种抽油机井动液面的测量方法 |
Non-Patent Citations (3)
| Title |
|---|
| 基于功图的油井动液面计算模型及应用;陈德春等;《科学技术与工程》;20151130;第15卷(第32期);32-35 * |
| 潜油电泵井油套环空泵上掺稀油井筒流体温度模型;陈德春等;《新疆石油地质》;20131031;第34卷(第5期);591-594 * |
| 电潜泵井憋压法计算动液面技术及其应用;杜林辉等;《石油钻采工艺》;20110531;第33卷(第3期);116-118 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777709A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777709B (zh) | 基于地面实测电参数的潜油电泵井动液面计算的方法 | |
| Ma et al. | A coupled heat transfer model of medium-depth downhole coaxial heat exchanger based on the piecewise analytical solution | |
| CN111400893B (zh) | 一种套管式地埋管换热器流体温度场分析方法 | |
| CN109086560B (zh) | 变工况下地源热泵竖直单u型地埋管流体温度分布预测方法 | |
| CN106968667A (zh) | 一种温度场预测方法及装置 | |
| CN105422084B (zh) | 一种热水循环加热降粘工艺的井筒温度场获取方法及装置 | |
| CN107939378B (zh) | 一种地面驱动螺杆泵井动液面实时获取方法 | |
| CN103726815B (zh) | 一种co2驱采出井井筒流态确定及参数优化方法 | |
| CN115293066B (zh) | 考虑地层渗流传热效应的气井温度场计算方法 | |
| CN104389589A (zh) | 基于空心杆确定井筒温度场分布的方法以及系统 | |
| CN105888646B (zh) | 毛细管测压电泵井在线流量实时计量系统与方法 | |
| CN114048657B (zh) | 一种水平井钻井井筒温度场分布模型的构建方法 | |
| CN109033012B (zh) | 一种空心抽油杆注热水循环井筒温度场的确定方法及装置 | |
| Dong et al. | Model calculations and factors affecting wellbore temperatures during SRV fracturing | |
| CN106934106B (zh) | 一种获取固井循环温度的方法及装置 | |
| CN115629102B (zh) | 一种基于地质分层的岩土热物性参数辨识方法 | |
| TWI734029B (zh) | 地熱井產能預估方法 | |
| CN110750918A (zh) | 一种二氧化碳压裂过程中井筒温度的预测方法 | |
| CN110566164A (zh) | 一种海上测试管柱最优产量的确认方法 | |
| CN113515863B (zh) | 一种基于数值反演的中深层套管式换热器取热量计算方法 | |
| CN105114060B (zh) | 一种计算注气井的井筒温度分布的方法及装置 | |
| CN109339765B (zh) | 一种海上油田油井动态清蜡周期预测方法 | |
| RU2677313C1 (ru) | Способ эксплуатации нефтяной скважины установкой электроцентробежного насоса | |
| Zhang et al. | Supercritical cyclic stream simulation of wellbore temperature and pressure distribution in Lukeqin oilfield | |
| CN113505454B (zh) | 一种中深层地热井套管式换热器取热量计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181106 |