导流能力

导流能力（共7篇）

导流能力 篇1

摘要：裂缝导流能力测定实验是石油工程主干专业课采油工程的必修实验内容。原实验仪器过于庞大,占用空间大,且操作不便,安全性差。针对教学实验的特点,对仪器进行了全面改进,改进后的仪器结构紧凑,占用空间小,操作方便,教学效果良好。

关键词：裂缝导流能力,采油工程,实验仪器,研制

裂缝导流能力是油层压裂效果评价的重要技术指标,裂缝导流能力测定实验是石油工程专业本科生主干专业课采油工程的必做实验内容。目前石油高校使用的裂缝导流能力实验仪存在以下5个缺点。(1)仪器体积庞大,一台仪器占用面积超过2.0 m2。(2)笨重,仪器重量超过300 kg。(3)操作费力,导流室是主要活动部件,做实验时需要频繁拆装、填砂。原仪器的导流室重约15 kg,搬动、拆装非常困难。(4)闭合压力的加载操作复杂、不安全。原仪器闭合压力的加载方式是采用油压机,加压时需要进油阀和泄油阀同时开启,且进、出油量达到完全平衡方能使压力值控制在稳定数值,加压控制比较困难。一旦操作中出现失误,很容易出现压力快速上升、损坏压力表的现象,甚至危及人身安全。(5)仪器结构不紧凑。供气源采用氮气瓶供气,流程较长,占用空间大。

基于原裂缝导流能力测定实验仪存在的问题,笔者进行了全面改进,研制出一台结构紧凑、占用空间小、操作安全方便的新型裂缝导流能力测定仪。新研制的裂缝导流能力测定仪总重量不超过50 kg,占用面积不超过0.5 m2,导流室重量不超过4 kg。操作方便、安全,仪器功能完全能够满足实验教学需要。

1 裂缝导流能力测量原理及实验流程

所谓裂缝导流能力就是测定人工压开的裂缝充填支撑剂以后,允许流体通过的能力,实验流程如图1所示。其主要部件为导流室(如图2所示)。

导流室中充填支撑剂,支撑剂上方加压,模拟地层的闭合压力,气体通过导流室,测定入口和出口两端的压差和气体流量,即可计算不同闭合压力下裂缝的导流能力。计算公式为:

式中:k Wf–裂缝导流能力;μ–通过裂缝的气体黏度;Q–通过裂缝的气体流量;L–导流室两测压点间的距离;W–支撑剂厚度(缝宽);∆p–导流室进出口压差。

2 仪器设计及改进

2.1 导流室的改进

导流室是实验操作的重要部件,实验过程中需要从流程中拆下,起出活塞,在导流槽中填入支撑剂。原仪器的导流室尺寸为24 cm×10 cm×8 cm,重量约为15 kg。拆装、搬动非常困难。为便于操作,在保持结构不变的前提下,缩小尺寸,新设计的导流室尺寸为15 cm×5 cm×5 cm,改进后的导流室重量只有4 kg,学生可以很容易地进行拆装和填砂操作。尺寸缩小后,功能不变。

2.2 闭合压力加载方式的改进

原导流能力实验仪的闭合压力加载方式是由油压机完成的,油压机设备庞大,控制困难,操作不当就会出现压力快速升高的现象。新型导流仪将闭合压力加载方式改成活塞式水力控制。轻轻转动活塞增压手柄,就可以很轻松的加载,操作轻便、安全,而且占用空间大大减小。改进前后的加载设备如图3所示。

2.3 供气源的改进

气源可以用氮气也可以用空气。原仪器供气是采用氮气瓶,体积较大,不仅增加了设备占用面积,也增加了氮气购买成本。由于实验过程中所需的供气压力不高,将气源由氮气瓶改成压力容器供气,压力容器体积为氮气瓶的1/5,最大承压20 MPa,容气量完全可以满足实验的需要。由于容器体积小,可以直接固定在仪器台上,不再额外占用空间,设备非常紧凑。利用一台增压泵可为多台仪器充气,节省了成本。

2.4 气体流量计的改进

原仪器的气体流量计为浮子流量计。浮子流量计的缺点是,流量一旦超过量程,浮子很容易卡在流量管的顶端,影响流量计正常使用。新研制的仪器改成了皂膜流量计,一方面皂膜流量计精度比浮子流量计高;另一方面,皂膜流量计耐用,不会出现因超量程引起损坏的问题,特别适合教学使用。改进前后的仪器对比图如图4所示。

3 仪器改进后教学效果

改进后的仪器已用于实验教学,达到了预想的教学效果。(1)仪器的操控性增强。以前只有男生才能够操作的导流室的拆装,现在每个学生都可以轻松自如地进行操作。(2)仪器的安全性增强。加载闭合压力时,完全避免了升压过快、过高的问题,消除了安全隐患。(3)每组人数减少。由于原设备过于庞大,笨重,根据空间大小,实验室只能购置两台,且只能放置在一层。改进后的仪器结构紧凑,且移动方便,目前设计加工了5台。每组人数由原来的6人/台降至2人/台,大大改善了教学效果。

4 结束语

新型裂缝导流能力测定仪占用空间大幅度减小,占用面积由2.0 m2降为0.5 m2,重量由300 kg降为50 kg,仪器操作更加轻便、省力,实验更加安全、可靠。新型裂缝导流能力测定仪已使用400余人次。教学效果良好。

导流能力 篇2

1 变导流能力研究

裂缝的导流能力是指油层条件下填砂裂缝渗透率与裂缝宽度的乘积[11],即

式(1)中:FCD为裂缝导流能力,10-3μm2·m;kf为填砂裂缝渗透率,10-3μm2;wf为填砂裂缝宽度,m。

1986年,Poulsen指出填砂裂缝的导流能力沿裂缝方向不是恒定的,而是不断变化的,其变化形式有直线变化、指数变化和对数变化三种趋势[7]。依据实验结果,通过曲线回归,可以得到裂缝导流能力随距离的指数变化形式[8,12]。图1为裂缝在井筒处的导流能力为200×10-3μm2·m时,裂缝导流能力沿缝长的变化情况示意图。

式(2)中:FCD(xD)为裂缝在无因次位置xD处的导流能力,10-3μm2·m;FCD0为裂缝在井筒处的导流能力,10-3μm2·m;xD为沿裂缝方向的无因次距离(xD=x/xf),m。

