压裂设计作业

人工裂缝支撑剂铺置优化设计

讲义

1 课程设计相关要求

1.1 专业硕士课程设计目的

目的：掌握水力压裂过程中支撑剂在裂缝中的铺置规律，并对压裂施工参数进行优化设计。

完成途径：理论学习、数学计算、软件应用

1.2 主要任务

每3人一组完成一个课题，每个课题参数不同。

主要任务：

（1）计算砂粒沉降匀速、阻力速度、平衡流速、平衡时断面高度、平衡时间、平衡砂堤高度；

（2）敏感性分析；

（3）参数优化设计

注意：计算参数由同学自行设计，可参考示例。同一组同学参数可一样，不同组同学参数必须完全不相同，否则两组同学成绩均为零。

2 支撑剂沉降数学模型

支撑剂在裂缝中的分布情况，决定了压裂后填砂裂缝的导流能力和增产效果。而支撑剂在裂缝中的沉降受到诸多因素的影响，包括砂浓度、壁颗粒形状等因素。

2.1 支撑剂在压裂液中的自由沉降

对于支撑剂在压裂液中的沉降现象，研究的起步是非常早的。早在50年代，人们就认识到了支撑剂在输送过程中，由于液体的粘性力不可能完全的大于支撑剂颗粒自身的重力，会引起支撑剂在随着压裂液的前进过程中产生沉降现象[7]；并对这种现象进行了理论上的研究，应用于水力压裂的设计中，取得了一定的效果。目前对于在牛顿液体中单颗粒的自由沉降或群粒的干扰沉降均有比较成熟的计算方法，对于非牛顿液体的幂律型液体中的沉降规律也有所认识、但对于带有粘弹性的冻胶压裂液中的沉降行为及其计算方法，只能说处于定性的、半定量的研究阶段，还有许多工作等待继续去做[8]。 2.1.1 单颗粒在牛顿流体中的自由沉降

质量为m 的颗粒在力1F 的作用下，在液体中以速度u 沉降。若2F 及3F 分别代表浮力及阻力，则在速度u 的方向上合力F 为：

123F F F F =-- (2-1)

将以上各种力的参数代入上式，则可写成下式：

2

(/)()2

s d A u F ma ma C ρρρ=-- (2-2)

式中 a ——加速度，m/s 2；

ρ——压裂液净液的密度，kg/m 3； ρs ——砂粒的密度，kg/m 3；

A ——垂直于沉降方向的颗粒面积，m 2； C d ——阻力系数，无因次； m ——颗粒的质量，kg ；

u ——单颗粒的重力沉降速度，m/s 。

设球形颗粒在重力下沉降，则上式中的a ，A ，m 分别为：

a g =，24p A d π=，36

p s m d π

ρ=。

式中 g ——重力加速度，m/s 2；

d p ——颗粒直径，m 。

又知：du

F ma m dt

==

将上面各式分别代入(2-2)式中并进行整理，可得到：

22/42

s d p du u m mg mg C d dt πρρρ=-- (2-3) 2

()34s d s p s

du g C u dt d ρρρρρ-=-

(2-4) 颗粒在自由沉降中先是不断地加速，其受到的液体阻力也在增大，当作用力与阻力达到平衡时，砂粒以均匀的速度下沉，即0du

dt

=。此时，可得到单个颗粒在牛顿流体中的沉降速度u=u p ，为：

1/2

4()3s p p d g d u C ρρρ⎡⎤-=⎢⎥⎣⎦

(2-5)

当使用(2-5)式求解颗粒的自由沉降速度时，首先应该确定阻力系数C d ，颗粒在沉降过程中，带动周围液体使之也产生运动[9]。当液体运动速度在层流范围内时，斯托克斯在1851年给出了线性纳维尔斯托克斯方程的解。牛顿液体的本构方程为：

τμw

du

dr

=() (2-6) 式中的τw 及du dr

分别为：

τπw p

F d =

3

2

及du dr u d p =3 代入(2-6)解出F 3为：

F d u p 33=πμ (2-7)

式(2-7)的F 3与式(2-2)的F 3是恒等的，即：

34222πμπ

ρd u C d u p d p

=⎛⎝

⎫

⎭

⎪ (2-8) 所以：

Re

2424

d p C d u N ρμ

=

= (2-9) 式中 N Re ——雷诺数，无因次；

μ——液体的粘度，mPa·s 。

上式中的阻力系数C d 只满足层流范围。作C d 与N Re 的关系曲线图，可以看出，式(2-9)只能在N Re <1时使用。将式(2-9)代入式(2-5)，得到斯托克斯的沉降等式：

2

()18s l p

p g d u ρρμ

-= (2-10)

近年来由于压裂工艺的发展，许多学者引用了各种条件下计算单颗粒自由沉降的计算公式，其中诺沃特尼[10]建议如下：

(1)当N Re ≤2时， Re

24

d C N =

2

()18s l p

p g d u ρρμ

-= (2-11)

(2)当2

Re 18.5

d C N =

0.71 1.44

0.29

0.43

20.34()s p p d u ρρρ

μ

-=

(2-12)

(3)当N Re ≥500时， 0.44d C =

p u = (2-13)

根据希勒与诺曼的建议[11]

，当雷诺数从0.2化到500-1000时，可以用下列通用方式求阻力系数C d ：

10.687

Re Re 24(10.15)d C N N -=+ (2-14)

由式(2-13)也可导出颗粒自由降落的阻力与雷诺数的关系：

(1)1Re 24d C N -=时，F d u 3315=.πμρ

(2)1Re 24d C N -=()10150687+.Re .N 时，F d u 3315=.πμρ()

1015

0687+.Re .N (3)d C =0.44时，F d u p 3220055=.πρ (4)d C =0.10时，F d u p 32200125=.πρ

2.1.2 单颗粒砂粒在幂律流体中的自由沉降

目前，在水力压裂中多数情况下使用水基植物冻胶或高分子聚合物配制的压裂液，这种压裂液的性能至今还没有很全面的认识，但在一定一条件下一般把它看作幂律流体来处理

[12]

。对于牛顿液体来说，作用在粒子上的剪切速率并不影响支撑剂的沉降速度，这是由

于牛顿液体的粘度相对于剪切速率来说是常数。然而对于幂律流体而言，由于液体比不流动时的剪切力大，支撑剂的沉降就较快。剪切速率对支撑剂的影响可用综合的剪切速率来表示如果幂律流体的流变指数与稠度系数分别为n 和k ，则其雷诺数N Re 可表示如下：

Re p p

a

d u N ρμ=

(2-15) 式中 N Re ——幂律流体的雷诺数，无因次；

ρ——压裂液的密度，kg/m 3；