Simulation of Stator and Rotor Meshing Rotation of Double Helix Single Screw Pump
-
摘要: 为了提高双头单螺杆泵的工作效率和使用寿命,研究了结构参数对其工作效率及使用寿命的影响。基于橡胶Mooney-Rivlin双参数模型理论,考虑橡胶衬套材料具有非线性变形特性、内压及转子惯性对橡胶衬套的作用,建立了双头单螺杆泵一个导程的定转子三维啮合模型,分析了过盈量、偏心距对螺杆泵定子衬套应力和摩擦阻力矩的影响。数值模拟结果显示,偏心距选择7.5 mm时,过盈量由0.2 mm增加到0.6 mm,螺杆泵摩擦阻力矩从88.52 N·m增大到466.51 N·m,衬套应力峰值由1.4 MPa增加到1.8 MPa;而过盈量取0.4 mm时,偏心距由7.3 mm增大到7.5 mm,摩擦阻力矩从240.00 N·m仅增加到264.77 N·m,衬套应力峰值由1.28 MPa增加到1.62 MPa。这表明,合理搭配过盈量和偏心距可以提高螺杆泵的工作效率和使用寿命,为螺杆泵设计及结构参数优化提供了理论依据。Abstract: In order to improve working efficiency and service time of double helix single screw pump,the effects of its structure parameters on working efficiency and service time were studied.Based on double parameter Mooney-Rivlin model,a three-dimensional meshing model taking material nonlinearity,internal pressure and inertia effect of the rotor on the rubber lining into account,was established for the double helix single screw pump.The effects of interference fit and eccentricity change on the stator rubber lining’s stress and fatigue failure were acquired from the analysis.Simulation results showed that frictional torque increased from 88.52 N·m to 466.51 N·m and stator lining stress increased from 1.4 MPa to 1.8 MPa,when interference fit increased from 0.2 mm to 0.6 mm under eccentricity 7.5 mm.Frictional torque increased only from 240.00 N·m to 264.77 N·m,when eccentricity increased from 7.3 mm to 7.5 mm under interference fit 0.4 mm,however,the peak stress on stator lining increased from 1.28 MPa to 1.62 MPa.This showed that the rational combination of interference fit and eccentricity could increase work efficiency and service time of the screw pump.The work in this paper established a theory foundation for the design and structure parameter optimization of the double helix single screw pump.
-
Keywords:
- screw pump /
- stator /
- rubber /
- lining /
- rotor /
- geometrical model /
- mathematic model
-
页岩油气储层致密,能量保持难、流动阻力大,大规模压裂改造是其高效开发的主要技术手段[1-5]。东营凹陷陆相页岩油资源量丰富,亟需采用新技术实现高效开发[6-7]。超临界CO2具有低黏度、高渗透、低表面张力的物理特性,被认为是非常规储层改造新方法之一[8-10]。然而,超临界CO2压裂的相态和物理性质复杂,造缝机理尚不清楚,学者们多采用试验和数值模拟方法开展相关研究[11-15]。试验表明,CO2可降低破裂压力[16],超临界CO2压裂形成的裂缝表面相对于滑溜水压裂产生的裂缝表面更为粗糙[17],使用超临界CO2在降低岩石破裂压力和造缝方面具有很大的优势[18]。目前,相关研究多针对CO2低黏度特性对压裂的影响,对CO2对岩石物性影响的关注较少。此外,CO2压裂裂缝形态复杂,传统方法模拟裂缝多为单缝结构,采用现的数值模拟手段难以精确刻画前置CO2压裂缝网[19-21]。
赵辉等人[22-24]提出了基于闪电模拟的油藏压裂模型,可以很好地表征裂缝复杂性,但仅适用于水力压裂,未考虑压裂过程中CO2对裂缝扩展形态的影响。借鉴该方法,笔者考虑CO2的相态变化特征及CO2对岩石性质的影响,结合物理试验和数值模拟,提出了基于节点连接方法(NCM)的CO2增能压裂裂缝扩展模拟模型,在此基础上进一步讨论了前置CO2比例、地质参数、物性参数及压裂施工参数对裂缝形态的影响,并在东营凹陷陆相页岩油X井开展分析应用。
1. CO2增能压裂裂缝扩展模拟模型
1.1 CO2岩石压裂起裂压力与应力变化表征
为准确表征胜利油田东营凹陷陆相页岩经过CO2浸泡后的岩石力学与破裂特征,选取2口典型井12块岩样,利用岩石力学仪测试其不同CO2浸泡时间下和不同围压下的强度、弹性模量、泊松比、内聚力等岩石力学参数,结合破裂压力剪切压力等计算公式,发现东营凹陷陆相页岩油岩体破裂压力、天然裂缝剪切及张开应力和地应力随CO2浸泡时间增长呈现线性降低的趋势(见图1)。
1.2 CO2物性参数变化表征
前置CO2注入过程中,随着温度和压力条件变化,CO2的密度和黏度均会发生显著变化,影响压裂裂缝扩展。采用Pen-Robinson方程刻画CO2密度变化,其状态方程可表示为[25]:
p=8.314TV−26.667−396 306.77[1+0.707 979(1−√Tr)]2V(V+26.667)+26.667(V−26.667) (1) 式中:T为CO2的绝对温度,K;V为CO2的摩尔体积,L/mol;Tr为对比温度,即CO2绝对温度T与临界温度Tc(Tc=304.13K)之比。
根据式(1)确定CO2的体积后,即可进一步获得CO2在该温度和压力条件下的密度。
CO2黏度与温度的关系[21]为:
μ=36.344[√1.259 3TrΩuFc(1G2+E6cVe6)+η**]√MTc3√V2c (2) 其中,部分系数定义为:
Fc=1−0.275 6ω+0.059 035μ4r+κ (3) η**=E7y2G2exp(E8+E9T∗+E10T∗2) (4) G2=E1[1−exp(−E4y)]y+E21−0.5y(1−y)3exp(E5y)−E31−0.5y(1−y)3E1E4+E2+E3 (5) Ei=fi(ω,κ) (6) 式中:Tr为临界温度,K;Fc为形状和极性因子;Ve为气体体积,m3;ω为偏心因子;y为摩尔分数;T*为绝对温度,K;M为摩尔质量,g/mol;Vc为CO2的临界体积,m3;Ωu为碰撞积分;μr为约化偶极矩;κ为高极性物质的关联因子;c为单位体积物质的量,mol/m3。
1.3 基于NCM的CO2增能压裂裂缝扩展数学模型
基于节点体系扩展模拟方法(Node Connection Method, NCM)[22-23],考虑前置CO2注入后对岩石力学、物性参数的影响,构建前置CO2增能压裂裂缝扩展模拟模型。
1.3.1 前置CO2水力压裂力学机制
假设岩石为多孔弹性介质,根据测井资料获取地层物性参数,以三维弹性理论为基础,可以得到地下水平主应力的计算方法。结合试验研究结果,储层中岩石的泊松比、弹性模量等参数会随着CO2浸泡时间增长而变化,在模型中需要考虑。因此,水平主应力与时间的关系可以表示为:
{σth=νt1−νtσv−νt1−νtαvpp+αhpp+Et1−(νt)2ξh+νtEt1−(νt)2ξHσtH=νt1−νtσv−νt1−νtαvpp+αHpp+Et1−(νt)2ξH+νEt1−(νt)2ξh (7) 式中:σth为t时刻的最小水平主应力,MPa;σtH为t时刻的最大水平主应力,MPa;σv为垂向应力,MPa;αv、αh为垂向和水平的Biot系数;ξH、ξh为最大水平主应力和最小水平主应力的应变系数;νt为t时刻泊松比;Et为t时刻的弹性模量,GPa;pp为孔隙压力,MPa。
压裂过程中,新生缝会改变岩石的力学状态,影响后续裂缝扩展,简称应力阴影效应。采用诱导应力解析解和叠加原理描述全局应力状态[26]:
{σt+1xx(i,j)=σtxx(i,j)+N∑m=1pn,m(Lm√L1,mL2,mcos2θm−θ1,m−θ2,m2−2)−N∑m=1pn,mLma(a2L1,mL2,m)3/2sinθmsin32(θ1,m+θ2,m)σt+1yy(i,j)=σtyy(i,j)+N∑m=1pn,m(Lm√L1,mL2,mcos2θm−θ1,m−θ2,m2−2)+N∑m=1pn,mLma(a2L1,mL2,m)3/2sinθmsin32(θ1,m+θ2,m)τt+1xy(i,j)=τtxy(i,j)+N∑m=1pnet,mLma(a2L1,mL2,m)3/2⋅sinθmcos32(θ1,m+θ2,m) (8) 式中:σt+1xx和σtxx分别为t+1和t时间步下(i, j)位置处x轴向主应力,MPa;σt+1yy和σtyy分别为t+1和t时间步下(i, j)位置处y轴向主应力,MPa;τt+1xy和τtxy分别为t+1和t时间步下(i, j)位置处剪切应力,MPa;N为裂缝微元数量;pn,m为裂缝m微元的缝内净压力,MPa;a为裂缝微元半长,m;Lm,L1,m和L2,m为位置(i, j)到裂缝m微元的距离,m;θm,θ1,m和θ2,m为位置(i, j)到裂缝m微元的角度。
对裂缝内流体流动简化处理,不考虑裂缝在缝宽方向上的流动过程,且裂缝壁面上的流速为0,则流体在裂缝内为一维流动。对流体压降方程进行处理,可得到缝内压力分布为:
pti,m=−12μtmqmxmhf,mw3m+p0 (9) 式中:pti,m为t时间步时裂缝m微元的净压力,MPa;μtm为t时间步时裂缝m微元中流体的黏度,mPa·s;qm为裂缝m微元的流量;xm为裂缝m微元与井筒的距离,m;hf,m为裂缝m微元的高度,m;p0为初始压力,MPa。
对于前置CO2压裂,裂缝内流体黏度随着CO2性质而改变,采用式(2)修正不同时间步下前置CO2的黏度。
将临界起裂应力作为起裂判据,该值与岩石断裂韧性有关,则裂缝扩展条件为[23]:
σfr=1√2πrcosθ2(K1cos2θ2−32K2sinθ)−σcr⩾ (10) 式中:σfr为残余起裂应力,MPa;σcr为临界起裂应力,MPa;K1为Ⅰ型应力强度因子,MPa·m0.5;K2为Ⅱ型应力强度因子,MPa·m0.5。
分支裂缝密度表示为在裂缝的改造范围内共有的裂缝节点数量,其计算公式为:
{\rho }_{\text{if}}=\frac{{n}}{{L}_{\text{f}}{w}_{\text{f}}} (11) 式中:ρif为分支缝密度,条/m2;n为裂缝节点数量;Lf为裂缝长度,m;wf为裂缝的带宽,m。
1.3.2 前置CO2压裂模拟流程
基于NCM的CO2增能压裂裂缝扩展模型计算流程如图2所示。具体步骤为:1)依据实际储层大小划分地质单元,作为计算的离散体;2)结合地质资料给地质单元赋值;3)根据实际压裂规模,设置射孔参数;4)计算现存裂缝单元缝内流体压力分布;5)结合应力阴影效应及缝内流体流动特征,校正储层应力分布;6)计算裂缝尖端周向应力和临界起裂应力;7)判断是否存在满足起裂条件的裂缝单元,不满足时输出裂缝形态;8)引入随机函数及概率分布,确定新的裂缝单元;9)根据CO2物性的变化,更新缝内流体物性参数;10)重复步骤4)—9),直至输出裂缝形态。
2. 前置CO2增能压裂裂缝扩展模拟
选用胜利油田东营凹陷陆相页岩油储层地质力学参数以及物性参数,基于构建模型,分析不同影响因素下裂缝扩展形态。模型尺寸大小为160 m×300 m,1段射孔3簇,射孔点在模型中间位置处,模拟裂缝排量为16 m3/min,页岩基质的弹性模量为26.36 GPa,泊松比为0.30。
2.1 前置CO2比例对缝网形态的影响
模拟前置CO2比例分别为0.1、0.2和0.3时的裂缝形态(见图3),分析不同前置CO2比例与裂缝长度缝网带宽及分支缝密度的关系(见图4)。
前置CO2比例对整体缝网扩展长度影响最大,前置CO2比列由0.1增至0.3时,缝网长度缩短近40%。压裂缝网带宽主要受压裂簇数影响,前置CO2比例对整体缝网扩展带宽影响较小。前置CO2比例对分支缝密度的影响很大,当前置CO2比例由0.1增加至0.3时,分支缝网密度由0.06 条/m2增加至0.13 条/m2,增加了117%。这主要是因为超临界CO2易于进入基质孔隙和微裂隙中,促进了基质的破裂和微裂隙的开启,增加了分支裂缝的数量,形成了复杂的网状裂缝形态。
2.2 水平应力差对裂缝形态的影响
模拟水平应力差分别为0,15和30 MPa时的裂缝形态(见图5),分析不同水平应力差与裂缝长度缝网带宽及分支缝密度的关系(见图6)。
应力差与缝长相关性很大,与缝宽相关性较小,与分支缝密度相关性大,应力差主要影响缝长及分支缝密度。水平应力差由5 MPa增至30 MPa时,裂缝长度由175 m增至266 m,缝网带宽由63 m降至58 m,降低了3.18%,分支缝密度由0.095 条/m2降至0.074 条/m2,下降了31.85%。相较于滑溜水压裂,前置CO2增能压裂在高应力差下依然可以形成复杂裂缝网。因此,高应力差储层采用前置CO2进行压裂,可以提高压裂改造效果,实现高效增产。
2.3 储层渗透率对裂缝形态的影响
模拟储层渗透率分别为0.05,0.5和5 mD时的裂缝形态(见图7),分储层渗透率与裂缝长度、缝网带宽及分支缝密度的关系(见图8)。结果发现,储层渗透率对裂缝长度和带宽的影响较小,对分支缝密度的影响较大。总体而言,储层渗透率影响压裂流体的滤失性,前置CO2压裂对于储层渗透率的敏感程度要远大于滑溜水压裂。其原因在于,储层渗透率影响CO2向基质的滤失量和波及范围,对于开启天然裂缝和储层增压作用明显。因此,对于高渗透储层,应适当增大排量和压裂液量,从而增大压裂改造范围。
3. 实例分析
东营凹陷陆相页岩油X井目的层为沙四上纯上3层组,水平段长度1 976 m,钻遇岩相以泥质灰页岩(970 m)、灰质泥页岩(550 m)为主;前置CO2增能压裂完成34段、120簇压裂施工,平均3.5 簇/段,加砂4 300.6 m3(设计4 204.7 m3)、注入携砂液73 166 m3,平均单段加砂126.5 m3、加砂强度2.4 m3/m。结合实际施工参数模拟裂缝扩展形态,X井36段反演裂缝平均缝长为277.80 m,平均缝宽为59.85 m。X井压裂时的微地震监测结果见图9,裂缝扩展模拟结果见图10。
X 井压裂过程中进行单井尺度模拟时,由于模拟布点精度及显示尺度的原因,裂缝整体上呈双翼对称复杂裂缝形态。对比微地震监测结果发现,水力裂缝缝长在微地震的75.09%~99.63%,水力裂缝缝宽在微地震的78.79%~98.84%,模拟结果与实际监测具有较好的一致性。此方法可用于压裂前裂缝形态模拟,根据模拟结果制定压裂方案,提高压裂效率。
4. 结 论
1)综合考虑CO2对岩体破裂压力、地应力的影响以及CO2物性参数的变化,基于节点连接方法的CO2增能压裂裂缝扩展数学模型可实现CO2增能压裂复杂裂缝网络形态的模拟。
2)前置CO2比例、应力差和渗透率对分支缝密度的影响较大,在高应力差储层增加缝网密度时,可适当增加前置CO2比例。
3)实际区块模拟结果与微地震监测结果吻合度较高,具有一定的可靠性,基于节点连接方法的CO2增能压裂裂缝扩展数学模型具有一定的现场应用潜力,可以为压裂方案的制定提供理论指导。
-
[1] 申小敏, 王玉堂, 王永站, 等.特高含水期油藏开发对策研究[J].内蒙古石油化工, 2010, 36(21):152-153. Shen Xiaomin, Wang Yutang, Wang Yongzhan, et al.Exploit countermeasures study of reservoirs under super-high water cut stage[J].Inner Mongolia Petrochemical Industry, 2010, 36(21):152-153. [2] 石志敏, 刘国生, 毛为成, 等.双河油田特高含水期油藏深度调剖技术试验研究[J].石油钻探技术, 2005, 33(4):63-65. Shi Zhimin, Liu Guosheng, Mao Weicheng, et al.Lab test on depth profile control technology for reservoirs under super-hish water cut stage in the Shuanghe Oilfield[J].Petroleum Drilling Techniques, 2005, 33(4):63-65. [3] 孟令尊, 许军.地面驱动大排量双头螺杆泵试验研究[J].石油机械, 2004, 32(7):37-38. Meng Lingzun, Xu Jun.Test research of the ground drive large displacement double screw pump[J].China Petroleum Machinery, 2004, 32(7):37-38. [4] 金红杰, 吴恒安, 王秀喜.螺杆泵腔室液压分布数值模拟研究[J].机床与液压, 2010, 38(1):92-94. Jin Hongjie, Wu Heng’an, Wang Xiuxi.Numerical simulation on cavity pressure of progressive cavity pump[J].Machine Tool Hydraulics, 2010, 38(1):92-94. [5] 叶卫东, 郭玉双, 杜秀华, 等.螺杆泵内部压力分布规律研究[J].科学技术与工程, 2009, 9(11):3069-3072. Ye Weidong, Guo Yushuang, Du Xiuhua, et al.Study on pressure law inside of the progressive cavity pump[J].Science Technology and Engineering, 2009, 9(11):3069-3072. [6] 李晓芳, 孙春龙, 李强.螺杆泵定转子过盈接触三维有限元分析[J].润滑与密封, 2008, 33(10):64-65. Li Xiaofang, Sun Chunlong, Li Qiang.3D-FEM analysis of stator-rotor interference fit contact of progressing cavity pump[J].Lubrication Engineering, 2008, 33(10):64-65. [7] 陈玉祥, 王霞, 周松, 等.提高螺杆泵定子橡胶材料寿命的分析与研究[J].排灌机械, 2005, 23(4):6-9. Chen Yuxiang, Wang Xia, Zhou Song, et al.Analyses on improving the life of the screw pump’s stator rubber[J].Drainage and Irrigation Machinery, 2005, 23(4):6-9. [8] 张劲, 张士诚, 帅国臣.常规螺杆泵定子有限元分析[J].力学季刊, 2003, 24(4):590-598. Zhang Jin, Zhang Shicheng, Shi Guochen.FEM analysis for a general screw pump stator[J].Chinese Quarterly of Mechanics, 2003, 24(4):590-598. [9] 陈次昌, 魏存祥, 单代伟, 等.双头单螺杆泵衬套有限元分析[J].西南石油大学学报, 2007, 29(3):1-4. Chen Cichang, Wei Cunxiang, Shan Daiwei, et al.Finite element analysis for a double-helix single screw pump bush[J].Journal of Southwest Petroleum University, 2007, 29(3): 1-4. [10] 祝效华, 石昌帅, 童华.静压、压差与泊松比对定子衬套变形的影响[J].西南石油大学学报:自然科学版, 2010, 32(6):175-179. Zhu Xiaohua, Shi Changshuai, Tong Hua.Influence of static press, pressure difference and Poisson’s ratio on deformation of stator lining[J].Journal of Southwest Petroleum University:Science Technology Edition, 2010, 32(6):175-179. [11] 韩传军, 邱亚玲, 刘清友, 等.螺杆钻具等壁厚衬套性能分析[J].西南石油大学学报:自然科学版, 2008, 30(4):163-165. Han Chuanjun, Qiu Yaling, Liu Qingyou, et al.Perormance analysis for uniform wall thickness bushing of PDM[J].Journal of Southwest Petroleum University:Science Technology Edition, 2008, 30(4): 163-165. [12] 王伟, 邓涛, 赵树高.橡胶Mooney-Rivlin模型中材料常数的确定[J].特种橡胶制品, 2004, 25(4):8-10. Wang Wei, Deng Tao, Zhao Shugao.Determination for material constants of rubber mooney-rivlin model[J].Special Purpose Rubber Products, 2004, 25(4):8-10. [13] Korunovic N, Trajanovic M, Stojkovic M.FEA of tyres subjected to static loading[J].Journal of Serbian Society for Computational Mechanics, 2007, 1(1):87-98.
[14] Zhu Xiaohua, Liu Shaohu, Shi Changshuai, et al.Successful application of the new double-effect protector in deep well drilling[R].SPE 130924, 2010.
[15] 葛占玉, 万邦烈.单螺杆式水力机械螺杆-衬套副摩擦阻力矩的计算[J].石油大学学报:自然科学版, 1992, 16(3):46-52. Ge Zhanyu, Wan Banglie.Calculation of friction resistance moment for screw-liner pair on single-screwed hydraul machinery[J].Journal of the University of Petroleum, China:Edtion of Natural Science, 1992, 16(3):46-52. -
期刊类型引用(9)
1. 徐智颢,杜博,李东东. 裂缝内CO_2压裂液相态分布特征的数值模拟研究. 中国资源综合利用. 2024(04): 41-43 . 百度学术
2. 孟宪波. 页岩油藏水平井分段多簇压裂全井段裂缝扩展数值模拟——以济阳坳陷沙河街组为例. 断块油气田. 2024(03): 459-466 . 百度学术
3. 李阳,祝仰文,李宗阳,蒋廷学,薛兆杰,沈子齐,肖朴夫,于洪敏,程紫燕,赵清民,张庆福. 济阳坳陷页岩油注二氧化碳开发技术. 石油勘探与开发. 2024(04): 855-864 . 百度学术
4. 孙丙向,成海,毕研涛. 石油工程“四提”技术进展与攻关方向. 石油科技论坛. 2024(03): 85-94 . 百度学术
5. 李阳,曹小朋,赵清民,刘祖鹏,薛兆杰,蒋龙. 济阳坳陷陆相断陷盆地页岩油开发的几点思考. 石油钻探技术. 2024(04): 1-7 . 本站查看
6. 吴壮坤,张宏录,池宇璇. 苏北页岩油二氧化碳强压质换技术. 石油钻探技术. 2024(04): 87-93 . 本站查看
7. LI Yang,ZHU Yangwen,LI Zongyang,JIANG Tingxue,XUE Zhaojie,SHEN Ziqi,XIAO Pufu,YU Hongmin,CHENG Ziyan,ZHAO Qingmin,ZHANG Qingfu. Shale oil recovery by CO_2 injection in Jiyang Depression, Bohai Bay Basin, East China. Petroleum Exploration and Development. 2024(04): 981-992 . 必应学术
8. 余前港,荣双,毛国扬,崔连文,彭根博,向传刚,王旭. 阳春沟区块页岩气超临界CO_2增能压裂研究与应用. 石油与天然气化工. 2024(05): 77-83 . 百度学术
9. 李阳,赵清民,杨勇,张峰,曹功泽,张潦源,杨峰,鲁明晶. 胜利油田陆相页岩油前置CO_2压裂技术研究与实践. 断块油气田. 2024(06): 945-954 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 2533
- HTML全文浏览量: 54
- PDF下载量: 3751
- 被引次数: 9