2 模型建立及求解

气藏中直井在压裂后,渗流模式将发生变化,可分为地层内的流动和裂缝内的流动两部分,气体首先由地层流向裂缝,再由裂缝流向井底,压裂裂缝(图2)的导流能力可分为无限导流能力与有限导流能力两种[13]。

2.1 地层内的流动

气井经过压裂后,气藏储层内部的气体渗流进入裂缝,可形成等压椭圆柱面,采用拟压力表示则有下式[14]。

式中:pe为地层压力,MPa;pw为油藏和裂缝区域衔接处压力,MPa;Qsc为标准状态下的气井压裂后产量,m3/d;T为气藏温度,K;k为地层渗透率,10-3μm2;h为地层厚度,m;re为油藏泄流半径,m;xf为压裂裂缝半长,m;μg为气体黏度,m Pa·s;Zg为气体偏差系数;m(p)为拟压力函数。

2.2 支撑剂有效渗透率

气井压裂后,气体由储层首先渗流进入裂缝,再由裂缝流到井筒。由于裂缝内气体高速流动,将产生高速非达西现象,使得裂缝内支撑剂的渗透率降低,可以通过雷诺数对支撑剂渗透率进行修正,计算得到其有效渗透率[15]。

式中,R为摩尔气体常数,;Bg为气体体积系数,无因次;hf为裂缝高度,m;wfp为裂缝宽度,m;β=a/kfb为非达西流动系数,m-1;ρg为气体密度,kg/m3;v为气体流速,m/s;μg为气体黏度,m Pa·s;kf为支撑剂渗透率,10-3μm2;kf-eff为支撑剂有效渗透率,10-3μm2;Re为雷诺数;a,b为常数,对于不同类型的支撑剂其数值不同,Cooke[16]通过实验确定常数a、b的值(表1)。

2.3 裂缝内的流动

气体进入裂缝后,由裂缝流向井底,气体在裂缝内为高速非达西流动,应该采用支撑剂的有效渗透率进行产能计算,假设裂缝为有限导流能力,则有[14]:

计算中采用裂缝的有效导流能力,即采用利用雷诺数修正后的有效渗透率进行计算。

式(13)中:pwf为气井井底压力,MPa。

由上述公式即可推得变导流能力情形下的垂直裂缝气井的非达西产能公式。

在气井井底流压和气藏压力范围内,可认为气体的黏度和气体偏差因子为常数,因此可移出积分号[17],故气井产量为

2.4 无因次采气指数

采气指数(productivity index)是单位生产压差下的气井产气量。对于拟稳态径向流情形下直井的采油指数可写为无因次采油指数JD的函数[18]。

式(16)中:PI为采油指数,m3/(d·MPa);α为单位换算系数;JD为无因次采油指数。

由式(15)按照式(16)对无因次采油指数的定义,可以得到变导流能力直井非达西无因次采气指数JDV。

式(17)中,JDV为压裂直井非达西无因次采气指数。

2.5 裂缝壁面伤害

在压裂过程中,由于压裂液的滤失等原因,地层不可避免地将受到伤害,导致裂缝壁面处储层渗透率下降,进而降低压裂井产能,这称之为裂缝的表皮效应,如图3所示[19]。1978年,Cinco-Ley等[20]提出利用裂缝壁面表皮系数Sf(fracture face skin)来定量表征裂缝的壁面污染效应,式(18),公式中假定沿裂缝方向储层污染深度ws(fracture face skin zone extent)恒定。

式(18)中:Sf为裂缝壁面表皮系数;ws为裂缝处储层污染深度,m;ks为受污染储层渗透率,10-3μm2。

实际计算中,参考相关实验结果,对参数进行了合理赋值。依据压裂液滤失实验取ws=0.1 m[21];依据压裂液岩心伤害率实验取k/ks=3[22],该值表示裂缝壁面储层污染程度,该值越大,储层污染越严重。

图4为储层在不同污染深度与污染程度情形下裂缝壁面表皮系数与裂缝半长的关系曲线。图中可以看出,裂缝壁面表皮系数随裂缝半长的增加而减小,储层的污染深度越深、污染程度越重,则裂缝壁面表皮系数越大。

当考虑裂缝壁面储层污染效应时,此时的无因次采气指数可由式(19)近似求得[19]。

式(19)中,JDV,D为压裂直井考虑裂缝壁面污染时的无因次采气指数。JDV|S=0为压裂直井裂缝壁面Sff=0时的无因次采气指数。

此时,气井产量为

2.6 模型求解

由于模型中采用雷诺数修正支撑剂有效渗透率的方法来表征裂缝内气体的非达西效应,因此气井产量无法直接求解,需要循环迭代进行计算,基本步骤如图5所示[10,15]。

(1)给定雷诺数的初值,计算裂缝内支撑剂的有效渗透率,式(6)。

(2)计算气井产量,式(20)。

(3)根据气井产量计算气体在裂缝内的流速,式(9)。

(4)用计算得到的气体流速计算新的雷诺数,式(7);如果初始雷诺数与计算的雷诺数在误差允许范围内,则计算终止;否则,继续循环求解。

3 算例

如图6,气井产量随着支撑剂渗透率的增加而增加,尤其当支撑剂渗透率在10~100 D范围内时增加较快;在支撑剂渗透率大于1 000 D后,产量随支撑剂渗透率的增加的增幅较小,且不同情形下的气井产量趋于接近。计算表明裂缝支撑剂渗透率较高时,裂缝导流能力的变化对产量影响较小;水力裂缝不仅提供了渗流通道,更减缓了裂缝内气体的非达西效应,这是压裂增产的原因之一。实际压裂施工中,支撑剂渗透率多集中在10~500 D之间,在此区间内支撑剂渗透率对产量的影响较大,不同情形下的产量差别也较大。考虑变导流能力时,产量会下降;同时裂缝内气体的非达西效应也使得产量急剧下降。

