本发明公开了一种水力压裂裂缝延伸与多相流体流动模拟方法，包括以下步骤：(1)收集工况和输入参数；(2)建立两相流应力平衡方程；(3)建立两相流体流动控制方程；(4)建立两相流相场演化方程；(5)采用有限元数值离散方法和Newton‑Raphson(NR)迭代法建立上述方程组的数值求解迭代格式；(6)将步骤(1)中的参数带入步骤(5)建立的求解迭代格式中模拟不同工况下的水力裂缝延伸轨迹。与现有的水力裂缝延伸模拟方法相比，考虑了压裂过程中压裂液从裂缝中滤失到地层后的两相流体流动特性，为水力压裂裂缝延伸轨迹的预测提供更准确的方法。