Flow through partially frozen pores in granular media containing ice or gas hydrate plays an essential role in diverse phenomena including methane migration and frost heave. As freezing progresses, the frozen phase grows in the pore space and constricts flow paths so that the permeability decreases. Previous works have measured the relationship between permeability and volumetric fraction of the frozen phase, and various correlations have been proposed to predict permeability change in hydrology and the oil industry. However, predictions from different formulae can differ by orders of magnitude, causing great uncertainty in modeling results. We present a floating random walk method to approximate the porous flow field and estimate the effective permeability in isotropic granular media with specified particle size distributions, without solving for the entire flow field in the pore space. In packed spherical particles, the method compares favorably with the Kozeny‐Carman formula. We further extend this method with a probabilistic interpretation of the volumetric fraction of the frozen phase, simulate the effect of freezing in irregular pores, and predict the evolution of permeability. Employing no adjustable parameters, our results can provide insight into the coupling between phase transitions and permeability change, which plays important roles in hydrate formation and dissociation, as well as in the thawing and freezing of permafrost and ice‐bed coupling beneath glaciers.
- Award ID(s):
- 1735862
- PAR ID:
- 10366289
- Date Published:
- Journal Name:
- Journal of Glaciology
- ISSN:
- 0022-1430
- Page Range / eLocation ID:
- 1 to 16
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract One of the major problems during gas injection in unconventional reservoirs is asphaltene precipitation and deposition. Asphaltenes can reduce the pore throat in the reservoir and plug the surface and subsurface equipment during the production process, thus, result in oil production reduction with significant financial consequences. The impact of carbon dioxide (CO2) gas injection on asphaltene deposition in unconventional reservoirs still poorly investigated. This research investigates the impact of CO2 gas injection on asphaltene aggregation in ultra-low-permeability pore structures, mainly present in unconventional shale resources. First, the minimum miscibility pressure (MMP) of crude oil with CO2 was determined using the slim tube technique. Then, several CO2 injection pressures were selected to conduct the filtration experiments using a specially designed filtration apparatus. All pressures selected were below the MMP. Various sizes of filter paper membranes were used to study the effect of pore structure on asphaltene deposition. The results showed that asphaltene weight percent was increased by increasing the pressure and a significant asphaltene weight percentage was observed on smaller pore size structures of the filter membranes. The visualization tests revealed the process of asphaltene precipitation and deposition and showed that asphaltene particles and clusters were precipitated after one hour and fully deposited in the bottom of the test tube after 12 hours. High-resolution photos of filter paper membranes were presented using microscopy imaging and scanning electron microscopy (SEM) analysis; these photos highlighted the asphaltene particles inside the filter paper membranes and pore plugging was observed. The study's findings will contribute to a better understanding of the main factors influencing the stability of asphaltene particles in crude oil under immiscible CO2 injection pressure, particularly in nano pores, which are predominant in shale unconventional resources.
-
more » « less
Abstract The electrokinetic transport mechanisms of multispecies ions through 3‐D nanoporous rocks with chemical reaction at the solid‐aqueous solution interfaces are investigated. We systematically study the multiphysics transport phenomena by considering either inhomogeneous (local surface charge based on local pH and ion concentrations) or prescribed homogeneous surface charge at solid‐aqueous solution interface while the pores are screened via electric double layers. We develop a lattice Boltzmann numerical framework to solve the set of governing equations (Poisson‐Nernst‐Planck plus Navier–Stokes). Our modeling results reveal that the averaged local electric potential of the nanoporous rock is significantly underestimated (about 83%) when a homogeneous surface charge is prescribed based on the bulk solution properties. It is shown that increasing the porosity of the nanoporous media considerably increases the absolute values and inhomogeneity of the surface charge, which means that while the electric double layers screened the pores, increasing the porosity enhances the ion selectivity of the porous medium. When the scenario with inhomogeneous charge is taken into account, the predicted electroosmotic permeability and tortuosity are higher in comparison with the prescribed homogeneous case. Moreover, we have studied the electrostatic tortuosity, coupling coefficient, and the effective excess charge density of the nanoporous rocks. The results demonstrate that ignoring the inhomogeneity of surface charges may cause erroneous prediction of the ion transport through porous rocks with chemically active surfaces.
-
Directional graphene aerogels (DGAs) are proposed as electrode materials to alleviate ionic and mass transport issues in organic redox flow batteries (ORFBs). DGAs with high pore directionality would provide low resistance channels for effective ionic charge and liquid electrolyte transport in these devices. DGAs’ porous and directional characteristics can be controlled by the growth of ice crystals during freeze casting, which is influenced by the self-diffusivity of water, phase change driving forces, water−ice graphene interactions, and convection in the water−graphene media. It is found that mass transport-related properties of DGAs, including pore size and directionality, show a significant dependence on freezing temperature, graphene oxide (GO) loadings, and synthesis vessel diameter-to-height ratio (D/H). For the freezing temperature change from −20 to −115 °C, the average pore size progressively decreased from 120 to 20 μm, and the pore directionality transitioned from lamellar to ill-defined structures. When GO loadings were increased from 2 to 10 mg/mL at a fixed freezing temperature, pore size reduction was observed with less defined directionality. Furthermore, the pore directionality diminished with an increased width-to-height aspect ratio of DGA samples due to the buoyancy-driven convective circulation, which interfered with the directional ice/pore growth. Understanding the comprehensive effects of these mechanisms enables the controlled growth of ice crystals, leading to graphene aerogels with highly directional microstructures.more » « less
-
null (Ed.)In this paper, we consider an important problem for modeling complex coupled phenomena in porous media at multiple scales. In particular, we consider flow and transport in the void space between the pores when the pore space is altered by new solid obstructions formed by microbial growth or reactive transport, and we are mostly interested in pore-coating and pore-filling type obstructions, observed in applications to biofilm in porous media and hydrate crystal formation, respectively. We consider the impact of these obstructions on the macroscopic properties of the porous medium, such as porosity, permeability and tortuosity, for which we build an experimental probability distribution with reduced models, which involves three steps: (1) generation of independent realizations of obstructions, followed by, (2) flow and transport simulations at pore-scale, and (3) upscaling. For the first step, we consider three approaches: (1A) direct numerical simulations (DNS) of the PDE model of the actual physical process called BN which forms the obstructions, and two non-DNS methods, which we call (1B) CLPS and (1C) LP. LP is a lattice Ising-type model, and CLPS is a constrained version of an Allen–Cahn model for phase separation with a localization term. Both LP and CLPS are model approximations of BN, and they seek local minima of some nonconvex energy functional, which provide plausible realizations of the obstructed geometry and are tuned heuristically to deliver either pore-coating or pore-filling obstructions. Our methods work with rock-void geometries obtained by imaging, but bypass the need for imaging in real-time, are fairly inexpensive, and can be tailored to other applications. The reduced models LP and CLPS are less computationally expensive than DNS, and can be tuned to the desired fidelity of the probability distributions of upscaled quantities.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jog.2022.1
