Fine particle migration is a crucial issue in hydraulic fracturing-based production from unconventional reservoirs and can significantly affect the effective conductivity of the packed proppant layers. In this study, to investigate the conveying behavior of fine particles, a Eulerian-Lagrangian method, namely the computational fluid dynamics–discrete element method (CFD-DEM), is adopted for 3D simulations. Two types of particles are involved in this process: larger proppant particles and fine particles. Proppant particles form a background skeleton, and a two-phase flow of fine particles and fluid is constrained in the connected pore networks. We adopt a unified CFD-DEM framework to resolve the non-slip boundary condition on the proppant particle surface and track fine particles using a four-way coupling strategy. In particular, a forcing term based on the fictitious domain method is embedded in a regular volume-averaged Navier-Stokes equation to solve the fluid motion. A perfectly packed skeleton of a face-centered cube is designed for subsequent numerical studies. The results show that the fluid and particle flow velocities are higher near the fracture surfaces than inside the interlayers. In addition, it is found that as the fine particle concentration increases, the average horizontal velocity increases, whereas the settling velocity tends to decrease.
Read full abstract