何右安，高武彬，路 敏

(1.长庆油田分公司勘探开发研究院，陕西 西安 710018；2.低渗透油气田勘探开发国家工程重点实验室，陕西 西安 710018)

采取水平井开发和大规模水力压力投产是非常规油气资源开发的重要手段[1-5]，然而水平井多级压裂过程中的裂缝扩展容易受到地应力、裂缝间的相互作用和施工参数等因素的影响[6-8]，因此，预测和模拟多簇压裂裂缝扩展动态规律对确保水平井多级压裂效果十分重要。经典裂缝理论模型包括KGD模型[9]和PKN糢型[10]，在此基础上学者们提出了拟三维模型[11-13]和全三维模型[14]，但这些模型主要用于描述单条裂缝扩展。此后，许多学者借助各种数值模拟方法如Cohesive单元法、扩展有限元法和离散元法等[15-18]来描述非常规储层中多裂缝同时扩展的复杂形态，这些研究模拟的裂缝多为二维形态，并且将注入地层中的压裂液视为纯液体，没有考虑实际压裂液中有支撑剂的存在，导致模拟结果与现场实际结果存在较大差异，难以直接应用于现场优化设计。三维位移不连续法[19]能够很好地描述岩体的三维变形，并且只需要对边界区域(即裂缝面)进行剖分，可极大减少计算量，同时省去了对远离裂缝区域的多余计算，节省大量计算空间，非常适用于描述致密非常规储层的三维裂缝扩展形态。因此，基于三维位移不连续法，分别考虑压裂液、支撑剂在井筒和水力裂缝中的流动，结合J井建立了三维分段多簇压裂数值模型，采用Newton-Raphson法求解数值模型，并分析了压裂液排量、黏度、砂比和簇间距等工程因素对多裂缝扩展的影响。……