The last decade has seen a strong increase of research into flows in fractured porous media, mainly related to subsurface processes, but also in materials science and biological applications. Connected fractures totally dominate flow-patterns, and their representation is therefore a critical part in model design. Due to the fracture's characteristics as approximately planar discontinuities with an extreme size to width ratio, they challenge standard macroscale mathematical and numerical modeling of flow based on averaging. Thus, over the last decades, various, and also fundamentally different, approaches have been developed. This paper reviews common conceptual models and discretization approaches for flow in fractured porous media, with an emphasis on the dominating effects the fractures have on flow processes. In this context, the paper discuss the tight connection between physical and mathematical modeling and simulation approaches. Extensions and research challenges related to transport, multi-phase flow and fluid-solid interaction are also commented on. arXiv:1805.05701v1 [physics.geo-ph]

Conventional vertically integrated models for CO 2 storage usually adopt a vertical equilibrium (VE) assumption, which states that due to strong buoyancy, CO 2 and brine segregate quickly, so that the fluids can be assumed to have essentially hydrostatic pressure distributions in the vertical direction. However, the VE assumption is inappropriate when the time scale of fluid segregation is not small relative to the simulation time. By casting the vertically integrated equations into a multiscale framework, a new vertically integrated model can be developed that relaxes the VE assumption, thereby allowing vertical dynamics to be modeled explicitly. The model maintains much of the computational efficiency of vertical integration while allowing a much wider range of problems to be modeled. Numerical tests of the new model, using injection scenarios with typical parameter sets, show excellent behavior of the new approach for homogeneous geologic formations.

Recent studies suggest the possibility of CO2 sequestration in depleted shale gas formations, motivated by large storage capacity estimates in these formations. Questions remain regarding the dynamic response and practicality of injection of large amounts of CO2 into shale gas wells. A two-component (CO2 and CH4) model of gas flow in a shale gas formation including adsorption effects provides the basis to investigate the dynamics of CO2 injection. History-matching of gas production data allows for formation parameter estimation. Application to three shale gas-producing regions shows that CO2 can only be injected at low rates into individual wells and that individual well capacity is relatively small, despite significant capacity variation between shale plays. The estimated total capacity of an average Marcellus Shale well in Pennsylvania is 0.5 million metric tonnes (Mt) of CO2, compared with 0.15 Mt in an average Barnett Shale well. Applying the individual well estimates to the total number of existing and permitted planned wells (as of March, 2015) in each play yields a current estimated capacity of 7200-9600 Mt in the Marcellus Shale in Pennsylvania and 2100-3100 Mt in the Barnett Shale.

