目前,我国油气资源开发逐渐从常规油气为主过渡到非常规油气为主,非常规油气储层渗透性低,开采难度大,需要进行储层改造以实现增产。而作为主要储层改造方法的水力压裂技术虽然已经成功应用70余年,但开采量离预期产量仍有很大差距,因此有效地模拟水力压裂是压裂设计及产量评估的关键。水力压裂是一个流固耦合过程,同时储层中含有的大量节理裂隙也为数值模拟增加了难度。为此本文在虚内键(VIB)理论基础上发展了一套水力压裂数值模拟方法,开展复杂储层水力压裂模拟研究,探讨水力裂纹扩展规律和水力缝网形成的影响因素,探索脉冲压裂下的裂纹扩展规律,为非常规油气资源的开采提供必要的理论与技术储备。为了模拟水力压裂中渗流场与应力场的耦合作用,首先建立了裂隙介质渗流模型,将离散虚内键(DVIB)模型扩展到裂隙介质渗流模拟。在该渗流模型中,采用键元胞来表征代表单元体内的微裂隙网络,宏观裂隙的渗流特性则由劈裂元胞(被裂缝穿过的键元胞)表征,元胞内每条键视为一条水力通道,将三维渗流问题转化为一维渗流。基于理想键元胞思想推导了“键”的等效渗透率和储水率与宏观渗流参数的理论关系。模型可以模拟含任意条裂缝的介质渗流问题,避免将裂缝设置为内边界和网格重构。模拟结果验证表明DVIB渗流模型可以很好地模拟复杂储层的渗流过程。建立了流固耦合的DVIB模型,采用同一套格子计算岩体力学变形和渗流过程,在微观“键”层面上统一了力学场和渗流场。通过建立劈裂元胞的水力耦合关系及接触本构关系,可以有效地反映裂缝面间的基本摩擦机制。这使得该模型能模拟水力裂缝与天然裂缝的交汇作用及复杂储层中的水力裂纹扩展,避免网格重构。为了模拟大尺度压裂过程,在改进劈裂单元方法基础上建立了全耦合水力压裂模型。在劈裂元胞中设置虚节点(裂纹与元胞边界交点),借助虚节点推导了劈裂元胞的渗流矩阵、力学矩阵及流固耦合矩阵。虚节点的水压和位移是多余自由度,为消除这些多余自由度,认为虚节点与其邻近且在裂纹同侧的实节点相关,采用最小二乘法技术将二者联系起来,从而消除了多余自由度。这样既允许裂纹从单元中穿过,又不引入额外自由度,为大尺度复杂储层压裂模拟提供了一条有效途径。全耦合模型同时考虑了流体的滤失和裂缝填充效应。与KGD模型和解析解对比验证表明该方法具有很高精度。对射孔转向模拟中发现:射孔角度对裂纹转向具有重要影响。地应力差相同时,射孔倾角越大,裂纹转向越明显。同时,地应力差越大,裂纹转向越快。脉冲压裂是当前尚在探索而远未成熟的一项压裂方法。通过脉冲压裂可以形成复杂缝网结构。为了探索脉冲压裂规律,本文采用DVIB模型模拟了脉冲压裂过程。发现:加压速率、加压时间、峰值压力、卸压速率共同影响了井筒上裂纹的起裂数量和裂纹扩展后的分叉程度。峰值压力的大小控制了裂纹分叉的程度;通过控制压力水平,可以调整多裂纹的形成位置;加压速率一定时,卸压速率越小,缝网结构越复杂,压裂区域越大,压裂效果就越好。另外,最大地应力控制了裂纹扩展的方向,地应力差越大,最大地应力方向的裂纹分叉越密集。射孔角度越大,裂纹转向越明显。大尺度下脉冲压裂裂纹扩展规律与小尺度下基本一致。在评估或设计脉冲压裂时,应综合考虑上述影响因素才能更准确地评估压裂效果。本文为复杂地层水力压裂模拟提供了一条有效的方法。