Graduation date: 1997

Geological carbon sequestration, as a method of atmospheric greenhouse gas reduction, is at the technological forefront of the climate change movement. During sequestration, carbon dioxide (CO���) gas effluent is captured from coal fired power plants and is injected into a storage saline aquifer or depleted oil reservoir. In an effort to fully understand and optimize CO��� trapping efficiency, the capillary trapping mechanisms that immobilize subsurface CO��� were analyzed at the pore-scale. Pairs of proxy fluids representing the range of in situ supercritical CO��� and brine conditions were used during experimentation. The two fluids (identified as wetting and non-wetting) were imbibed and drained from a flow cell apparatus containing a sintered glass bead column. Experimental and fluid parameters, such as interfacial tension, fluid viscosities and flow rate, were altered to characterize their relative impact on capillary trapping. Computed x-ray microtomography (CMT) was used to identify immobilized CO��� (non-wetting fluid) volumes after imbibition and drainage events. CMT analyzed data suggests that capillary behavior in glass bead systems do not follow the same trends as in consolidated natural material systems. An analysis of the disconnected phases in both the initial and final flood events indicate that the final (residual) amount of trapped non-wetting phase has a strong linear dependence on the original amount of non-wetting phase (after primary imbibition), which corresponds to the amount of gas or oil present in the formation prior to CO��� injection. More importantly, the residual trapped gas was also observed to increase with increasing non-wetting fluid phase viscosity. This suggests that CO��� sequestration can be optimized in two ways: through characterization of the trapped fluid present in the formation prior to injection and through alterations to the viscosity of supercritical CO2. / Graduation date: 2013

Nonlinearity in the Navier-Stokes equations can originate from a variety of sources, such as contributions stemming from the advective term, constitutive closure models or external factors such as chemical reactions and capillarity. Needless to say, a combination of any of the above sources has the potential to exasperate the problem significantly. This dissertation explores cases that predominantly feature advective and/or capillary effects. In particular, we first consider the inertia-dominated problem of single-phase flow past a confined square cylinder, followed by a study focused on the low-Re dynamics of rigid particles straddling non-planar interfaces. The first part of the thesis investigates transient, three-dimensional, incompressible and isothermal flow of a Newtonian fluid past a symmetrically confined obstacle at zero incidence. Results from both Laser Doppler Velocimetry (LDV) experiments and direct simulations upto Re = 250 have been reported. Beyond the onset of instability (Recr ≈ 58), an inflexion point around Re ≈ 115 is detected for the Strouhal number with no evidence of hysteresis in any of the measurements. Furthermore, incommensurate frequencies observed in the range 127 ≤ Re ≤ 175 suggest a quasi-periodic transition to three-dimensionality. This is shown to be followed by an intermediate periodic window starting around Re ≈ 180. Fourier analysis and spanwise velocity correlations are then used to characterize the observed phenomena. Subsequent analysis of consolidated data suggest that only a parametric variation of transverse and spanwise blockage ratios can bring closure to the subject of bluff-body wake transitions. The second part of the thesis implements and validates a physically consistent continuum model for the Moving Contact Line (MCL) through direct simulations. After elaborately discussing the MCL conundrum, a fundamental framework for the simulations is outlined in a theoretical orientation which combines the Level set method with a Fictitious domain approach in a finite-element scheme. The thesis objectives are then realized through simulation of various case studies that show favorable comparisons with theoretical and/or published experimental data. In short, the current work successfully illustrates the potential of novel boundary conditions (such as the GNBC) to accurately describe MCL dynamics. / Chemical Engineering

Thin liquid films driven by surface tension gradients are studied in diverse applications, including the spreading of a droplet and fluid flow in the lung. The nonlinear partial differential equations that govern thin films are difficult to solve analytically, and must be approached through numerical simulations. We describe the development of a numerical solver designed to solve a variety of thin film problems in two dimensions. Validation of the solver includes grid refinement studies and comparison to previous results for thin film problems. In addition, we apply the solver to a model of surfactant spreading and make comparisons with theoretical and experimental results.

