Due to the coexistence of multiple types of reservoir bodies and widely distributed aquifer support in karst carbonate reservoirs, it remains a great challenge to understand the reservoir flow dynamics based on traditional capacitance–resistance (CRM) models and Darcy’s percolation theory. To solve this issue, an improved injector–producer-pair-based CRM model coupling the effect of active aquifer support was first developed and combined with the newly-developed Stochastic Simplex Approximate Gradient (StoSAG) optimization algorithm for accurate inter-well connectivity estimation in a waterflood operation. The improved CRM–StoSAG workflow was further applied for real-time production optimization to find the optimal water injection rate at each control step by maximizing the net present value of production. The case study conducted for a typical karst reservoir indicated that the proposed workflow can provide good insight into complex multi-phase flow behaviors in karst carbonate reservoirs. Low connectivity coefficient and time delay constant most likely refer to active aquifer support through a high-permeable flow channel. Moreover, the injector–producer pair may be interconnected by complex fissure zones when both the connectivity coefficient and time delay constant are relatively large.

The flow field is the fluid dynamic flow path in strongly heterogeneous reservoirs, and its evolution significantly affects the distribution of remaining oil, showing a disordered and highly dispersed mode caused by long-time water injection. By combining traditional flow-field evaluation with flow diagnostics, this paper proposes a methodology to quantitatively characterize and adjust the flow field in real time during the late stage of waterflooding in strongly heterogeneous reservoirs. In the study, the fluid velocity, abundance of predominant remaining oil, and Lorenz coefficient are preferred as the characteristic parameters to evaluate the effect of reservoir and flow heterogeneity on the flow field. Taking the minimization of the Lorentz coefficient as the objective function, the optimal injection and production parameters are obtained by dynamic adjusting the flow field. The results show that as water injection continues, the shape and variation of the flow field are jointly influenced by the reservoir rhythm, vertical permeability contrast, and lateral permeability distribution. The larger the permeability contrast, the greater the influence of the lateral permeability distribution. When the permeability contrast is large, the Lorenz coefficient strongly depends on the lateral permeability distribution. Finally, this method is applied to an actual heterogeneous reservoir, and a better effect of oil increase and water reduction is achieved.

Coalbed Methane (CBM) has become an important gas resource in recent decades. The brittle property of coal matrix and overactive operation make the migration of coal fines inevitable. Blockage by coal fines that plugs flow paths is a non-negligible issue that results in a significant decline in gas production. By setting different experimental conditions with the following factors—coal fines concentration of the mixture displacing fluids, constant flow pump rate, inlet pressure, outlet pressure and confining pressure—six experimental schemes were designed to investigate the two-phase water and coal fines flow in natural core samples. When the differential pressure and flooding flow reach a pseudo-steady status, the equivalent permeability of coal samples can be approximately calculated considering coal fines migration. Furthermore, the influences of coal fines migration on the cleat opening and permeability variation are analyzed in the porous coal medium. The study will benefit CBM development and save pump maintenance costs. In this work, we found that maintaining the differential pressure for a longer period may result in new cleat openings and severe coal rock damage during the single-phase water flooding process. While coal fines may plug some natural cleats and pores, especially in the core samples with micro-cleats during the two-phase flooding stage, coal fines migration significantly reduces the equivalent permeability and dewatering ability of the coal rock in the earlier flooding. While enlarging the differential pressure in two-phase water and fines flooding, breakthrough of coal fines from the samples contributes to widened cleats. While coal fines are difficult to flood into the core pores for low-permeability core samples, coal fines gather in the inlet, and it is also difficult to reach the pseudo-steady status even under higher differential pressure. The damage to permeability mainly occurs in the early stage of coal fines migration, and an abrupt increase in the flow velocity can damage reservoirs and induce substantial coal fines generation. Thus, maintaining a stable effective strength and a controlled depressurization rate during drainage can effectively constrain coal fines output and decrease permeability damage within coal reservoirs.

Cyclic water huff and puff (CWHP) has proven to be an attractive alternative to improve oil production performance after depletion-drive recovery in fractured-vuggy carbonate reservoirs. However, due to the impact of strong heterogeneity, multiple types of fractured-vuggy medium, poor connectivity, complex flow behaviors and oil-water relationship, CWHP is merely suitable for specific types of natural fractured-vuggy medium, usually causing a great difference in actual oil-yielding effect. It remains a great challenge for accurate evaluation of CWHP adaptability and quantitative prediction of production performance in fractured-vuggy carbonate reservoir, which severely restricts the application of CWHP. For this study, we firstly enable the newly developed fuzzy grey relational analysis to quantify the adaptability of CWHP. With production history of several targeted producers, the accuracy of the proposed method is validated. Based on the traditional percolation theory and waterflood mechanisms in various types of fractured-vuggy medium, a quantitative prediction model for cyclic water cut f_wp and increased recovery factor ΔR is presented. The CWHP production performance is discussed by using the Levenberg-Marquardt algorithm for history matching. With a better understanding of the f_wp ~ ΔR curve characteristics in different types of fractured-vuggy medium, proper strategies or measures for potential-tapping remaining oil are provided. This methodology can also offer a good basis for engineers and geologists to develop other similar reservoirs with high efficiency.