Processing math: 100%

层理页岩水平井井周剪切失稳区域预测方法

马天寿, 陈平

马天寿, 陈平. 层理页岩水平井井周剪切失稳区域预测方法[J]. 石油钻探技术, 2014, 42(5): 26-36. DOI: 10.11911/syztjs.201405005
引用本文: 马天寿, 陈平. 层理页岩水平井井周剪切失稳区域预测方法[J]. 石油钻探技术, 2014, 42(5): 26-36. DOI: 10.11911/syztjs.201405005
Ma Tianshou, Chen Ping. Prediction Method of Shear Instability Region around the Borehole for Horizontal Wells in Bedding Shale[J]. Petroleum Drilling Techniques, 2014, 42(5): 26-36. DOI: 10.11911/syztjs.201405005
Citation: Ma Tianshou, Chen Ping. Prediction Method of Shear Instability Region around the Borehole for Horizontal Wells in Bedding Shale[J]. Petroleum Drilling Techniques, 2014, 42(5): 26-36. DOI: 10.11911/syztjs.201405005

层理页岩水平井井周剪切失稳区域预测方法

基金项目: 

国家重点基础研究发展计划("973"计划)课题"页岩气水平井钻完井关键基础研究"(编号:2013CB228003)资助。

详细信息
    作者简介:

    马天寿(1987—),男,四川绵阳人,2009年毕业于西南石油大学石油工程专业,2012年获西南石油大学油气井工程专业硕士学位,在读博士研究生,主要从事页岩气钻井、深水钻井与随钻测量技术方面的研究。

  • 中图分类号: TE21

