Abstract. The study is devoted to modeling of multiphase flows of immiscible viscoplastic fluids in a hydraulic fracture. In the framework of the lubrication approximation, threedimensional Navier-Stokes equations are reduced to hyperbolic transport equations for the fluid tracers and a quasi-linear elliptic equation in terms of the fluid pressure. The governing equations are solved numerically using the finite-difference approach. A parametric study of the displacement of Bingham fluids in a Hele-Shaw cell is carried out. It is found that fingers developed through the pillar of a yield-stress suspension trigger the development of unyielded zones. An increase in the Bingham number leads to an increase in the so-called finger shielding effect, which manifests itself via an increase in the overall finger penetration zone and a decrease in the total number of fingers. The effect of flow parameters on the displacement of hydraulic fracturing proppant-laden suspension by a clean fluid in the vicinity of the perforation zone is carried out. This particular case is considered in application to overflush at the end of a stimulation treatment, when a small portion of a thin clean fluid is injected to wash out the particles from the wellbore into the fracture. It is found that an increase in the yield stress and the viscosity contrast between the fracturing and the overflush fluids typically reduces the area of the cavity thus mitigating the risk of loosing the conductive path between the wellbore and the fracture after the fracture closure.
IntroductionProppant transport models incorporated into existing hydraulic fracturing simulators describe the flow of particle-laden suspension inside a hydraulic fracture in the framework of the lubrication approximation using the power-law rheological model (see [1]). Hydraulic fracturing suspensions with large concentration of solids or fibers show a yield-stress behavior in rheological experiments [2,3], which is not taken into account in the proppant transport models implemented into commercial simulators of hydraulic fracturing.A state of the art in the modeling of injection of particle-laden suspensions into hydraulic fractures is the family of 2D width-averaged models based on the lubrication approximation to Navier-Stokes equations [1,4,5,6]. In the regime of non-inertial settling, the momentum conservation equation for particles is reduced to an algebraic relation for the particle velocity slip in the vertical direction given by the Stokes formula with a correction for hindered-settling effects due to a finite particle volume fraction. In the case of Newtonian suspension rheology, the total momentum conservation equation for the suspension is reduced to the linear expression of the fluid (or mixture) velocity through the pressure gradient (similar to the Darcy law, hence a well-known analogy between filtration and a flow in a Hele-Shaw cell). The dependence of the