5. Conclusion
To summarize this work, we have shown that the lattice Boltzmann method can be accurately applied to the case of porous media flows at the REV scale within a fracture network environment. To do so, the original lattice Boltzmann method is extended through the Guo et al. method to incorporate damping terms present within the flow field due to the presence of porous media. A modified Zou–He boundary condition is then developed to model the interface between the well and fracture system with the surrounding porous media. The treatment of the boundary is important because the largest pressure drop exists immediately around the wellbore and fracture system and therefore the majority of the producible fluid is in this region. Lastly, a system of non-dimensional governing equations are derived for the case of the FRS such that the conversion between the lattice and physical system is consistent