A novel multiphase and multicomponent model for simulating molecular diffusion in shale oil reservoirs with complex fracture networks

Molecular diffusion is critical for enhanced oil recovery (EOR) in shale oil reservoirs with complex fracture networks. Understanding the influence of fractures on diffusive mass transfer is crucial for predicting oil recovery and remaining oil distribution. Diffusive mass transfer between fractures and matrix is critical in comprehensively and effectively simulating molecular diffusion. Resolution of matrix cells significantly affects diffusion accuracy at the fracture–matrix interface. Low resolution results in multiple fractures in the same matrix cell, leading to decreased precision in calculating mass transfer by conventional methods. To address this, a novel multiphase and multicomponent model is proposed. The new model integrating the consideration of fracture spacing modifies molecular diffusion transmissibility between fracture and matrix in an embedded discrete fracture model. The discretization employs the two-point flux approximation in the finite-volume method. Validation compares the coarser mesh to the finest grid as a reliable reference. Results show the proposed model accurately captures diffusive mass transfer in a coarser mesh. Modified models study molecular diffusion's effects on EOR in shale oil reservoirs with complex fracture networks by CO2 huff and puff. Results indicate that increasing injection rates cannot improve oil recovery under extremely low porosity and permeability. Molecular diffusion facilitates CO2 penetration into the formation. This expands the swept CO2 volume and increases both volume expansion and formation energy. In addition, the light and heavy components of the crude oil are diffused into the fractures and eventually produced, which reduces gas production in the case of diffusion.