It is widely known in oil industry that changes in fluid flow conditions such as water breakthrough or unsteady flow due to well shut-in can lead to sand destabilization, with a possible consequent sand production. In this research, different flow situations are incorporated into stress and stability analysis for the region around a wellbore producing oil from weak or unconsolidated sands, and the analyses involve strength weakening, stress redistribution, and decrease of rock stiffness. Two main mechanisms, chemical reactions of rock with formation water and variations of rock capillary strength, are identified and analyzed to study strength weakening after water breakthrough, both qualitatively and quantitatively. Using theories from particle mechanics, rock mechanics, and interfacial science, four novel capillarity models are developed and verified to analytically capture the physical behaviors of capillary strength at the grain scale. Based on model calculations, significantly better understanding of strength behavior in two-phase fluid environments is achieved. Based on a simplified model that can conservatively but efficiently quantify capillary strength with only two input parameters (i. e. particle radius and water saturation), a verified new method that physically calculates pore pressure in a multiphase environment, and a coupled poro-inelastic stress model, the redistributions of effective stresses with water saturation around a wellbore are solved. In terms of stress changes and growth of a plastic radius defining shear-failure zone, the effects of different stability factors, including capillarity through water-oil menisci, pore pressure changes due to the variations of fluid relative permeabilities, and loss of strength through chemical reactions of water-sensitive cementation materials, are quantified and compared in order to clarify when and how they contribute to sand production after water breakthrough. The nonlinearities of rock elastic properties in stressed and biphasic fluid environments is analytically addressed, based on an improved nonlinear theory that considers both a failure-based mechanism and a confining-stress-based mechanism, the strength model, and the coupled stress model. The calculations demonstrate the redistributions of stress-dependent rock stiffness around a wellbore and its evolution with increase of water saturation, clarify the relative importance of each mechanism in reducing rock stiffness, and fundamentally explain why current predictive technologies are invalid when water appears in a flowing wellbore. To quantify the effect of well shut-down on rock stability, the redistributions of fluid pressure in reservoir are analytically solved and coupled with the stress model, while the water hammer equations provide a boundary condition for the bottom-hole pressure. This approach allows direct solution of the relationships among fluid properties, rock properties and production parameters, within the context of rock stability. The proposed new approaches and models can be applied to evaluate sand production risk in multiphase and unsteady fluid flow environment. They can also serve as points of departure to develop more sophisticated models, or to develop more useful constitutive laws for numerical solutions.

It is widely known in oil industry that changes in fluid flow conditions such as water breakthrough or unsteady flow due to well shut-in can lead to sand destabilization, with a possible consequent sand production. In this research, different flow situations are incorporated into stress and stability analysis for the region around a wellbore producing oil from weak or unconsolidated sands, and the analyses involve strength weakening, stress redistribution, and decrease of rock stiffness. Two main mechanisms, chemical reactions of rock with formation water and variations of rock capillary strength, are identified and analyzed to study strength weakening after water breakthrough, both qualitatively and quantitatively. Using theories from particle mechanics, rock mechanics, and interfacial science, four novel capillarity models are developed and verified to analytically capture the physical behaviors of capillary strength at the grain scale. Based on model calculations, significantly better understanding of strength behavior in two-phase fluid environments is achieved. Based on a simplified model that can conservatively but efficiently quantify capillary strength with only two input parameters (i. e. particle radius and water saturation), a verified new method that physically calculates pore pressure in a multiphase environment, and a coupled poro-inelastic stress model, the redistributions of effective stresses with water saturation around a wellbore are solved. In terms of stress changes and growth of a plastic radius defining shear-failure zone, the effects of different stability factors, including capillarity through water-oil menisci, pore pressure changes due to the variations of fluid relative permeabilities, and loss of strength through chemical reactions of water-sensitive cementation materials, are quantified and compared in order to clarify when and how they contribute to sand production after water breakthrough. The nonlinearities of rock elastic properties in stressed and biphasic fluid environments is analytically addressed, based on an improved nonlinear theory that considers both a failure-based mechanism and a confining-stress-based mechanism, the strength model, and the coupled stress model. The calculations demonstrate the redistributions of stress-dependent rock stiffness around a wellbore and its evolution with increase of water saturation, clarify the relative importance of each mechanism in reducing rock stiffness, and fundamentally explain why current predictive technologies are invalid when water appears in a flowing wellbore. To quantify the effect of well shut-down on rock stability, the redistributions of fluid pressure in reservoir are analytically solved and coupled with the stress model, while the water hammer equations provide a boundary condition for the bottom-hole pressure. This approach allows direct solution of the relationships among fluid properties, rock properties and production parameters, within the context of rock stability. The proposed new approaches and models can be applied to evaluate sand production risk in multiphase and unsteady fluid flow environment. They can also serve as points of departure to develop more sophisticated models, or to develop more useful constitutive laws for numerical solutions.