Prediction Method of Shear Instability Region around the Borehole for Horizontal Wells in Bedding Shale

  • 摘要: 为了揭示层理页岩强度各向异性对水平井井周剪切失稳的影响规律,基于直剪试验研究了川南龙马溪组页岩的强度各向异性特征,并与国外典型层理页岩进行了对比;建立了层理页岩井周坐标转换关系模型,结合井周应力弹性解、Mohr-Coulomb弱面强度理论、Mohr-Coulomb准则,建立了层理页岩地层井周剪切失稳区域预测方法;分析了钻进方位对井周剪切失稳区域分布规律的影响。结果表明:弱面与本体强度参数基本呈线性关系,本体内聚力约为弱面内聚力的2.10倍,本体内摩擦角约为弱面内摩擦角的1.20倍,其中川南龙马溪组页岩强度各向异性特征更加显著;考虑层理等弱面影响计算出的失稳区域的形状不是传统的"狗耳朵"形状,而是近似于长方形,这种形状与典型的物模试验结果基本吻合;若页岩地层中存在显著的层理,则井周地层更容易沿弱面剪切滑动,导致水平井坍塌失稳的风险显著增大。利用该方法可以分别计算出井周地层沿弱面剪切、本体剪切和全部剪切的失稳区域分布,为层理发育地层井眼稳定分析、成像测井反演地应力方位等研究提供基础理论支撑。
    Abstract: In order to understand the effects of strength anisotropy of bedding shale on shear instability region around the borehole of a horizontal well,the strength anisotropy characteristics for Longmaxi shale in the southern of Sichuan were researched based on direct shear tests,and the test result comparison was made with typical bedding shale data of some other countries,so as to establish the coordinate transformation model around the borehole in the bedding shale.In combination with stress solution around a borehole,Mohr-Coulomb weak plane strength theory,Mohr-Coulomb criterion,the prediction method of shear instability region around the borehole of a horizontal well in bedding shale was established,and the influence of borehole azimuth on shear instability region was also analyzed.It showed that strength parameters of weak plane and rock mass are basically in a linear relationship as the cohesion and internal friction angle of rock mass are about 2.10 and 1.20 times of those of weak planes,respectively.Strength anisotropy was more evident in Longmaxi shale.The shear instability region around the borehole was not in the commonly seen "dog ear" shape due to the weak planes.Instead,it had a square shape,which agreed with typical physical experimental results.If bedding planes in shale reservoir are significantly developed,the formation around borehole may experience shear sliding along the weak planes,which would increase the risk of instability in horizontal wells.The prediction method in this paper could be used to calculate the shear instability regions from weak planes shearing and rock mass shearing providing the basis in borehole stability analysis of bedding formation and the in-situ stress orientation inversion from imaging logging.
  • 目前,全球已建成715座储气库。中国地下储气库建设起步于20世纪70年代,至今已建成25座储气库[1-2],其中呼图壁储气库是目前国内规模最大的储气库,总库容达107×108 m3。随着储气库日采气量增大,储气库调峰周期性的采–关(阀门由采气状态向关井状态调整)过程中会产生较大的波动压力[3-6],导致发生阀芯刺漏、油管柱密封失效等问题,会造成天然气窜漏或泄漏,可能引发火灾、爆炸等重大事故。挪威科技工业研究也证实,储气库采–关气井的泄漏情况是一般生产井的2~3倍[7-10]

    国内外研究者已开展了储气库储层应力敏感、气体速敏等实验,揭示了储气库注气压力动态变化诱发的储层微粒运移规律[11]。同时,采用数值模拟、有限元方法等方法,揭示了储气库井筒内的动力特性和频域变化特征,探究了储气库注采过程管柱非线性振动特性[12-13]。目前,大量相关研究集中在数值模拟分析方面,基于储气库采–关作业转换过程的波动压力变化规律尚不明确,且少有系统开展阀门、油管系统以及储气库采–关作业油管气相波动压力计算分析的研究[14-18]。为此,笔者考虑气相压缩性、气相压力波速[19]、开关阀门速度等因素,建立了呼图壁储气库采–关作业油管气相波动压力数学模型,提出了具体特征线计算方法,以呼图壁储气库H2井为例,分析了储气井日产气量、油管深度、油管内径、关阀时间、气体流量等因素对波动压力的影响规律,以期为现场阀门调节提供理论指导。

    根据动量守恒方程,建立油管气相波动压力运动方程:

    pA+ρg(A+Asds2)dssinβ+(p+psds2)Asds[pA+(pA)sds]τ0Xds=ρ(A+Asds2)dsdvdt (1)

    式中:p为压力,MPa;A为油管内横截面积,m2ρ为液相密度,kg/L;g为重力加速度,m/s2s为管柱长度,m;β为油管与地平面夹角,(°);τ0为湿相周切向力,N/m2X为管柱周长,m;v为气相流动速度,m/s;t为时间,s。

    通过化简式(1),约去微元项dsA等参数,得到:

    psρgsinβ+ρdvdt+τ0R=0 (2)

    考虑水力半径、油管深度、摩擦阻力等参数,根据连续守恒方程得到:

    zs+1ρgps+1g(vt+vvs)+fv|v|8Rg=0 (3)

    式中:R 为油管水力半径,m;f为油管摩阻系数。

    油管气相波动压力运动方程为:

    1ρps+vt+vvs+fv|v|8R=0 (4)

    根据质量守恒方程,可得到油管气相波动压力连续方程:

    [ρAvdt+(ρAvdt)sds]+(ρAds)tdt=ρAvdt (5)

    整理式(5),得到:

    (ρAv)s+(ρA)t=0 (6)

    将式(6)变形,可得:

    1ρ(vρs+ρt)+1A(vAs+At)+vs=0 (7)

    化简式(7),可得:

    (1ρdρdp+1AdAdp)dpdt+vt=0 (8)

    油管气相波动压力连续方程为:

    pt+vps+ρc2vs=0 (9)

    式中:c为压力波速,m/s。

    式(5)、式(6)的微分变量需满足:

    vtdt+vsds=dv (10)
    ptdt+psds=dp (11)

    根据特征线法求解原理,考虑式(10)与式(11)唯一解条件,得到如下行列式:

    |dsvdtρc2dt1ρdtdsvdt|=0 (12)

    对气相波动压力特征线第1族函数差分,可得:

    (pi,tpi1,tΔt)+ρcA(Qi,tQi1,tΔt)+ρcf8mΔtvi1,tΔt|vi1,tΔt|=0 (13)

    对气相波动压力特征线第2族函数差分,可得:

    (pi,tpi+1,tΔt)ρcA(Qi,tQi+1,tΔt)ρcf8mΔtvi + 1,tΔt|vi + 1,tΔt|=0 (14)

    式中:pi,ti节点在t时刻的压力,MPa;Qi,ti节点在t时刻的流量,m3/s;Δt为时间间隔,s;i为网格节点。

    储气库采–关工况条件下,油管气相波动压力求解步骤如下:

    1)将油管深度划分 n 个单元网格,当 t =0 时,波动压力为初始时间,此时井筒流体处于稳定流动状态,作为网格求解初始条件。波动压力传播过程中,将第 i 个( t 时刻)油管网格深度的气相波速、密度、黏度和压力视为恒定常数(网格内长度波动压力数值不发生变化),滑移计算不同网格的波动压力,得到沿着井筒不同网格变化的波动压力。

    2)确定每个网格初始稳态气相偏差因子、气体密度、气体黏度、油管压力和温度等物性参数,可通过气体状态方程等方法求取稳态油管气相压力。

    3)结合每个网格初始稳态气相物性参数,采用双流体数学模型,求解油管气相压力波速。

    4)输入t=0时刻(第1个网格)气相密度、温度、井口压力、气相黏度、油管内径等参数,迭代计算t时刻(第i个网格)的气相波动压力。

    实测波动压力数据与所建模型计算波动压力数据对比结果如图1所示。采用井筒波动压力测试结果[20],数学模型输入数据与测试条件一致。采用ϕ215.9 mm内径管体,测试井深3273.75 m、流体平均密度0.87 kg/L,流性指数0.725、测试时间40 s、井口压力0.101 MPa等关键参数,按照波动压力模型求解方法计算波动压力。与实测数据对比,最大误差≤5%,由于采–关作业油管与实验井筒所采用的波动压力模型及求解方法相同,因此采用实验井实测数据验证模型准确性具有一定可靠性。

    图  1  实测数据与本模型计算数据对比[20]
    Figure  1.  Comparison between measured data and calculated data by the model[20]

    呼图壁储气库H2井井深为3191 m,日产气量为5×104m3/d,油管内径为88.9 mm,气体成分为甲烷,阀门关阀方式为线性,井口温度为25 ℃,关阀时间为50 s,井口气相动力黏度为10.318 μPa·s,井口气相密度为1.138 kg/m3,采用的是储气井ϕ50.8 mm阀门,电压频率为50 Hz,现场调阀模式采用无模型触探式的关阀方式,即直接通过远程控制,一次性线性全关至0%。

    阀门阻力系数曲线由厂家交付阀门时提供,储气库采–关常用ϕ50.8 mm和ϕ 38.1 mm阀门,采用三次样条插值方法,将阀门阻力系数离散点平滑成曲线,计算波动压力模型时,首先计算相应开度下的阀门阻力系数,作为模型求解的基础条件;呼图壁储气库井区气井阀门的有效阻力系数区间约为20~80,当阀门开度在有效阻流区间内时,适当延长阀门调节时间,达到减小波动压力峰值的目的;当阀门开度在无效阻流区间内时,通过减少调阀时间,提高调阀效率,结果如图2所示。随着井深及日产气量的增大,油管内气相压力波速增大,井底压力响应时间缩短。

    图  2  H2 井 ϕ50.8 mm 阀门和 ϕ38.1 mm 阀门阻力系数关系曲线
    Figure  2.  Resistance coefficient curves of ϕ50.8 mm and ϕ38.1 mm valve in the Well H2

    不同日产气量条件下的井深与油管内气相压力波速关系如图3所示。井深不大于500 m时,随着井深增大,日产气量为5×104 m3/d的压力波速趋于恒定,而日产气量为20×104 m3/d的压力波速急剧增大。井深大于500 m时,随着井深增大,压力波速增大趋势显著。压力波速受气体可压缩性的影响较大,井深为1 000 m时,日产气量从5×104 m3/d变化至20×104 m3/d,压力波速从88.95 m/s增至729.81 m/s,增大720.5%。准确求解油管气相的压力波速,不仅能提高波动压力求解精度,还可以指导阀门调节时间间隔的制定。

    图  3  不同日产气量条件下井深与压力波速的关系
    Figure  3.  Relationship between well depth andpressure wave velocity under differentdaily gas productions

    日产气量对井口波动压力的影响如图4所示。随着日产气量增大,油管内气体动能增大,波动压力呈现增大趋势。由于摩擦阻力消耗了波动压力传播能量,井底波动压力同井口比较,井底波动压力同时受气体流速和油管沿程摩擦阻力的影响,因此井口比井底受波动压力的影响更大,更易发生油管疲劳损坏,过大的井口波动压力会缩短阀芯的使用寿命。

    图  4  储气井日产气量对油管波动压力的影响
    Figure  4.  Influence of daily production of gas storage well on the fluctuation pressure for the oil pipe

    油管深度对油管气相有效压力的影响如图5所示,油管内气相有效压力为波动压力与稳态压力代数和,井底的波动压力相比井底初始稳态压力较小,井底波动压力对有效压力影响不显著。当波动压力从井口向井底传播时,波动压力受油管内壁摩擦阻力影响,井口波源处的波动压力远大于井底波动压力,距井口距离越小,油管柱受波动压力冲击载荷影响越大。

    图  5  油管深度对油管波动压力的影响
    Figure  5.  Influence of depth on the fluctuation pressure for the oil pipe

    油管内径对井口油管气相波动压力的影响如图6所示。油管内径减小使相同流量的气体流动速度增大,单位时间内转换的气体动能增大,从而使油管气相波动压力呈现增大趋势。当油管内径从88.9 mm变化至50.8 mm,油管内气相波动压力峰值从1.22 MPa增大至1.43 MPa,油管内气相波动压力峰值增大了17.21%,油管内径的变化对波动压力影响不大,可以将其作为次要因素考虑。

    图  6  油管内径对井口波动压力的影响
    Figure  6.  Influence of inner diameter of oil pipe on the wellhead fluctuation pressure

    随着阀门关闭时间延长,波动压力峰值减小,且达到峰值的时间增长(见图7)。相比于液相,气相可压缩性大,更易使波动压力能量耗散,波动压力周期峰值衰减更快。关阀时间为4 s时,气相波动压力峰值在第3个周期(60 s)时衰减至0,波动压力对油管影响消失。当关阀时间从4 s变化至119 s时,油管气相波动压力峰值从1.58 MPa降至1.31 MPa,油管气相波动压力峰值降低了17.1%,关阀时间对波动压力影响显著,延长关阀时间可减小波动压力峰值。在日产气量较大的情况下,适当延长调阀时间,有利于减小阀芯的疲劳损坏程度。

    图  7  关阀时间对油管波动压力的影响
    Figure  7.  Influence of valve closing time on the fluctuation pressure for the oil pipe

    呼图壁储气库H2 水平井高产气量采–关转换过程中井口压力数据见表1。油管内气相波动压力受气体流量和关阀时间影响较大,采用15 s线性关阀方案,日产气量从50×104 m3/d增至600×104 m3/d,井口波动压力从9.21 MPa增大至19.26 MPa,井底最大波动压力峰值从0.61 MPa增大至4.62 MPa。结合阀门型号、油管承压能力、阀门阻流系数、井身结构等基础参数,利用油管气相波动压力数学模型,按照阀门开度区间绘制阀门开度曲线,可得到适用于现场的最优阀门控制方法。

    表  1  高产气量注–采转换井口压力分析
    Table  1.  Analysis of wellhead pressure for injection and production conversion for high gas production
    气体流量/
    (104m3·d−1
    井口波动
    压力/MPa
    井口初始
    稳压值/MPa
    有效井口
    压力/MPa
    井底波动
    压力/MPa
    井底初始
    稳压值/MPa
    有效井底
    压力/MPa
    501.218.009.210.6121.5522.16
    1501.668.009.660.8321.5522.38
    2504.938.0012.932.3721.5523.92
    3507.258.0015.253.2621.5524.81
    60011.268.0019.264.6221.5526.17
    下载: 导出CSV 
    | 显示表格

    1)针对储气库采–关作业工况,建立了油管气相波动压力数学模型,该模型精度高,可适用于关–采、关–注以及注–关等注采转换过程中的波动压力计算。

    2)相比于低产气井,波动压力冲击载荷对高产气井影响更为显著,需要适当延长调阀时间。随着油管深度增大,井底波动压力呈现减小趋势,浅井井底波动压力冲击载荷相比于深井更大。

    3)现场储气库采–关作业时,可根据阀门阻流区间确定阀门调节速度。储气库阀门开度在有效阻流区间内时,可适当延长阀门调节时间减小波动压力峰值;阀门开度在无效阻流区间内时,可通过减少调阀时间提高调阀效率。

  • [1]

    Li Yufei,Fu Yongqiang,Tang Geng,et al.Effect of weak bedding planes on wellbore stability for shale gas wells[R].SPE 155666,2012.

    [2] 刘忠飞,何世明,全家正,等.川西首口页岩气水平井钻井技术[J].石油钻采工艺,2014,36(1):18-21. Liu Zhongfei,He Shiming,Quan Jiazheng,et al.Drilling technology for the first shale gas horizontal well in West Sichuan[J].Oil Drilling Production Technology,2014,36(1):18-21.
    [3] 王红岩,刘玉章,董大忠,等.中国南方海相页岩气高效开发的科学问题[J].石油勘探与开发,2013,40(5):574-579. Wang Hongyan,Liu Yuzhang,Dong Dazhong,et al.Scientific issues on effective development of marine shale gas in Southern China[J].Petroleum Exploration and Development,2013,40(5):574-579.
    [4] 马天寿,陈平.基于CT扫描技术研究页岩水化细观损伤特性[J].石油勘探与开发,2014,41(2):227-233. Ma Tianshou,Chen Ping.Studies of meso-damage characteristics of shale hydration based on CT scanning technology[J].Petroleum Exploration and Development,2014,41(2):227-233.
    [5]

    Aadnøy Bernt Sigve.Stresses around boreholes drilled in sedimentary rocks[R].SPE 17119,1987.

    [6]

    Dusseault M B,Gray K E.Mechanisms of stress-induced wellbore damage[R].SPE 23825,1992.

    [7]

    Ong S H,Roegiers J C.Influence of anisotropies in borehole stability[J].International Journal of Rock Mechanics and Mining Sciences Geomechanics Abstracts,1993,30(7):1069-1075.

    [8]

    Okland D,Cook J M.Bedding-related borehole instability in high-angle wells[R].SPE 47285,1998.

    [9]

    Al-Bazali T M,Zhang J,Wolfe C,et al.Wellbore instability of directional wells in laminated and naturally fractured shales[J].Journal of Porous Media,2009,12(2):119-130.

    [10]

    Lee H,Ong S H,Azeemuddin M,et al.A wellbore stability model for formations with anisotropic rock strengths[J].Journal of Petroleum Science and Engineering,2012,96/97:109-119.

    [11] 金衍,陈勉.井壁稳定力学[M].北京:科学出版社,2012:13-66. Jin Yan,Chen Mian.Wellbore stability mechanics[M].Beijing:Science Press,2012:13-66.
    [12] 刘向君,叶仲斌,陈一健.岩石弱面结构对井壁稳定性的影响[J].天然气工业,2002,22(2):41-42. Liu Xiangjun,Ye Zhongbin,Chen Yijian.Influence of rock weak plane texture on sidewall stability[J].Natural Gas Industry,2002,22(2):41-42.
    [13] 袁俊亮,邓金根,蔚宝华,等.页岩气藏水平井井壁稳定性研究[J].天然气工业,2012,32(9):66-70. Yuan Junliang,Deng Jingen,Yu Baohua,et al.Wellbore stability of horizontal wells in shale gas reservoirs[J].Natural Gas Industry,2012,32(9):66-70.
    [14]

    Lu Yunhu,Chen Mian,Jin Yan,et al.Influence of porous flow on wellbore stability for an inclined well with weak plane formation[J].Petroleum Science and Technology,2013,31(6):616-624.

    [15] 刘志远,陈勉,金衍,等.裂缝性储层裸眼井壁失稳影响因素分析[J].石油钻采工艺,2013,35(2):39-43. Liu Zhiyuan,Chen Mian,Jin Yan,et al.Analysis on wellbore stability of open hole in fractrued formation[J].Oil Drilling Production Technology,2013,35(2):39-43.
    [16] 刘志远,陈勉,金衍,等.多弱面地层水平井裸眼井壁垮塌量分析[J].石油勘探与开发,2014,41(1):102-107. Liu Zhiyuan,Chen Mian,Jin Yan,et al.Calculation model for bore-hole collapse volume of a horizontal openhole in multiple-weak-plane formation[J].Petroleum Exploration and Development,2014,41(1):102-107.
    [17] 温航,陈勉,金衍,等.钻井液活度对硬脆性页岩破坏机理的实验研究[J].石油钻采工艺,2014,36(1):57-59. Wen Hang,Chen Mian,Jin Yan,et al.Experimental research on brittle shale failure caused by drilling fluid activity[J].Oil Drilling Production Technology,2014,36(1):57-59.
    [18]

    Chenevert M E,Gatlin C.Mechanical anisotropies of laminated sedimentary rocks[J].SPE Journal,1965,5(1):67-77.

    [19]

    Lee H,Chang C,Ong S H,et al.Effect of anisotropic borehole wall failures when estimating in situ stresses:a case study in the Nankai accretionary wedge[J].Marine and Petroleum Geology,2013,48:411-422.

    [20] 蔡美峰,何满潮,刘东燕.岩石力学与工程[M].北京:科学出版社,2009:104-109. Cai Meifeng,He Manchao,Liu Dongyan.Rock mechanics and engineering[M].Beijing:Science Press,2009:104-109.
  • 期刊类型引用(1)

    1. 郭海伟. 基于增量学习的储气库井底压力快速计算方法. 断块油气田. 2025(02): 292-299 . 百度学术

    其他类型引用(0)

计量
  • 文章访问数:  3520
  • HTML全文浏览量:  110
  • PDF下载量:  3701
  • 被引次数: 1
出版历程
  • 收稿日期:  2014-03-16
  • 修回日期:  2014-03-16
  • 刊出日期:  1899-12-31

目录

/

返回文章
返回