由图7可知在本算例中雷诺数随着支撑剂渗透率先增加后降低,存在一个最大值;且变导流能力时的雷诺数总小于定导流能力时的雷诺数。变导流能力时,裂缝的导流能力较定导流能力时低,气体在裂缝内的渗流阻力大,使得气井产量降低,裂缝内气体流速降低,因此雷诺数较小。支撑剂渗透率较低时(小于100 D),与气藏的供气能力相比,裂缝的导流能力不足;裂缝导流能力的增加,将减小裂缝内气体流动阻力,使得产气量大幅增加,裂缝内气体流动速度也大幅增加,此时气体流速是控制雷诺数的主导因素;当支撑剂渗透率较高时(大于100 D),随着裂缝内渗透率的增加,雷诺数呈现递减趋势,且两种情形下的雷诺数趋于重合;由产量曲线可知,此时气井产量增加缓慢且趋于稳定,气藏的供气能力是限制产量增加的主因,即裂缝内气体流速趋于定值,支撑剂渗透率是雷诺数的主导因素,渗透率增加,雷诺数降低。

4 敏感性分析

设定裂缝渗透率Kf=200 000 m D,由表2中基本参数计算并绘制不同裂缝长度及裂缝宽度情形下的产量曲线。

4.1 裂缝长度的影响

在该算例中,四种情形下气井产量都随裂缝半长增加而增加(图8),但曲线形态却存在一定差别。定导流能力、达西情形下气井产量最高,变导流能力、非达西情形下气井产量最低,且二者间的产量差值随着裂缝半长的增加而增加。算例中,达西情形下的产量均高于非达西情形下的产量,表明裂缝内气体的非达西效应对产量影响严重。

4.2 裂缝宽度的影响

裂缝宽度较小(Wf=0.002 m、0.004 m)时,气井产量随着裂缝长度先增加后降低(图9),此时的雷诺数也具有相似的规律(图10)。裂缝宽度较小时,气体在裂缝内高速流动、具有较大阻力,随着裂缝半长增加,初期产量上升;但裂缝过长时,地层内的气体大量流入裂缝,裂缝内气体的高速非达西流动消耗大量能量,反而使得产能降低,也表现为此时雷诺数随裂缝半长增加而减小。裂缝宽度较宽时,产量则是随着裂缝长度增加而增加,但前期产量增加较快,后期产量增加较缓慢;雷诺数也随裂缝半长增加而增加,但后期基本趋于稳定;此时裂缝的导流能力充足,而气藏的供气能力是限制产能的主导因素。

压裂气井的产能由气藏的供气能力和裂缝的导流能力共同决定,只有二者相协调时产能才能最大。对于给定的气藏,其供气能力一定,因此需要优化裂缝的导流能力;在裂缝导流能力较低时,会存在最优缝长,使得产量最大,缝长过长或过短都会使产量降低。

5 结论

裂缝的导流能力沿缝长方向不是固定不变的,而是不断递减的,可采用指数、对数、线性递减形式来表示导流能力的递减,在产能公式中也应考虑裂缝导流能力的变化。

裂缝变导流能力时的产能较定导流能力时的产能低,在支撑剂渗透率较低或较高时,两种情形下的产能差值较小,但在实际的支撑剂渗透率范围内,两种情形的产能有明显差异。

由于气井裂缝内气体的高速非达西流动,会降低支撑剂的渗透率,可采用雷诺数对支撑剂渗透率进行修正,得到有效渗透率,从而进行产能计算。考虑气体非达西效应时,气藏产量也将降低。

裂缝内支撑剂渗透率对产量有明显影响,增加支撑剂渗透率,可提高裂缝导流能力,降低气体在裂缝内的渗流阻力,减弱非达西效应,进而提高产量。裂缝宽度对产量影响也有很大影响,气井产量随着裂缝宽度增加而增加。因此实际压裂施工中,应尽量选用高渗透率支撑剂,且增加缝宽。

导流能力 篇3

1 剪切裂缝形成条件

闭合的天然裂缝(或弱面)受远场水平最大主应力σH和水平最小主应力σh作用,与水平最大主应力夹角为θ(0<θ≤π/2),如图1所示。根据二维线弹性理论,作用在天然裂缝面上的法向正应力σn和剪切应力τ分别为

当天然缝内压力较小时,作用于裂缝面上的法向正应力为压应力,裂缝处于闭合状态;同时较大的剪切应力促使缝面发生剪切滑移[5]。

式(3)中,τ0为天然裂缝面的黏聚力;Kf为天然裂缝面的摩擦系数。

缝面有效剪切应力为

压剪裂缝的破裂准则为:

式(5)中,KΠC为岩石Ⅱ型断裂韧性;KΠ为应力强度因子,。

将式(1)、式(2)和式(4)代入式(5),可得天然裂缝发生剪切滑移,并形成分支裂缝的临界压力

考虑水力裂缝仅诱导天然裂缝剪切滑移但不扩展,即,所需的裂缝净压力为

当θ=arccot Kf/2时净压力为最小值,天然裂缝最易发生剪切滑动,最小净压力pmin为

式中,Δσ=σH-σh。而当θ=π/2时净压力达最大值,最大净压力pmax为

如果天然裂缝面非黏结,即τ0=0,则剪切滑移所需净压力为水平主应力差Δσ。由式(7)可知,当水力裂缝与天然裂缝相交后,决定天然裂缝是否发生剪切滑移的影响因素包括水平主应力差、天然裂缝角度、天然裂缝面摩擦系数和黏聚力。在低应力差、低角度或低摩擦系数的条件下,天然裂缝发生剪切滑移所需的净压力较小,从而更易形成错位裂缝。

2 错位裂缝表征

2.1 缝面错位程度

裂缝面在剪切应力作用下发生滑移,根据断裂力学中Westergaard函数,无限大介质中Ⅱ型裂缝面(单面)剪切位移表达式为[11]

式中:l为裂缝半长;x为沿缝长方向任意点坐标;κ为Kolosov常数,κ=3—4ν;G为剪切模量,;E,v分别为岩石杨氏模量和泊松比。

在裂缝中心处(x=0)有最大错位程度

将式(4)带入式(8)得压剪情况裂缝错位程度

当p=pnet+σh,τ0=0时,上式简化为

压剪裂缝与张开裂缝的临界条件是p=σn,此时缝内净压力达到,缝面之间没有摩擦阻力,可有最大错位程度

由式(9)可知,剪切裂缝的错位程度主要受岩石的杨氏模量、泊松比、地应力差、夹角、天然裂缝长度及缝面摩擦系数等因素影响。当地层应力差较大,杨氏模量、泊松比和缝面摩擦系数较低时,天然裂缝易形成较大的错位,同时天然裂缝(或弱面)较长也有利于形成较大错位,天然裂缝在压剪状态下形成错位的有利角度是arccot Kf/2。此外,提升天然裂缝内净压力,可降低天然裂缝面之间的摩擦阻力,增加了剪切滑移发生机会。因此,页岩储层压裂施工过程中,高排量或高净压力有利于天然裂缝剪切滑移,形成较大位移的错动。

