水力压裂力学.pptxVIP

会计学;水力压裂力学是对压裂工艺和压裂机理的简单描述。;开发和利用水力压裂施工的重要原因;6.2 早期水力压裂模拟;对于缝高hf不变和无限大（即平面应变）裂缝其最大宽度为： ;; 6.2.2 水力压裂二维模拟; KGD模型 ;6.2.2.1 垂向裂缝的Perkins 和 Ken模型的推导; 实际的裂缝宽度：;6.2.2.2 模型中考虑流体滤失;假设裂缝在空间和时间上都保持恒定，上式变为：;Nordgren（1972）;压力用缝宽表示代替，方程（6.18）写为：;6.2.2.3 Khristianovich-geerssma-de Klerk 模型的导出; Geertsma 和 de Klerk（1969）对于缝端区域很小这个问题给出了解。; 裂缝宽度方程：; 井壁裂缝宽度：; 假设流体滤失对裂缝形态或压力分布没有影响，将模型推广到包括流体滤失的情况下：;6.2.2.4 PKN 和 KGD 模型的假设;前面简单模型的局限性： ? 需要给定缝高或假设产生的是径向缝;包括缝高增长的三种主要水力压裂模型; 平面三维模型 ; 拟三维模型 ; 6.3.1 平面三维模型;Hirth 和Lother（1968）以及Bui（1977）;简单的裂缝形态和缝内压力分布如定压下的椭圆形裂缝，破裂准则：方程（6.3）;6E 拟三维模型中的水平耦合; 质量守恒方程（描述流体流动）：; 方程（6.34）的左边为动量改变速率；右边分别为压力、粘滞力和重力，它可解释为小的流体单元在力的作用下而加速。该方程可以扩展并根据压裂地层的不同形状而简化（见旁注6F）。;应力与流速之间的本构方程：;6F 水力压裂中的动量守恒; 对于非渗透介质，z方向分量可以忽略不计；并假设为稳态流，（6F.3）简化为：;牛顿流体; 对于不可压缩流体，正应力方程中的第三项为零;将方程（6F.6）代入（6F.4）得到：;对上式积分，得整个通道上的平均流速：;6G 非牛顿流体的动量平衡和本构方程; μa可能仅是?的函数，而且两者之间存在某种函数关系：; 稠度指数K’取决于流形，并与流体基本性质——综合稠度系数K有关。对于平板流动，关系为：;; 方程（6.31）到（6.38）总结了牛顿流体的平面模型，非牛顿流体也可得到相似的结果（见旁注6G）。; 平面模型的数值模拟; Barree给出的另外一种数值??拟方法;; 6.3.2 以单元为基础的拟三维模型; 这个假设使得可将固体力学和断裂力学解与流体流动分开：;6.3.2.1 固体力学解;（6.40）;式中： pf —— 流体密度 hcp —— 射孔中心高度 hi —— 由底部缝端至第i 层底层顶部的高度，如图 6.7所示; 通过迭代求解该非线性方程，如果pcp取某一值时的解（两个垂向缝端的位置以及压力）已知，假设缝高增量，则两个垂向的缝高就可以计算出来而且满足方程（6.40）和（6.41），裂缝端达到这些位置所需的pcp可计算出来。最后与该解相关的缝宽剖面也可得到。;现如考虑如图6.8所示的一个对称的三个地层的情况;图6.9为由（6.43）计算得到的缝高和净压力的关系曲线;6.3.2.2 流体力学解;采用平均流速和宽度（宽度以横截面积除以高度代替） 质量守恒方程简化为：;对于具有性质n和K的幂律流体：;求沿缝高方向的所有横截面的总流速 ? 将方程（6.45）中换成平均流速 ? 将方程（6.48）从横截面的缝底到缝顶积分;6.3.2.2.1 层流和紊流;6.3.2.2.2 压裂液的流变性;6.3.2.3 模型的数值解; 显式有限差分方法;将导数以中心有限差分近似代替得到：; 这就是拉格朗日动坐标方法，它的局限性在于： ? 网格点坐标随流体移动，滤失引起每个单元收缩甚至随着深入裂缝中而消失。 ? 井筒必须有新的单元不断加入，这就使得控制任意时刻运用多少单元或单元的大小很困难。;6H.1 坐标系拉伸; 假设压力梯度为：; 为了进行误差分析，所有变量必须为绝对值，D定义为：; 为将A用A（1+ε）代替， A（1+ε）可近似为（对于ε小量）：; 显式有限差分的一个主要局限：在计算中所用的时间步长不能