No description available.

No description available.

Thesis (M.S.)--Oregon State University, 2010. / Printout. Includes bibliographical references (leaves 66-69). Also available on the World Wide Web.

The Steam-Assisted Gravity Drainage (SAGD) Process involves the injection of steam through a horizontal well and the production of heavy oil through a lower horizontal well. Several authors have tried to model this process using analytical, semi-analytical and fully numerical means. In this dissertation, we improve the predictive ability of previous models by accounting for the effect of anisotropy, the effect of heat transfer on capillarity and the effect of water-in-oil (W/O) emulsion formation and transport which serves to enhance heat transfer during SAGD. We account for the effect of anisotropy during SAGD by performing elliptical transformation of the resultant gravity head and resultant oil drainage vectors on to a space described by the vertical and horizontal permeabilities. Our results, show that unlike for the isotropic case, the effect of anisotropy is time dependent and there exists a given time beyond which it ceases to have any effect on SAGD rates. This result will impact well spacing design and optimization during SAGD. Butler et al. (1981) derived their classical SAGD model by solving a 1-D heat conservation equation for single phase flow. This model has excellent predictive capability at experimental scales but performs poorly at field scales. By assuming a linear saturation -- temperature relationship, Sharma and Gates (2010b) developed a model that accounts for multiphase flow ahead of the steam chamber interface. In this work, by decomposing capillary pressure into its saturation and temperature components, we coupled the mass and energy conservation equations and showed that the multi-scale, multiphase flow phenomenon occurring during SAGD is the classical Marangoni (or thermo-capillary) effect which can be characterized by the Marangoni number. At low Marangoni numbers (typical of experimental scales) we get the Butler solution while at high Marangoni numbers (typical of field scales), we approximate the Sharma and Gates solution. The Marangoni flow concept was extended to the Expanding Solvent SAGD (ES-SAGD) process and our results show that there exists a given Marangoni number threshold below which the ES-SAGD process will not fare better than the SAGD process. Experimental results presented in Sasaki et al. (2002) demonstrate the existence of water-in-oil emulsions adjacent to the steam chamber wall during SAGD. In this work we show that these emulsions enhanced heat transfer at the chamber wall and hence oil recovery. We postulate that these W/O emulsions are principally hot water droplets that carry convective heat energy. We perform calculations to show that their presence can practically double the effective heat transfer coefficient across the steam chamber interface which overcomes the effect of reduced oil rates due to the increased emulsified phase viscosity. Our results also compared well with published experimental data. The SAGD (and ES-SAGD) process is a short length-scaled process and hence, short length-scaled phenomena (typically ignored in other EOR or conventional processes) such as thermo-capillarity and in-situ emulsification should not be ignored in predicting SAGD recoveries. This work will find unique application in predictive models used as fast proxies for predicting SAGD recovery and for history matching purposes. / text