2.2 缝面粗糙度

在非支撑剂充填裂缝中,裂缝面接触后的残存缝隙是疏导流体的主要通道,裂缝面的粗糙度是决定残余孔隙空间大小的一个重要因素,但如果没有对缝面粗糙程度的准确描述,很难去确定何种剪切裂缝易提供高导流能力。实际地层裂缝的凸起大小和错位程度难以预测,缺少直接的现场数据。室内采用机械劈裂加工而成的岩石裂缝形态及其表面轮廓与实际水力裂缝有很大差别[12]。在页岩气藏中,清水压裂主要是开启并沟通天然裂缝,而通过描述页岩体内部天然裂缝表面特征,沿天然裂缝或层理加工形成的导流岩心如图2所示,岩心长17.7 cm(长)×3.8 cm(宽)×2—3 cm(厚)。测试其在整合或不同错位程度时的导流能力,有助于准确评估地层中缺少支撑剂的裂缝的实际导流能力。

裂缝表面轮廓可以用分形维数来描述[13]。两点间凸起增量平方均值函数

式(11)中,r为测点间距;Z(xi)为xi点凸起高度;N(r)为测点间距为r时的计算次数。以lg[V(r)]为纵坐标,lg(r)为横坐标绘图,在r值较小时,存在近似线性段,线性段的斜率为β,则分形维数D与β的关系为:D=2-β/2。粗糙度JRC与缝面分形维数之间的经验关系式为

通过上述方法,首先对页岩中天然裂缝面或弱面上的凸起进行统计,可利用激光扫描仪测量(如图3所示的裂缝表面轮廓),然后计算裂缝面的分形维数D,再利用式(12)就可以得到表征页岩剪切裂缝表面凹凸不平程度的粗糙度系数JRC。

3 剪切裂缝导流能力

3.1 页岩性质和试验方法

根据上述理论对四川盆地某区块页岩气层进行了测试分析。地层水平最大主应力43 MPa、水平最小主应力38 MPa,天然裂缝角度7°—15°、长度5—10 m。页岩静态杨氏模量在11.4—18.4 GPa之间,平均14.2 GPa,泊松比0.3。矿物组成(X光衍射法)包括黏土39.9%(高岭石5.4%、绿泥石7.2%、伊利石21.7%和伊蒙混层5.6%),石英41.3%,长石5.7%和碳酸盐13.1%。

将地层参数带入式(10)得裂缝可形成的最大剪切错位程度在0.39—1.60 mm。并对13块岩心裂缝表面的凸起高度进行测试统计,如表1所示,凸起最高在1.94—6.12 mm,平均高度1.88 mm,偏差17.9%—54.7%。再利用式(12)计算得到粗糙度在6.2—54.9,呈现出裂缝表面轮廓较强的随机性。测试的页岩缝面凸起高度与其它岩性存在差异,Fredd[12]室内砂岩劈裂形成的裂缝面凸起为3.05—8.38 mm,平均2.29 mm,而Warpinski报道的矿场挖掘结果为均质砂岩1.02 mm,而含煤泥夹层非均质砂岩达4.57 mm。

剪切错位裂缝的导流能力测试利用FCES-100导流仪(API标准)完成,裂缝承受的有效闭合应力在10—60 MPa,每个压力点24 h。由于页岩中含有较多的伊利石和高岭石,容易脱落形成可移动的粘土颗粒,堵塞剪切裂缝的残余缝隙,降低导流能力,需要进行防膨胀和长效黏土稳定处理,因此应用了含有黏土稳定剂(1%KCl+1%CS—2)的清水作为测试流体,流量2—5 m L/min。

3.2 实验结果与分析

图4为2条整合裂缝和5条错位裂缝的导流能力综合图。整合裂缝在10 MPa闭合应力作用下,导流能力最高为0.26μm2·cm,而当应力超过20MPa后,导流能力微乎其微,整合裂缝表现出较强的应力敏感性。而一般地层应力均超过20 MPa,因此整合的非支撑裂缝并不能提供稳定、足够的导流能力。当裂缝发生剪切错位后,导流能力有了显著提高,但在每个闭合应力下,导流能力均表现出较大的变化范围,可达1—2个数量级,在10 MPa闭合应力下最高可达56.70μm2·cm,最低为1.57μm2·cm,部分裂缝在60 MPa闭合应力下仍有0.018μm2·cm。从整合裂缝到剪切错位裂缝,导流能力提高2个数量级以上,但导流能力并不是随着错位程度增大而持续增加,如图5所示,从低错位程度到高错位程度,导流能力整体显著提高,而局部存在波动性变化,主要与缝面凸起的随机分布相关。综上所述,缝中缺少支撑剂时,裂缝需要通过剪切错位来形成残余缝宽,从而提供导流能力,但错位程度的影响存在随机性。

剪切裂缝导流能力值表现出较大的变化范围,主要受裂缝表面粗糙程度(凸起大小和分布)的影响,实验结果表明粗糙度较大的裂缝更易形成高导流能力,但导流能力与粗糙度并没有绝对的相关性,仅闭合应力低于30 MPa时,导流能力与粗糙度表现出一定的相关性。低闭合应力下,缝面如有较高的凸起,相应的粗糙度也较大,容易形成宽度较大的连通通道,但随着闭合应力增大,凸起逐渐发生破碎,裂缝闭合程度不断加剧,不同粗糙度裂缝的宽度越接近。因此,粗糙度对剪切裂缝导流能力仅在低闭合应力条件下影响显著。

随着闭合应力增大,导流能力降低幅度与缝面凸起的破碎程度是一致的,实验后裂缝表面凸起破碎平均达0.5 mm,约降低25.2%(如图3所示)。凸起破碎和表面黏土矿物脱落易形成软化的泥层,如图6所示,实际可支撑裂缝的有效凸起高度要远远低于统计的剩余高度。而硬度高的砂岩(杨氏模量24.8—48.2 GPa)表面凸起破碎平均仅为0.254mm,约降低9%[12]。因此,在评价剪切裂缝导流能力时,除了缝面的粗糙度外,还要充分考虑岩性、矿物成分和岩石力学性质的影响,如图7所示,软硬页岩(杨氏模量分别为14.2 GPa、30.0 GPa)剪切错位裂缝的导流能力可相差一个数量级以上。

