压裂设计作业

1人工裂缝支撑剂铺置优化设计讲义21课程设计相关要求1.1专业硕士课程设计目的目的：掌握水力压裂过程中支撑剂在裂缝中的铺置规律，并对压裂施工参数进行优化设计。完成途径：理论学习、数学计算、软件应用1.2主要任务每3人一组完成一个课题，每个课题参数不同。主要任务：（1）计算砂粒沉降匀速、阻力速度、平衡流速、平衡时断面高度、平衡时间、平衡砂堤高度；（2）敏感性分析；（3）参数优化设计注意：计算参数由同学自行设计，可参考示例。同一组同学参数可一样，不同组同学参数必须完全不相同，否则两组同学成绩均为零。32支撑剂沉降数学模型支撑剂在裂缝中的分布情况，决定了压裂后填砂裂缝的导流能力和增产效果。而支撑剂在裂缝中的沉降受到诸多因素的影响，包括砂浓度、壁颗粒形状等因素。2.1支撑剂在压裂液中的自由沉降对于支撑剂在压裂液中的沉降现象，研究的起步是非常早的。早在50年代，人们就认识到了支撑剂在输送过程中，由于液体的粘性力不可能完全的大于支撑剂颗粒自身的重力，会引起支撑剂在随着压裂液的前进过程中产生沉降现象[7]；并对这种现象进行了理论上的研究，应用于水力压裂的设计中，取得了一定的效果。目前对于在牛顿液体中单颗粒的自由沉降或群粒的干扰沉降均有比较成熟的计算方法，对于非牛顿液体的幂律型液体中的沉降规律也有所认识、但对于带有粘弹性的冻胶压裂液中的沉降行为及其计算方法，只能说处于定性的、半定量的研究阶段，还有许多工作等待继续去做[8]。2.1.1单颗粒在牛顿流体中的自由沉降质量为m的颗粒在力1F的作用下，在液体中以速度u沉降。若2F及3F分别代表浮力及阻力，则在速度u的方向上合力F为：123FFFF(2-1)将以上各种力的参数代入上式，则可写成下式：2(/)()2sdAuFmamaC(2-2)式中a——加速度，m/s2；ρ——压裂液净液的密度，kg/m3；ρs——砂粒的密度，kg/m3；A——垂直于沉降方向的颗粒面积，m2；Cd——阻力系数，无因次；m——颗粒的质量，kg；u——单颗粒的重力沉降速度，m/s。设球形颗粒在重力下沉降，则上式中的a，A，m分别为：ag，24pAd，36psmd。式中g——重力加速度，m/s2；dp——颗粒直径，m。又知：duFmamdt将上面各式分别代入(2-2)式中并进行整理，可得到：422/42sdpduummgmgCddt(2-3)2()34sdspsdugCudtd(2-4)颗粒在自由沉降中先是不断地加速，其受到的液体阻力也在增大，当作用力与阻力达到平衡时，砂粒以均匀的速度下沉，即0dudt。此时，可得到单个颗粒在牛顿流体中的沉降速度u=up，为：1/24()3sppdgduC(2-5)当使用(2-5)式求解颗粒的自由沉降速度时，首先应该确定阻力系数Cd，颗粒在沉降过程中，带动周围液体使之也产生运动[9]。当液体运动速度在层流范围内时，斯托克斯在1851年给出了线性纳维尔斯托克斯方程的解。牛顿液体的本构方程为：wdudr()(2-6)式中的w及dudr分别为：wpFd32及dudrudp3代入(2-6)解出F3为：Fdup33(2-7)式(2-7)的F3与式(2-2)的F3是恒等的，即：34222duCdupdp(2-8)所以：Re2424dpCduN(2-9)式中NRe——雷诺数，无因次；μ——液体的粘度，mPa·s。上式中的阻力系数Cd只满足层流范围。作Cd与NRe的关系曲线图，可以看出，式(2-9)只能在NRe1时使用。将式(2-9)代入式(2-5)，得到斯托克斯的沉降等式：2()18slppgdu(2-10)5近年来由于压裂工艺的发展，许多学者引用了各种条件下计算单颗粒自由沉降的计算公式，其中诺沃特尼[10]建议如下：(1)当NRe≤2时，Re24dCN2()18slppgdu(2-11)(2)当2NRe500时，0.6Re18.5dCN0.711.440.290.4320.34()sppdu(2-12)(3)当NRe≥500时，0.44dC()1.74sppgdu(2-13)根据希勒与诺曼的建议[11]，当雷诺数从0.2化到500-1000时，可以用下列通用方式求阻力系数Cd：10.687ReRe24(10.15)dCNN(2-14)由式(2-13)也可导出颗粒自由降落的阻力与雷诺数的关系：(1)1Re24dCN时，Fdu3315.(2)1Re24dCN10150687.Re.N时，Fdu3315.10150687.Re.N(3)dC=0.44时，Fdup3220055.(4)dC=0.10时，Fdup32200125.2.1.2单颗粒砂粒在幂律流体中的自由沉降目前，在水力压裂中多数情况下使用水基植物冻胶或高分子聚合物配制的压裂液，这种压裂液的性能至今还没有很全面的认识，但在一定一条件下一般把它看作幂律流体来处理[12]。对于牛顿液体来说，作用在粒子上的剪切速率并不影响支撑剂的沉降速度，这是由于牛顿液体的粘度相对于剪切速率来说是常数。然而对于幂律流体而言，由于液体比不流动时的剪切力大，支撑剂的沉降就较快。剪切速率对支撑剂的影响可用综合的剪切速率来表示如果幂律流体的流变指数与稠度系数分别为n和k，则其雷诺数NRe可表示如下：ReppaduN(2-15)式中NRe——幂律流体的雷诺数，无因次；ρ——压裂液的密度，kg/m3；6dp——支撑剂粒径，m；up——支撑剂在幂律流体中的沉降速度，m/s；μa——幂律流体的视粘度，Pa·s；丹尼什[13]取颗粒在匀速沉降时的剪切速率为3ppud，则对于幂律流体而言，其视粘度为：1naaK(2-16)其中21()3nanKKn式中Ka——幂律流体在裂缝中流动时的稠度系数，Pa·sn；γ——压裂液流动剪切速率，s-1；n——压裂液的流变指数(反映偏离牛顿流体的程度)，无因次。将式(2-16)代入式(2-15)得：2Re13nnppnaduNK(2-17)式中K——压裂液的稠度系数，Pa·sn；利用当ReN1时，Re24dCN，1/24()3sppdgduC及式(2-17)得到：11n1()18(3)nnSppngduK(2-18)因为在裂缝中流动时：21()3nanKKn代入式(2-18)，整理后有：1(21)()90006000npsppanddunK(2-19)式中各变量如以上各式。诺沃尼特采用的剪切速率不同于丹尼什，Dudpp，此时斯托克斯定律为：1()18npsppndgudK(2-20)7比较式(2-5)与(2-19)或式(2-20)可以看出，砂粒在牛顿液体中的沉降速度与粒径的平方成正比，而在幂律液体中与粒径的(1+n)/n成比例。当n=1时，二者相同，液体的非牛顿性愈强，差别也愈大。上述计算砂粒在幂律液体中沉降的方法，只考虑到由于砂粒下降时在液体上所产生的剪切作用。除此之外在流动的液体中，液体的流变性质还与流动所产生剪切速率有关。因此诺沃特尼认为在流动的幂律液体中由于颗粒的沉降所引起的总剪切速率，应该是二者之代数和，即：2222122ppuDDDDd(2-21)式中D——总剪切速率，s-1。1D——液体的流变性质所产生的剪切速率，s-1。2D——液体流动所产生的剪切速率，s-1。在层流范围内(斯托克斯定律)，砂粒的沉降匀速pu应为22211()()()181818pspspsnpnnnadgdgdguDKDK(2-22)将式(2-21)代入式(2-22)：122222()18npsppnpdguuDKd(2-23)上式要用试算法求解，式中的D与液体在缝中的流速分布有关。对幂律液体在缝中的流速分布为：12111/nnnyuunWZ(2-24)式中y——缝中心线至壁面的任意距离，cm；W——缝宽，cm；ū——缝中的平均流速，cm/s。式(2-24)对y求导数，1221()1/nduunydyWnWz(2-25)所以：12221()()1/nduunyDdyWnWz(2-26)在缝的中心线上(y=0)，D2=0，液体没有受到剪切如同在静止的液体中一样，砂粒在此处沉降得慢。在壁面附近，由于剪切速率的增加粘度降低，沉降速度应快一些。82.2颗粒沉降速度的影响因素压裂施工中，压裂液和支撑剂是按照一定的砂比注入裂缝的，支撑剂在这种条件下的沉降与单颗粒的自由沉降是不同的。多颗粒在沉降时会相互干扰，并且支撑剂在沉降过程中还受到裂缝壁面的影响[14]。这些影响因素必然会影响支撑剂的沉降速度。2.2.1砂浓度对沉降速度的影响压裂液中多颗粒沉降时，由于粒间相互的干扰作用，使得支撑剂的沉降速度低于单粒的自由沉降速度。这种相互干扰作用包括以下两个方面：(1)单颗粒的沉降引起周围液体的向上流动，阻碍了周围流体的下沉，砂比越高，阻尼作用越大；(2)混有支撑剂的液体混合物在比重、粘度上都有所增加，其结果是增加了支撑剂的浮力和沉降的阻力，这都使沉降速度变慢。通常以浓度校正系数fc反映干扰沉降的影响，fc的表达式如下：(1)mcpff(2-27)式中fc——砂浓度校正系数，无因次；fp——混砂液中支撑剂所占的体积百分数，无因次；m——与雷诺数相关的经验常数，无因次。有关文献建议的修正公式如下[15]：(1)当NRe≤2时，m=5.5；(2)当2NRe≤500时，m=3.5；(3)当NRe≥500时，m=2。由此可见，支撑剂的体积浓度对其沉降速度的影响是很大的。2.2.2裂缝壁面效应混砂液不是在无限大容器中自由沉降，而是在有限缝高的裂缝内沉降。存在的裂缝壁而将增大颗粒沉降阻尼从而降低沉降速度[16]。垂直裂缝中的壁面校正系数fw按下式计算：当NRe1时，34510.65260.1470.1310.0644ppppwddddf(2-28)当NRe100时，1.5=12pwdfW(2-29)若雷诺数介于1与100之间，则可进行简单的直线内插法取值。2.3单颗粒修正后的速度9混砂液沿着裂缝流动，由于裂缝横剖面尺寸、剪切速率和砂浓度都是变化的，导致支撑剂在裂缝中各点的沉降速度各不相同。最终沉降速度为：ppcwuuff(2-30)式中up'——修正后的支撑剂沉降速度，m/s；up——单颗粒匀速沉降速度，m/s。则不同雷诺数的沉降速度分别如下：(1)当NRe1时，34525.510.65260.1470.1310.0644118psppppppdgdddduf(2-31)(2)当1≤NRe≤2时，25.5118pspwpdguff(2-32)(3)当2≤NRe100时，0.711.443.50.290.4320.341spppwduff(2-33)(4)当100≤NRe≤500时，1.50.711.443.50.290.4320.34112sppppddufW(2-34)(5)当NRe500时，1.521.74112sppppgddufW(2-35)2.4支撑剂在裂缝中分布的计算方法水力压裂是在地层中造出一定导流能力及一定长度的填砂裂缝，以此改变地层中