Sharif Digital Repository

SUMMARY: In this paper, a numerical model is developed for the fully coupled hydro-mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non-wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two-phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid... 

A genetic-based pore network extraction method from micro-computed tomography (micro-CT) images is proposed in this paper. Several variables such as the number, radius and location of pores, the coordination number, as well as the radius and length of the throats are used herein as the optimization parameters. Two approaches to generate the pore network structure are presented. Unlike previous algorithms, the presented approaches are directly based on minimizing the error between the extracted network and the real porous medium. This leads to the generation of more accurate results while reducing required computational memories. Two different objective functions are used in building the... 

It is well known that the displacement efficiency of EOR processes is mainly affected by wettability of porous medium; however, the role of nanoparticles on wettability alteration of pores surfaces remains a topic of debate in the literature. Furthermore, a little is known about how the dispersed silica nanoparticles affect the microscopic/macroscopic recovery efficiency of heavy oils during common immiscible EOR processes such as water flooding. In this study, a series of injection experiments was performed on five-spot glass micromodel which is initially saturated with the heavy oil. Distilled water and dispersed silica nanoparticles in water (DSNW) at different values of weight percent... 

Spontaneous water imbibition into the matrix blocks is known as the main mechanism for increased oil recovery from naturally fractured oil reservoirs. The rate of oil recovery and its ultimate value is mostly affected by wettability of the rocks and their pore structure. Oil viscosity also greatly influences the rate of oil recovery. A novel experimental model was utilized to study the imbibition mechanism under different wettability conditions. Matrix blocks made from different grain types and size distributions of glass beads were saturated with two different types of synthetic oil, to mimic the oil-saturated matrixes. The wetting characteristic of the models used in this study were... 

There are many fundamental problems with the injection of nano-zero-valent iron (NZVI) particles to create permeable reactive barrier (PRB) treatment zone. Among them the loss of medium porosity or pore blocking over time can be considered which leads to reduction of permeability and bypass of the flow and contaminant plume up-gradient of the PRB. Present study provides a solution for such problems by confining the target zone for injection to the gate in a funnel-and-gate configuration. A laboratory-scale experimental setup is used in this work. In the designed PRB gate, no additional material from porous media exists. NZVI (d 50 = 52 ± 5 nm) particles are synthesized in water mixed with... 

In this paper, a mathematical model is presented for the analysis of dynamic fracture propagation in the saturated porous media. The solid behavior incorporates a discrete cohesive fracture model, coupled with the flow in porous media through the fracture network. The double-nodded zero-thickness cohesive interface element is employed for the mixed mode fracture behavior in tension and contact behavior in compression. The crack is automatically detected and propagated perpendicular to the maximum effective stress. The spatial discretization is continuously updated during the crack propagation. Numerical examples from the hydraulic fracturing test and the concrete gravity dam show the... 

It is well known that the oil recovery is affected by wettability of porous medium; however, the role of nanoparticles on wettability alteration of medium surfaces has remained a topic of debate in the literature. Furthermore, there is a little information of the way dispersed silica nanoparticles affect the oil recovery efficiency during polymer flooding, especially, when heavy oil is used. In this study, a series of injection experiments were performed in a five-spot glass micromodel after saturation with the heavy oil. Polyacrylamide solution and dispersed silica nanoparticles in polyacrylamide (DSNP) solution were used as injected fluids. The oil recovery as well as fluid distribution in... 

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic... 

A comprehensive study was performed to analyze the unsteady laminar flow characteristics around a porously covered, a fully porous, and a solid squared section cylinder located in the middle of a plane channel. In order to simulate fluid flow inside porous media and porous–fluid interface accurately (minimizing modeling error), the porous region was analyzed in pore scale, using LBM. Additionally, to minimize the LBM-related compressibility error through the porous region, a multi-block multiple relaxation time lattice Boltzmann method (MRT-LBM) was used. Also, to decrease CPU time, a Navier–Stokes flow solver, based on finite volume method and SIMPLE algorithm, was coupled with MRT-LBM to... 

The majority of literature regarding optimized Latin hypercube sampling (OLHS) is devoted to increasing the efficiency of these sampling strategies through the development of new algorithms based on the combination of innovative space-filling criteria and specialized optimization schemes. However, little attention has been given to the impact of the initial design that is fed into the optimization algorithm, on the efficiency of OLHS strategies. Previous studies, as well as codes developed for OLHS, have relied on one of the following two approaches for the selection of the initial design in OLHS: (1) the use of random points in the hypercube intervals (random LHS), and (2) the use of... 