Naturally Fractured Reservoirs (NFR's) have received little attention as potential CO2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO2 tends to be nonwetting in target formations and capillary forces will keep CO2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO2 plume. Numerical simulations are a powerful tool to understand the physics behind brine‐CO2 flow in NFR's. Dual‐porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual‐porosity models are based on crude approximations of the matrix‐fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO2 transfer to the matrix. This work presents an assessment of CO2 storage potential in NFR's using dual‐porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO2 in the matrix blocks and the maximum CO2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual‐porosity simulations. We then analyze how injection rates should be limited in order to avoid early spill of CO2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO2.

A recent macroscopic theory of biphasic flow in porous media ͓R. Hilfer, Phys. Rev. E 73, 016307 ͑2006͔͒ has proposed to treat microscopically percolating fluid regions differently from microscopically nonpercolating regions. Even in one dimension the theory reduces to an analytically intractable set of ten coupled nonlinear partial differential equations. This paper reports numerical solutions for three different initial and boundary value problems that simulate realistic laboratory experiments. All three simulations concern a closed column containing a homogeneous porous medium filled with two immiscible fluids of different densities. In the first simulation the column is raised from a horizontal to a vertical orientation inducing a buoyancy-driven fluid flow that separates the two fluids. In the second simulation the column is first raised from a horizontal to a vertical orientation and subsequently rotated twice by 180°to compare the resulting stationary saturation profiles. In the third simulation the column is first raised from horizontal to vertical orientation and then returned to its original horizontal orientation. In all three simulations imbibition and drainage processes occur simultaneously inside the column. This distinguishes the results reported here from conventional simulations based on existing theories of biphasic flows. Existing theories are unable to predict flow processes where imbibition and drainage occur simultaneously. The approximate numerical results presented here show the same process dependence and hysteresis as one would expect from an experiment.

[1] Several models for two-phase flow in porous media identify trapping and connectivity of fluids as an important contribution to macroscale hysteresis. This is especially true for hysteresis in relative permeabilities. The trapping models propose trajectories from the initial saturation to the end saturation in various ways and are often based on experiments or pore-network model results for the endpoints. However, experimental data or pore-scale model results are often not available for the trajectories, that is, the fate of the connectivity of the fluids while saturation changes. Here, using a quasi static pore-network model, supported by a set of pore-scale laboratory experiments, we study how the topology of the fluids changes during drainage and imbibition including first, main and scanning curves. We find a strong hysteretic behavior in the relationship between disconnected nonwetting fluid saturation and the wetting fluid saturation in a water-wet medium. The coalescence of the invading nonwetting phase with the existing disconnected nonwetting phase depends critically on the presence (or lack thereof) of connected nonwetting phase at the beginning of the drainage process as well as on the pore geometry. This dependence involves a mechanism we refer to as ''reversible corner filling.'' This mechanism can also be seen in laboratory experiments in volcanic tuff. The impact of these pore-network model results on existing macroscopic models is discussed.

The concepts of relative permeability and capillary pressure are crucial for the accepted traditional theory of two phase flow in porous media. Recently, a theoretical approach was introduced that does not require these concepts as input (Hilfer, Physica A, 359:119-128, 2006a; Phys. Rev. E, 73:016307, 2006b). Instead it was based on the concept of hydraulic percolation of fluid phases. This paper presents the first numerical solutions of the coupled nonlinear partial differential equations introduced in Hilfer (Phys. Rev. E, 73:016307, 2006b). Approximate numerical results for saturation profiles in one spatial dimension have been calculated. Long time limits of dynamic time-dependent profiles are compared to stationary solutions of the traditional theory. The initial and boundary conditions are chosen to model the displacement processes that occur when a closed porous column containing two immiscible fluids of different density is raised from a horizontal to a vertical position in a gravitational field. The nature of the displacement process may change locally in space and time between drainage and imbibition. The theory gives local saturations for nonpercolating trapped fluids near the endpoint of the displacement.

Spontaneous countercurrent imbibition into a finite porous medium is an important physical mechanism for many applications, included but not limited to irrigation, CO 2 storage, and oil recovery. Symmetry considerations that are often valid in fractured porous media allow us to study the process in a onedimensional domain. In 1-D, for incompressible fluids and homogeneous rocks, the onset of imbibition can be captured by self-similar solutions and the imbibed volume scales with ffiffi t p . At later times, the imbibition rate decreases and the finite size of the medium has to be taken into account. This requires numerical solutions. Here we present a new approach to approximate the whole imbibition process semianalytically. The onset is captured by a semianalytical solution. We also provide an a priori estimate of the time until which the imbibed volume scales with ffiffi t p . This time is significantly longer than the time it takes until the imbibition front reaches the model boundary. The remainder of the imbibition process is obtained from a selfsimilarity solution. We test our approach against numerical solutions that employ parametrizations relevant for oil recovery and CO 2 sequestration. We show that this concept improves common first-order approaches that heavily underestimate early-time behavior and note that it can be readily included into dual-porosity models. Key Points: New physically based model for spontaneous imbibition Model captures transition from earlytime to late-time imbibition Model validated for different applications Correspondence to: R. March, rafael.march@pet.hw.ac.uk Citation: March, R., F. Doster, and S. Geiger (2016), Accurate early-time and latetime modeling of countercurrent spontaneous imbibition, Water Resour.