4 结论

(1)水力裂缝诱导天然裂缝发生剪切滑移的条件是裂缝内净压力超过地层水平主应力差。地层应力差、天然裂缝角度或摩擦系数越小,发生剪切滑移所需的净压力越小。高应力差,低杨氏模量和低泊松比地层中的天然裂缝易形成较大错位,压剪闭合状态下最有利的天然裂缝角度是arccot Kf/2。

(2)裂缝中缺少支撑剂时,需要通过剪切错位来形成残余缝宽,从而提供导流能力,但错位程度的影响存在随机性,且错位裂缝导流能力变化范围较大,可达1—2个数量级。粗糙度较大的裂缝更易形成高导流能力,但导流能力与粗糙度并没有绝对的相关性,仅闭合应力低于30 MPa时,导流能力与粗糙度表现出了一定的相关性。

导流能力 篇4

本文利用实际碳酸岩岩心分别进行硫沉积后和酸压后裂缝导流能力的测试实验, 通过对比硫沉积前后和酸化前后裂缝导流能力的变化, 主要探讨酸压过程对硫沉积堵塞后的裂缝导流能力改善机理。

1 实验装置与测试方法

1. 1 实验装置

酸蚀裂缝导流能力实验装置是在原有的裂缝导流能力仪的基础上发展而来的, 主要用来测定酸压或复合酸压施工中的几个关键参数, 图l为导流能力测量系统流程图。

①为ISCO泵;②为高压中间容器;③为高压岩心加持器;④为压差传感器;⑤为流量计量器

1. 2 测试方法

室内实验应用岩心流动的方法, 模拟硫沉积以及酸蚀硫沉积对裂缝导流能力的影响。

