This article, written by Special Publications Editor Adam Wilson, contains highlights of paper SPE 164837, ’Multiscale Simulation of WAG Flooding in Naturally Fractured Reservoirs,’ by Mohamed Ahmed Elfeel, Adnan Al-Dhahli, Sebastian Geiger, SPE, and Marinus I.J. van Dijke, Heriot-Watt University, prepared for the 2013 EAGE Annual Conference and Exhibition/SPE Europec, London, 10-13 June. The paper has not been peer reviewed. Naturally fractured reservoirs (NFRs) contain a significant amount of remaining petroleum reserves and are now being considered for water-alternating-gas (WAG) flooding as secondary or tertiary recovery. The authors face the challenge of reservoir simulation of WAG by building models at various scales, starting with pore scale and expanding to an intermediate scale and then to reservoir scale. They show how pore-network modeling and fine-grid modeling where the fractures and matrix are represented explicitly can be used to increase the accuracy of numerical simulations at field scale. Introduction A significant portion of the world’s remaining petroleum resources is located in NFRs, including supergiant fields in the Middle East. A detailed understanding of the recovery processes involved in extracting hydrocarbons from NFRs by use of enhanced-oil-recovery techniques is key to increasing ultimate recovery for such reservoirs. Waterflooding has been used with various degrees of success in NFRs. For unfavorable (i.e., mixed- to oil-wet) matrix wettability, however, waterflooding can be ineffective. Gas/oil gravity drainage (GOGD) provides an important drive mechanism that can be effective irrespective of the rock wettability. In NFRs, fractures increase the exposure of the injected gas to oil in reservoir rock, which renders GOGD more effective than it is in unfractured reservoirs. Consequently, gas injection has been applied in many NFRs. However, because the gas mobility is very high compared with that of water and oil, so is the risk of bypassed oil and gravity override, which can lead to very early gas breakthrough. This is particularly true for NFRs. WAG flooding combines the merits of the two injection fluids on macroscopic and microscopic scales while stabilizing the injection front, delaying breakthroughs, and, therefore, leading to increased oil recovery compared with continuous water or gas injection. Reservoir simulation of WAG injection is very challenging because a representative three-phase saturation model is required to predict relative permeability and capillary pressure as water and gas saturations increase and decrease alternately. Three-phase relative permeability and capillary pressure data are extremely difficult to measure experimentally. Even if experimental evaluation becomes feasible, an infinite number of saturation paths can occur in the reservoir. This necessitates the use of empirical, or interpolation, models to predict three-phase relative permeability and capillary pressure from two-phase experiments. In NFRs, capillary pressure and relative permeability functions have a major effect on fluid exchange between matrix blocks and fractures. Predicting the effects of the interplay of viscous, capillary, and gravity forces is challenging because fluid flow is viscous-dominated in the fractures while transfer between fractures and matrix blocks is dominated by capillary and gravity forces. Because most of the oil is contained inside the matrix, capillary and gravity forces can be more important in NFRs than in conventional reservoirs.