To avoid or mitigate the unwanted water and gas content, inflow control devices (ICDs) are designed and installed in the well to disturb the water and gas breakthrough which are trying to overtake the oil inflow, water and gas coning and sand production. Smart wells with permanent downhole valves such as ICDs are used to balance production and injection in wells. A paramount issue regarding using downhole control devices is determining the required cross-sectional area of them for control of the imposed pressure drop across the device to stabilize the fluid flow. Current methods for calculating the opening size of the ICDs are mainly based on sensitivity analysis of the ICD flow area or optimization algorithms coupled with simulation models. Although these approaches are quite effective in oil field cases, they tend to be time-consuming and require demanding system models. This paper presents a fast analytical method to determine the ICD flow area validated by a genetic algorithm (GA). Analytically, a closed-form expression is introduced by manipulating Darcy’s law applicable to multi-layer injection wells with different layer properties to balance the injection profile in the reservoir pay zone, based on equalizing injected front velocity in layers with different permeability. Considering various scenarios of analytical technique, GA optimization, and sensitivity analysis scenarios for ICD cross-sectional area determination, results for oil recovery, water production, water breakthrough time, and net present value (NPV) are discussed and compared. NPV values obtained by both analytical and GA approaches are virtually identical and greater than those of other scenarios. Compared to the base field case, the analytical method improved the oil recovery by almost 1%, reduced water production by almost 91%, and synchronized the water breakthrough time of high- and low-permeability layers (from a ratio of 1.76–1.06). The proposed analytical solution proved to be capable of providing desirable results with only one reservoir simulation run in contrast to GA and sensitivity analysis scenarios which require iterative simulation runs. The proposed analytical solution outperformed the GA as it is less computationally demanding in addition to its success in case of lowering water production for the field data. The findings of this study can help for a better understanding of the situation where water injection into the oil reservoir is problematic as the layers present different permeabilities which can induce problems such as early water breakthrough from the more permeable layer and hinder the success of the water injection process. Using ICDs and a faster and more accurate approach to calculate its cross-sectional area such as the analytical method that was used in this study can greatly increase the success rate of water injection in case of oil recovery and lower the amount of the produced water.