取两块造缝处理后的岩心 ( 编号分别为1#和2#) , 在进行实验时, 在酸液中加入一定浓度 ( 浓度大小不同的) 的固相微颗粒 ( 颗粒大小与硫颗粒相当, 可直接选择硫磺颗粒) , 经过酸蚀后, 测试固相颗粒存在对裂缝导流能力的影响。

( 1) 为了模拟碳酸岩类气藏气体的流动, 采用注氮气的方式, 并使用碳酸岩岩心。

( 2) 为了模拟硫颗粒的沉积, 在实验流程中加一段管线, 放置一定量的、研磨成20 ~40 目的硫磺粉。

( 3) 为了对比酸蚀裂缝与未酸蚀裂缝中硫颗粒沉积对导流能力的影响, 在实验过程中, 围压上升到2 MPa之后进行酸洗, 并静置一段时间。

1. 3 实验步骤

( 1) 设定温度28 ℃, 测定造缝岩心的直径、长度等基本静态参数。

( 2) 测试初始状态裂缝导流能力: 将待测岩样装入岩心夹持器中, 加2 MPa环压, 让干燥的氮气在一定的压差下通过样品, 待流量稳定后测得其进出口压力和氮气的流速。

( 3) 测试加入硫粉后裂缝导流能力变化: 将岩样依次加2 MPa环压, 用管线注入一定量的硫粉, 后依次升高围压至5 MPa、10 MPa、15 MPa测定其进口压力和出口处氮气的流速, 后将压力按原来压力值降低, 并测定进口压力和出口处氮气的流速。

( 4) 测试酸蚀硫沉积裂缝导流能力变化: 将岩样依次加2 MPa环压, 注入一定量15% 的盐酸, 其他步骤同上。

( 5) 实验数据整理, 设备关闭。

2 实验结果及处理

利用以上实验得到1#和2#岩心测得不同围压下、加入硫粉以及酸蚀前后的进出口压力和氮气的流速, 结果如表1 所示。

气体的体积流量随压力和温度的变化而不同。由于流体在岩石中渗流时, 在岩石长度L的每一断面的压力不同, 因而, 流过岩石的气体体积流量在岩石内各点上是变化的, 是沿着压力下降的方向不断膨胀, 增大。

考虑这一特征, 气体在岩石中任一点的流动状态必须用达西定律的微分形式表达:

式 ( 1) 中, Q为P压力下气体的体积流量 ( cm3/ s) ;Kg为气测渗透率 ( μm2) ; μ 为黏度 ( MPa·s) ; A为岩石样品的截面积A = W × b ( cm2) ; L为岩石样品的长度 ( cm) 。

即认为在微元dL上, 气体体积流量可变。由于dL增量为正时, dP增量为负, 为保证渗透率为正值, 在公式右端取负号。

如果气体在岩石孔隙中的渗流为稳定流, 气体流过的各断面上的质量流量不变, 那么根据波义尔-马略特定律, 等温条件下气体体积流量随压力的变化关系可表为

式 ( 2) 中: P0为大气压力 ( MPa) ; Q0为P0压力下的气体体积流量。

将 ( 2) 式代入 ( 1) 式后得到

对式 ( 3) 分离变量并积分

得到气测渗透率计算公式

式 ( 4) 中, P1为进口压力 ( MPa) ; P2为出口压力 ( MPa) ; μg为气体黏度 ( MPa·s) ; l为缝长 ( cm) ; W为缝宽 ( cm) 。

由此导出气测岩心裂缝导流能力计算公式

式 ( 5) 中: l为缝长 ( cm) ; W为缝宽 ( cm) 。

由于气体间黏度差较小, 定义气测裂缝导流能力DFG:

根据式 ( 6) , 将表1 实验结果及记录计算出不同实验条件下的裂缝导流能力 ( 表2) 。

注: P进为注气端注入压力, MPa; Q为P压力下的注入速度, mL/min。

3 实验数据分析

3. 1 硫沉积对裂缝导流能力影响

对1#和2#岩心加入硫粉后裂缝导流有能力测试结果均表明:

( 1) 当存在硫沉积时, 裂缝导流能力严重降低。例如, 1#岩心在加压前围压2 MPa时, 加入硫粉前裂缝导流能力为6 765. 97μm2·cm/ ( mPa·s) , 加入硫粉之后, 裂缝导流能力降低到1 312. 67 μm2·cm / ( mPa· s) , 接近原来的1 /5, 平均降低了3. 54倍。在15 MPa下, 1#岩心的裂缝导流能力降低了47%。

( 2) 裂缝导流能力随围压变化弹性增强, 塑性降低, 并且在围压上升过程中围压较低时裂缝导流能力对围压更为敏感, 在围压上升阶段裂缝导流能力明显高于围压降低阶段 ( 如图1 和图2) 。

2. 2 酸蚀对硫沉积后裂缝导流能力影响

由于实验过程中存在酸液的干扰, 液体的流动性较气体的差, 酸液在一定程度上阻碍了注入氮气的流动, 导致气测渗透率的大幅度降低, 再利用气测裂缝导流能力进行比较已失去意义, 酸蚀后裂缝导流能力的变化一般用裂缝能力降低系数来评价。

当注入酸液时, 气测岩心裂缝导流能力已不再准确, 故引入堵塞后裂缝能力降低系数概念, 即初始设定围压下的裂缝导流能力与最终设定围压下的裂缝导流能力的比值

式 ( 7) 中, ΔDFGaf为裂缝能力降低系数, 无量纲量;DFGiaf、DFGfaf分别为初始围压下和最终围压下的裂缝导流能力。

根据式 ( 7) 中的裂缝能力降低系数来评价酸蚀后裂缝导流能力的变化如表3 所示。

对存在硫沉积的1#和2#岩心进行酸蚀后裂缝导流有能力测试结果均表明:

( 1) 当进行酸蚀时, 对于存在硫沉积的裂缝, 其导流能力在围压的变化中表现出基本保持稳定的特征, 说明了酸蚀裂缝能够在一定程度上降低硫沉积对裂缝的堵塞, 进一步改善裂缝的导流能力; 例如硫沉积产生堵塞后的1#岩心在进行酸蚀后, 裂缝导流能力降低系数由3. 23 降低到了1. 0, 酸蚀后围压变化没有导致裂缝导流能力的进一步降低, 反映出裂缝的堵塞得以改善;

( 2) 对硫沉积造成堵塞后的裂缝进行酸蚀的结果也表现出裂缝导流能力随围压变化弹性增强, 塑性降低, 但裂缝导流能力比围压上升过程中裂缝导流能力受压力影响近于直线关系, 围压上升酸蚀裂缝导流能力降低 ( 如图3 和图4) 。

3 结论

( 1) 当存在硫沉积时, 裂缝导流能力严重降低。

( 2) 裂缝导流能力随围压变化弹性增强, 塑性降低, 并且在围压上升过程中围压较低时裂缝导流能力对围压更为敏感。

( 3) 酸蚀裂缝能够在一定程度上减少硫沉积对裂缝的堵塞, 进一步改善裂缝的导流能力, 所以在实际生产中遇到硫沉积堵塞裂缝时, 可进行酸化的方法改善裂缝的导流能力。

导流能力 篇5

在油气田开发过程中, 准确预测水平井产量对于水平井规划是至关重要的。用于水平井产能预测的数学模型可分为三类:①20世纪80年代后期至90年代初期, 基于油井无限导流能力的假设而推导出的简单解析模型;②20世纪90年代后推导出的关于有限导流能力油井的复杂解析模型;③考虑井筒水力因素的数学模型。

第一类模型由于使用方便而在工业中被广泛使用, 但由于没有考虑井内的摩擦压降, 所以对油井的生产能力预测值偏高。根据Hill和Zhu (2006) 的研究, 在一个厚度为100 ft (1 ft=30.38 cm) 、渗透率为100~200 mD (1 mD=1.02×10-2 μm2) 、长度为5 000 ft的油藏内, 一口井筒直径为4 in (1 in=25.4 mm) 的油井的摩擦压降可达到整个压降的2%~8%, 该值会随井筒长度的增加而增加, 并且在高渗层小井筒中会变得更加突出。

对第二类模型进行推导可解释井筒水力因素对油井产量的影响, 尽管该类模型将油层流入与油井流出紧密结合起来并可获得精确的预测结果, 但由于需要复杂的计算方法而阻碍了其在该领域的应用。

第三类模型可以给出比较精确的产量预测结果, 但由于该类模型中隐式计算程序和算法的迭代特性需要较长的CPU时间, 因此对于大多数石油工程师来说, 很难实现对该类模型较好的应用。

尽管有许多数学模型可进行水平井产能预测, 但由于它们共有的数学复杂性和数值处理特性而使得应用起来很不方便。因此, 人们非常渴望获得一个简单且能够精确预测水平井产能的数学模型。

本文将油藏流入和油井流出的实际情况紧密结合起来, 推导出一个简单的数学模型。该模型与Dikken (1990) 的模型类似, 不同点在于:①使用更接近实际的边界条件;②考虑从层流到紊流的过渡;③对于有限导流能力的井筒推导出一个解析解法。该模型较其他模型更加精确且使用方便。

2 数学模型

图1为水平井示意图。推导这个新模型的方法与推导Dikken模型所用的方法类似。Dikken模型的边界条件是在水平段尖端压降导数都为零, 而本文采用更加实际的边界条件, 即在任何压降为零的地方压降导数都为零。新模型的推导在附录A和B中 (本文略) 。

考虑稳流状态。当流体从水平段趾部流向跟部时, 即流体从油层流向井筒, 采出液的体积流量增加。当在临界距离 (从趾部到跟部) 时, 流率也达到临界值, 经过一个过渡后流态从层流变为紊流。假设临界雷诺数是2 000, 可以得到临界流率的表达式, 如下:

undefined

临界距离的表达式如下:

undefined

其中undefined (3)

当临界距离大于水平井段的长度时 (xc>L) , 流体处于层流状态, 井内总流率表达式为:

undefined

当xc

undefined

其中undefined (6)

undefined

其中范氏摩擦系数ff通过Chen (1979) 提出的相关公式计算得出。由于雷诺数是与流率有关的量, 因此, 在公式 (`5) 中采用误差方法来解QoH。图2给出了该方法的流程图, 该方法在电子表格中进行编码并运行迭代算法。

3 模型比较和讨论

将新模型与Joshi (1988) 、Butler (1994) 和Furui等人 (2003) 提出的常用模型进行比较, 所用数据来源于Chauvel等人 (1994) 提出的以下油田实例。

在澳大利亚的一个油田中钻1口水平井。这口井是8.5 in的裸眼井, 采用5.5 in预填砂筛管完井。井筒轨迹几乎是水平的, 总垂直深度将近6 ft, 水平段长2 438 ft, 位于131.2 ft厚的油层的中部。油层特性和井的参数在表1中作了简述。生产测井数据表明, 该井产量为5 677 bbl/d (1 bbl/d=0.159 m3/d) 。这个值与在地面上测得的产量 (5 660 bbl/d) 相差不大, 而在井底条件下的产量为5 990 bbl/d。当油井在略低于泡点压力条件下生产时, 几乎不存在自由气体。在井底和地上所测得的流率的一致性表明, 在衬管和裸眼井段之间的环形空间内实际上是没有流动的。压降从水平段趾部的3.5 psi (1 psi=6.895 kPa) 均匀变化到跟部附近的7.5 psi。

由于Chauvel等人 (1994) 的文章中没有给出有关油层渗透率和液流边界信息, 因此无法根据公式 (5) 直接计算油井产量。通过拟合井筒压力剖面, 根据附录A中公式A.33计算, 可估算出油层渗透率和液流边界条件。图1给出了拟合井筒压力剖面, 并可看出井筒压力从水平段趾部的930 psi降到跟部的925 psi。拟合使用的边界距离yb=500 ft, 水平渗透率kH=770 mD, 渗透率比kV/kH=1/3, 表皮系数s=0。

根据表1中的参数, 公式 (5) 给出1口油井的产量为6 823 bbl/d。井筒的液体流量可以通过附录B中的公式B.12计算出来。图3显示出由生产测井得到的流量剖面与计算得到的流量剖面相互比较的结果。由图可知, 计算得到的井筒流率在水平段趾部为0, 跟部为6 823 bbl/d, 计算结果除井底附近数值较高外, 与测井得到的剖面基本一致。事实上, 测井得到的流率在井底的影响值 (14%) 比在整个油井的平均影响值 (14.29%) 低, 导致这一误差的原因可能是由于原油开采之前井筒已经受到污染。

根据Joshi (1988) 、Butler (1994) 和Furui等人 (2003) 提出的模型及公式 (5) 所计算出的油

井生产率的比较在表2中作了简述, 除了公式 (5) , 其他模型得出的产量预测值均显著偏高, 公式 (5) 计算得出的结果误差为20.54%。

4 结论

本文提出了一个简单且能够精确预测有限导流能力水平井产能的解析模型。该模型与现有解析模型相比具有以下特点:①使用更加接近实际情况的边界条件;②考虑了井筒中从层流到紊流的过渡期。利用实际水平井数据对该模型与现有的3个数学模型进行了比较, 得出以下结论:

◇ 对高渗层内水平井产能进行预测时, 该模型预测结果较现有的3个模型更加精确。在实例计算过程中, 该模型所得产量预测值偏高20.54%, 而其他模型计算结果偏高达125.49%。

◇ 该模型采用闭合式公式, 使用非常方便。

◇ 该模型将井筒压力和/或流量剖面与相应的测量剖面拟合来对储层特性进行评估。

符号说明

C__根据公式A.25定义的常数;

c__公式A.15定义的常数;

dh__井眼等效直径, in;

ff__范氏摩擦系数, 无因次;

Jsp__比采油指数, stb/ (d5psi5ft) ;

L__井筒长度, ft ;

pr__油层压力, psi;

pwH__井底压力, psi;

Qoc__临界产量, stb/d;

QoH__井的总产量, stb/d;

xc__到水平段趾部的临界距离, ft;

yb__从井筒到边界的距离, ft;

μo__油的黏度, mPa5s;

导流弯管数值模拟分析 篇6

在实际弯管中,用于改善管内流体流动的方法之一是在所述弯管的内壁两侧对称布置两导流片,其中导流片的型线与弯管中心线在同一曲面内; 另外,方法二是沿弯管外侧内壁面设置一弯曲的折流片以强制干扰流体流动,以缩短流速达到平衡的出口直管段长度。

1计算模型

采用空气作为流动介质,取入口速度为20 m/s,弯管模拟计算区域为: 弯管入口直径为120 mm,入口直管段长度为150 mm,出口直管段长度为300 mm; 以此建立模型并计算。 首先,在GAMBIT中生成相应的弯管三维数值模拟模型。该三维模型为包括弯管及折流片在内的计算区域。其次,在GAMBIT三维建模软件中对建立的模型进行网格划分及边界条件设置。最后,输出FLUENT求解器能够计算的Mesh网格文件。将已建立好的计算流场三维物理模型和网格导入FLUENT6. 3中,对读入的网格质量进行检查并将FLUENT求解器的计算单位修改为米。FLUENT求解器选用基于压力的分离式求解器( Pressure Based) 。对于气体的流动,采用有限体积方法来离散控制方程。压力与速度之间的耦合通过Simple算法来处理。用可实现k - ε 双方程湍流模型来处理气体的湍流流动, 该模型能够有效的处理包括旋转均匀剪切流、包含有射流和混合流的自由流动、管道内流动等。流场计算方法采用SIMPLEC算法,对于基本方程离散差分主要采用如下格式: 压力项方程采用Standard格式,动量方程、湍动能方程和湍动能耗散率方程均采用二阶迎风格式。求解压力项时松弛系数为0. 3,密度项为1,体积力项为1,动量项为0. 7,湍动能项为0. 8,湍动能耗散率项为0. 8。所有CFD问题都需要有边界条件,对于瞬态问题还需要有初始条件,本计算模型中弯管入口采用速度入口边界条件,出口采用压力出口边界条件。

2计算结果分析

2. 1速度云图分析

图1 ~ 图3依次为弯管、带双导流片及带折流片的弯管流速分布图。从图1可看出,在入口直管段外径弧面处流体的流速减小,而在内径弧面处流体的流速相应的增大; 流体流经弯曲段过渡到出口直管段时,由于弯管段流体的惯性而流向外径弧面,造成流体层从内径至外径的分离从而引起出口流速不均匀。从图2可看出,由于弯管处两导流片的分流作用,不仅流体在弯管出口处的流速较均衡,流体在整个弯管内的流速变化都不大。从图3可看出,由于折流片对流体的强制转向作用, 在一定程度上改善了流体出口处的流速均衡性,但使得折流片后部管段的流体出现了新的分层。因此,双导流片能较好的调节弯管内的流速,使整体流速处于一较均衡的状态。

2. 2出口截面速度分析

图4速度曲线图的横坐标表示出口截面上的直径坐标值, 单位为m; 纵坐标表示出口截面上直径方向的速度变化曲线, 单位为m/s。从三条速度变化曲线可看出,无任何导流措施的

弯管出口截面最大速度为24 m/s,发生在出口外径侧,从出口截面外径至内径方向速度递减至14 m/s; 带双导流片的弯管出口截面上的速度较均衡,速度大小基本保持在20 m/s,即恢复至入口速度大小; 带折流片的弯管出口截面上的速度最大为22 m / s,最小值为15 m / s。三条速度曲线相比较而言,带双导流片的弯管出口截面上的速度最均衡且恢复至入口速度大小, 带折流片的弯管出口截面速度波动小于无任何导流措施的弯管。因此,双导流片对弯管内流体的流动特性改善效果明显, 强制折流片对弯管内流体的流动特性具有一定的改善效果。

3结论

河道工程导流施工技术探讨 篇7

关键词：围堰,河道工程,施工技术

1 围堰选型及布置

围堰断面型式选择需综合考虑河道特征水位、围堰使用期限及使用时段、使用功能、河道现状断面条件及地质情况等因素。

1.1 城市河道工程施工导流围堰选择的类型

1) 河道断面较宽, 河道净宽大于50 m时, 河道两岸均有建筑物施工需要时, 可选择在河道两岸分别设置纵向围堰。围堰断面如图1所示。

2) 河道断面较窄, 河道净宽小于30 m时, 可选择在河道中部设置一条纵向围堰, 采用分期导流的方式进行施工。围堰断面如图2所示。

3) 河道断面较窄时, 即河道净宽小于15 m的城区河道, 可选择小断面的袋装土围堰, 进行分期导流的方式分别进行左右岸施工。围堰断面如图3所示。

1.2 围堰的选型及布置主要考虑因素分析

1) 就地取材。水利工程施工围堰常常有混凝土围堰、土石围堰、袋装土围堰等型式, 由于城市河道工程单位长度内施工工期较短, 往往不选择混凝土围堰。同时, 考虑土方平衡及提高两岸永久建筑物基坑土方的利用率, 较小工程造价, 通常采用土围堰。

2) 场内外施工条件。城市河道两岸堤顶经常存在较密集民宅、工业厂房等建筑物, 或有较复杂的市政道路等设施, 两岸巡河路常常并未贯通, 因此将导致河道工程施工过程中场内外交通疏解难度大。在这种情况下, 可考虑利用临时围堰兼做场内临时施工道路的原则, 加宽围堰堰顶宽度, 并增加堰顶硬化即可解决场内施工道路问题, 以此解决工程施工交通部署。

3) 工程地质条件。根据工程本身不同的地形地质条件, 考虑不同的围堰断面型式, 遇到特殊的地质地形条件时针对性的考虑。在城市河道工程中较常见的异常地质条件是河床较厚淤泥层, 当河床淤泥层较厚时, 宜考虑选择围堰结构稳定性强、断面结构尺寸相对较小的围堰型式, 例如钢板桩围堰、增加钢管骨架的袋装土围堰等。

2 围堰表面处理及防渗

在城市河道工程中导流围堰工程往往采用土围堰, 因此在水土保持、防止扬尘等环境保护方面将有较高的要求, 同样在遇到透水层地质或采用渗透系数较大的堰体材料时需考虑增加防渗措施, 确保围堰安全。

2.1 围堰顶面处理

对于有兼做场内临时施工道路要求的围堰堰顶需进行硬化处理, 通常采用造价合理的泥结石铺筑, 同样在遇到不可预见的高水位时过水及防冲刷功能。对于枯水期围堰设计不需具备过水要求, 或不需兼做临时施工道路时, 通常采用沙袋或土工膜等进行临时防护, 防止水土流失。围堰迎水侧通常采用防渗土工膜+袋装土护面进行表面处理, 起到堰体防渗及防止堰体水土流失的作用, 特别对例如袋装土围堰堰体填充材料本身孔隙率较大时, 围堰迎水面防渗土工膜将发挥巨大的防渗作用。

2.2 围堰基础防渗处理

对于围堰基础以下河床存在砂砾层地质情况, 并且砂砾层高程高于基坑底部高程时, 通常需考虑围堰基础防渗处理。一般采用高压旋喷灌浆或自重灌浆等灌浆处理, 防止围堰基础砂砾层形成过水通道存在安全隐患。

3 纵横围堰布置及错车

河道工程中围堰工程量主要是纵向围堰, 合理适宜的横向围堰布置能有效地将建筑物基坑进行合理的分区, 通过横向围堰的分区可以最大限度地提高建筑物基坑开挖土方的重复利用率, 最大限度地减少外借土方工程量, 以减少工程造价, 并加快工程施工进度。对于兼做临时施工道路的围堰顶宽一般为4 m, 同时可结合横向围堰布置错车平台, 横向围堰及错车平台布置间距一般100~200 m。横向围堰及错车平台如图4所示。

4 永久建筑物导流

河道工程中新建截污系统通常采用箱涵或管道截流方式, 截流建筑物往往沿河道纵向布置, 在利用围堰导流施工完成一侧截污建筑物以后, 有河道另一侧截污建筑物施工要求或后期河道干地施工需要时, 可考虑利用已完建建筑物进行施工导流。在具备条件的河道上下游填筑全断面横向围堰利用永久建筑物导流时, 可以创造全河段干地施工条件。利用永久建筑物导流应具备以下条件。

1) 完建建筑物具备导流过水断面要求, 利用建筑物导流时对永久建筑物本身及周边相关构筑物排水不会带来影响。

2) 利用建筑物导流可以带来原河道干地施工条件, 为建筑物完建后围堰拆除、河道清淤等工程施工带来极大的便利。

3) 保证全河段干地条件后, 对河道原本的生态环境带来较小的影响, 或者通过全河段干地施工为河道岸坡恢复工程创造更好施工条件, 已达到更好的岸坡恢复效果。

5 支流及排放口导流

在河道工程中往往存在河道支流及两岸雨水、污水管道导流的问题, 当支流河道流量较大、水位较深时, 需考虑采用分期导流的方式, 对于河道两侧雨污管道可选择采用一次拦断架设管道导流或者直接采用抽水泵抽排的方式进行导流, 在永久建筑物施工完毕具备通水条件后拆除围堰直排其中。

6 结语