Abstract
In low-permeability reservoirs, such as shale and tight sandstone, imbibition is an important mechanism for enhancing oil recovery. After hydraulic fracturing treatment, these reservoirs create a network of fracture pathways for fluid flow. Therefore, understanding the imbibition mechanisms in fractured porous media and quantitatively characterizing oil–water distribution are crucial for the development of low-permeability reservoirs. In this study, a mathematical model of two-phase flow in porous media with branching fractures was established. The phase-field method was employed to track the oil–water interface, and quantitative characterization of imbibition was conducted based on fractal theory, and the effects of wetting phase injection rate, the number of disconnected fractures, fracture spacing, and fracture morphology on imbibition in branched fracture porous media were discussed. The research findings indicate that in branched fracture porous media, both co-current and countercurrent imbibition processes occur simultaneously, and there exists a diffusion interface layer with a certain thickness at the oil–water interface. The hydraulic pressure generated by the wetting phase injection rate provides the driving force for imbibition oil recovery, but it also affects the contact time between the wetting and non-wetting phases. The presence of disconnected fractures hinders the propagation of hydraulic pressure, reducing the effectiveness of imbibition. The imbibition displacement zone is limited and occurs only within a certain range near the fractures. As the number of branching fractures increases, the channels for the wetting phase to enter matrix pores are enhanced, resulting in higher efficiency of imbibition displacement of the oil phase. The results of this research can provide guidance for the design of fracturing programs and recovery prediction in low-permeability reservoirs.