The sloshing phenomena in a partially filled tank can affect its stability. Modifications of tank instability due to the movement of the tank carrier, are key design points for the stability of a carrier. Even though the sloshing phenomenon has already been investigated using the BEM-FDM technique, the research in this paper covers this phenomenon in a porous media, which is new in 2-D coordinates. For this purpose, a Laplace equation has been used for potential flow, and kinematic and dynamic boundary conditions have been applied to the free surface. Also, a formulation has been developed for a free surface in porous media. BEM has been used for solving the governing equation and FDM... 

A Green Element numerical formulation is used to solve the time-dependent nonlinear one-dimensional counter-current spontaneous imbibition diffusion equation in which water enters a water wet rock spontaneously while oil escapes by flowing in the opposite direction. The Green Element Method (GEM) is an element by element approach of the boundary element method. In this new method, by generating large sparse global matrices and yet taking advantage of properties of Green's function, solution of more complicated physical problem is achievable while at the same time much less computational effort is needed rather than boundary element method (BEM). By discretizing both the boundary and problem... 

Including gravity and wettability effects, a full analytical solution for the frontal flow period for 1D counter-current spontaneous imbibition of a wetting phase into a porous medium saturated initially with non-wetting phase at initial wetting phase saturation is presented. The analytical solution applicable for liquid-liquid and liquid-gas systems is essentially valid for the cases when the gravity forces are relatively large and before the wetting phase front hits the no-flow boundary in the capillary-dominated regime. The new analytical solution free of any arbitrary parameters can also be utilized for predicting non-wetting phase recovery by spontaneous imbibition. In addition, a new... 

In this work, hybrid of silica nanoparticles (NPs) with sulfonated copolymer has been prepared in order to improve the copolymer properties for Enhanced Oil Recovery (EOR). Some tests are done to find the effectiveness of injecting fluid on Recovery Factor (RF), such as Contact Angle (CA), interfacial tension (IFT), inherent viscosity, and eventually Micro-Model (MM) flooding. In CA test, wettability alteration from Oil-Wet (OW) to Water-Wet (WW) is reached by sedimentation and adsorption of NPs on the rock slice. In addition, IFT reduction is obtained by increasing the NaCl concentration. The viscosity change is investigated for the hybrid and copolymer under simulated high-temperature oil... 

Diffusivity equation which can provide us with the pressure distribution, is a Partial Differential Equation (PDE) describing fluid flow in porous media. The quadratic pressure gradient term in the diffusivity equation is nearly neglected in hydrology and petroleum engineering problems such as well test analysis. When a compressible liquid is injected into a well at high pressure gradient or when the reservoir possess a small permeability value, the effect of ignoring this term increases. In such cases, neglecting this parameter can result in high errors. Previous models basically focused on numerical and semi-analytical methods for semi-infinite domain. To the best of our knowledge, no... 

A sequential implicit numerical method based on discrete-fracture model and the Galerkin Finite Element method, for time-dependent coupled fluid flow and geomechanics problems in fractured subsurface formations is presented. Discrete-fracture model has been used to explicitly represent the fracture network inside porous media. The Galerkin Finite Element method with adaptive unstructured gridding is implemented to numerically solve the spatially three-dimensional and time-dependent problem. The presented method is validated with previously obtained solutions. Two problems are numerically solved by applying the presented methodology in a three-dimensional fractured petroleum reservoir under... 

Accurate determination of relative permeability (kr) curves and their hysteresis is needed for reliable prediction of the performance of oil and gas reservoirs. A few options (e.g., Carlson, Killough and Jargon models) are available in commercial reservoir simulators to account for hysteresis in kr curves under two-phase flow. Two-phase kr curves are also needed for estimating kr hysteresis under three-phase flow during water-alternating-gas (WAG) injection. Although, most oil reservoirs are mixed-wet, the existing hysteresis predictive approaches have been developed based on water-wet conditions. Experimentally measured data are needed to assess the performance of these methodologies under... 

Generally, in most areas, groundwater level is deep and heat sources (e.g., energy piles)are embedded in unsaturated soil media. Therefore, in order to accurately analyze the soil response close to heat sources, both heat and moisture transport in unsaturated soil domain should be considered. Thermal loading changes the moisture content in the porous media. In this study, the energy conservation and mass fluid continuity equations derived from hydrothermal analysis of a partially saturated soil medium are considered in cylindrical coordinate system. To make the analytical solution possible, partial differential equations (PDEs)are turned into ordinary differential equations (ODEs), through... 

Practical designs of non-aqueous phase liquids (NAPLs) remediation strategies require reliable modeling of interphase mass transfer to predict the retraction of NAPL during processes such as dissolution. In this work, the dissolution process of NAPL during two-phase flow in heterogeneous porous media is studied using pore-network modeling and micromodel experiments. A new physical-experimental approach is proposed to enhance the prediction of the dissolution process during modeling of interphase mass transfer. In this regard, the normalized average resident solute concentration is evaluated for describing the dissolution process at pore-level. To incorporate the effect of medium... 

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase...