The study of phase behavior is important for the oil and gas industry. Many approaches have been proposed and developed for phase behavior calculation. In this thesis, an alternative method is introduced to study the phase behavior by means of minimization of Helmholtz free energy. For a system at fixed volume, constant temperature and constant number of moles, the Helmholtz free energy reaches minimum at the equilibrium state. Based on this theory, a stochastic method called Particle Swarm Optimization (PSO) algorithm, is implemented to compute the phase diagrams for several pure component and mixture systems. After comparing with experimental and the classical PT-ash calculation, we found the phase diagrams obtained by minimization of the Helmholtz Free Energy approach match the experimental and theoretical diagrams very well. Capillary effect is also considered in this thesis because it has a significant influence on the phase behavior of reservoir fluids. In this part, we focus on computing the phase envelopes, which consists of bubble and dew point lines. Both fixed and calculated capillary pressure from the Young-Laplace equation cases are introduced to study their effects on phase envelopes. We found that the existence of capillary pressure will change the phase envelopes. Positive capillary pressure reduces the dew point and bubble point temperatures under the same pressure condition, while the negative capillary pressure increases the dew point and bubble point temperatures. In addition, the change of contact angle and pore radius will affect the phase envelope. The effect of the pore radius on the phase envelope is insignificant when the radius is very large. These results may become reference for future research and study. Keywords: Phase Behavior; Particle Swarm Optimization; Capillary Pressure; Reservoir Fluids; Phase Equilibrium; Phase Envelope.

No description available.

No description available.

No description available.

The primary objective of this research is to develop a mathematical framework that could be used to model or predict the rate of fluid absorption and release in fibrous sheets made up of solid or porous fibers. In the first step, a two-scale two-phase modeling methodology is developed for studying fluid release from saturated/unsaturated thin fibrous media made up of solid fibers when brought in contact with a moving solid surface. Our macroscale model is based on the Richards’ equation for two-phase fluid transport in porous media. The required constitutive relationships, capillary pressure and relative permeability as functions of the medium’s saturation, are obtained through microscale modeling. Here, a mass convection boundary condition is considered to model the fluid transport at the boundary in contact with the target surface. The mass convection coefficient plays a significant role in determining the release rate of fluid. Moreover the release rate depends on the properties of the fluid, fibrous sheet, the target surface as well as the speed of the relative motion, and remains to be determined experimentally. Obtaining functional relationships for relative permeability and capillary pressure is only possible through experimentation or expensive microscale simulations, and needs to be repeated for different media having different fiber diameters, thicknesses, or porosities. In this concern, we conducted series of 3-D microscale simulations in order to investigate the effect of the aforementioned parameters on the relative permeability and capillary pressure of fibrous porous sheets. The results of our parameter study are utilized to develop general expressions for kr(S) and Pc(S). Furthermore, these general expressions can be easily included in macroscale fluid transport equations to predict the rate of fluid release from partially saturated fibrous sheets in a time and cost-effective manner. Moreover, the ability of the model has been extended to simulate the radial spreading of liquids in thin fibrous sheets. By simulating different fibrous sheets with identical parameters but different in-plane fiber orientations has revealed that the rate of fluid spread increases with increasing the in-plane alignment of the fibers. Additionally, we have developed a semi-analytical modeling approach that can be used to predict the fluid absorption and release characteristics of multi-layered composite fabric made up of porous (swelling) and soild (non-swelling) fibrous sheets. The sheets capillary pressure and relative permeability are obtained via a combination of numerical simulations and experiment. In particular, the capillary pressure for swelling media is obtained via height rise experiments. The relative permeability expressions are obtained from the analytical expressions previously developed with the 3-D microscale simulations, which are also in agreement with experimental correlations from the literature. To extend the ability of the model, we have developed a diffusion-controlled boundary treatment to simulate fluid release from partially-saturated fabrics onto surfaces with different hydrophilicy. Using a custom made test rig, experimental data is obtained for the release of liquid from partially saturated PET and Rayon nonwoven sheets at different speeds, and on two different surfaces. It is demonstrated that the new semi-empirical model redeveloped in this work can predict the rate of fluid release from wet nonwoven sheets as a function of time.

A superhydrophobic coating can be produced using a hydrophobic material textured with surface roughness on the micro-/nano-scale. Such a coating on the outside of a submersible body may result in reduced skin-friction drag due to a trapped layer of air in the coating. However, this layer may become unstable when subjected to elevated hydrostatic pressures, and a coating’s performance is compromised beyond a certain threshold (critical pressure). This thesis presents a numerical model for predicting the pressure tolerances of superhydrophobic coatings comprised of randomly deposited hydrophobic particles or fibers. We have also derived a set of force-balance-based analytical equations for predicting critical pressure in surfaces with ordered roughness, and compared our numerical model against it, observing reasonable agreement. The numerical model was then applied in a large parameter study, predicting critical pressure for coatings with a given set of microstructure properties.

Non-Aqueous Phase Liquids (NAPLs) include commonly occurring organic contaminants such as gasoline, diesel fuel and chlorinated solvents. When released to subsurface environments their spreading is a complex process of multi-component, multi-phase flow. This work has strived to develop new models and methods to describe the spreading of NAPLs in heterogeneous geological media. For two-phase systems, infiltration and immobilisation of NAPL in stochastically heterogeneous, water-saturated media were investigated. First, a methodology to continuously measure NAPL saturations in space and time in a two-dimensional experiment setup, using multiple-energy x-ray-attenuation techniques, was developed. Second, a set of experiments on NAPL infiltration in carefully designed structures of well-known stochastic heterogeneity were conducted. Three detailed data-sets were generated and the importance of heterogeneity for both flow and the immobilised NAPL architecture was demonstrated. Third, the laboratory experiments were modelled with a continuum- and Darcy’s-law-based multi-phase flow model. Different models for the capillary pressure (Pc) – fluid saturation (S) – relative permeability (kr) constitutive relations were compared and tested against experimental observations. A method to account for NAPL immobility in dead-end pore-spaces during drainage was introduced and the importance of accounting for hysteresis and NAPL entrapment in the constitutive relations was demonstrated. NAPL migration in three-phase, water-NAPL-air systems was also studied. Different constitutive relations used in modelling of three-phase flow were analysed and compared to existing laboratory data. To improve model performance, a new formulation for the saturation dependence of tortuosity was introduced and different scaling options for the Pc-S relations were investigated. Finally, a method to model the spreading of multi-constituent contaminants using a single-component multi-phase model was developed. With the method, the migration behaviour of individual constituents in a mixture, e.g. benzene in gasoline, could be studied, which was demonstrated in a modelling study of a gasoline spill in connection with a transport accident. / Flerfasföroreningar innefattar vanligt förekommande organiska vätskor som bensin, dieselolja och klorerade lösningsmedel. Spridningen av dessa föroreningar i mark är komplicerad och styrs av det samtidiga flödet av organisk vätska, vatten och markluft samt utbytet av komponenter (föroreningar) mellan de olika faserna. Detta arbete syftade till att utveckla nya metoder och modeller för att studera spridningen av flerfasföroreningar i mark: (i) En metodik utvecklades för att i laboratorium noggrant och kontinuerligt mäta hur en organisk vätska är rumsligt fördelad i en tvådimensionell experimentuppställning. Metoden baserades på röntgenutsläckning för olika energinivåer. (ii) Infiltration av organisk vätska i vattenmättade medier studerades för olika konfigurationer av geologisk heterogenitet. I experimentuppställningen packades olika sandmaterial noggrant för att konstruera en välkänd, stokastiskt heterogen struktur. Spridningsprocessen dokumenterades i tre detaljerade mätserier och heterogenitetens påverkan på flöde och kvarhållning av den organiska vätskan påvisades. (iii) Experimenten simulerades med en numerisk modell. Olika modeller prövades för att beskriva de grundläggande relationerna mellan kapillärtryck (Pc) vätskehalt (S) och relativ permeabilitet (kr) för detta tvåfassystem av vatten och organisk vätska. En relation infördes för att beskriva partiell orörlighet hos den organiska vätskan i porer vars halsar tillfälligt blockeras av vatten då mediet avvattnas. Vikten av att i de grundläggande relationerna ta hänsyn till hysteresis och kvarhållning av organisk fas visades. (iv) Olika Pc-S-kr relationer för trefassystem av vatten, organisk vätska och markluft testades mot befintliga experimentella data. En ny relation för hur slingrigheten (eng. tortuosity) beror av vätskehalten infördes i kr-S relationen och olika möjligheter för att skala Pc-S relationen analyserades. (v) En modelleringsmetodik utvecklades för att studera spridningen av flerkomponentsföroreningar. Med metoden kunde spridningsbeteendet hos enskilda, särskilt skadliga komponenter som t.ex. bensen särskiljas då ett bensinutsläpp i samband med en transportolycka simulerades.

Producing since 1995, Bozova Field is a mature oil field to be re-evaluated. When evaluating an oil field, the common approach followed in a reservoir simulation study is: Generating a geological model that is expected to represent the reservoir / building simulation models by using the most representative dynamic data / and doing sensitivity analysis around a best case in order to get a history-matched simulation model. Each step deals with a great variety of uncertainty and changing one parameter at a time does not comprise the entire uncertainty space. Not only knowing the impact of uncertainty related to each individual parameter but also their combined effects can help better understanding of the reservoir and better reservoir management. In this study, uncertainties associated only to fluid properties, rock physics functions and water oil contact (WOC) depth are examined thoroughly. Since sensitivity analysis around a best case will cover only a part of uncertainty, a full factorial experimental design technique is used. Without pursuing the goal of a history matched case, simulation runs are conducted for all possible combinations of: 19 sets of capillary pressure/relative permeability (Pc/krel) curves taken from special core analysis (SCAL) data / 2 sets of pressure, volume, temperature (PVT) analysis data / and 3 sets of WOC depths. As a result, historical production and pressure profiles from 114 (2 x 3 x 19) cases are presented for screening the impact of uncertainty related to aforementioned parameters in the history matching of Bozova field. The reservoir simulation models that give the best match with the history data are determined by the calculation of an objective function / and they are ranked according to their goodness of fit. It is found that the uncertainty of Pc/krel curves has the highest impact on the history match values / uncertainty of WOC depth comes next and the least effect arises from the uncertainty of PVT data. This study constitutes a solid basis for further studies which is to be done on the selection of the best matched models for history matching purposes.

Numerical simulations indicate that hydrocarbon migration and entrapment in stacked fault-bounded reservoirs are mainly affected by the following factors: charge time, faults, pressure and geological structures. The charge time for commercial hydrocarbon accumulation is much longer in oil-water systems than in oil-gas-water systems. Faults are classified into charging faults and 'back doors' faults other than charging faults in stacked fault-bounded reservoirs. The lower the displacement pressure of a fault, the higher its updip oil transportation ability. The downdip oil transportation ability of a fault is usually low and cannot cause commercial downdip oil accumulation. Back doors affect both hydrocarbon percent charge and hydrocarbon migration pathways. Updip back doors improve updip oil charge. The lower the displacement pressure of an updip back door, the more efficient the updip oil charge before 3,000 years. Back doors whose displacement pressure is equal to or higher than 28.76 psi are effective in sealing faults in oil-water systems. On the contrary, only sealing faults result in commercial gas accumulations in stacked fault-compartmentalized reservoirs. Otherwise gas is found over oil. Downdip back doors generally have few effects on downdip hydrocarbon charge. Geopressure enhances the updip oil transportation of a fault and improves the positive effects of updip back doors during updip oil charge. Geopressure and updip back doors result in more efficient updip oil charge. A physical barrier is not necessarily a barrier to oil migration with the aid of geopressure and updip back doors. The chance for hydrocarbon charge into reservoirs along growth faults is not equal. Any one of the above controlling factors can change the patterns of hydrocarbon charge and distribution in such complex geological structures. Generally, lower reservoirs and updip reservoirs are favored. Reservoirs along low-permeability charging faults may be bypassed. Gas can only charge the updip reservoirs. Both updip and downdip back doors can facilitate oil penetrating a barrier fault to charge reservoirs offset by the barrier fault. Interreservoir migration among stacked fault-compartmentalized reservoirs is an important mechanism for hydrocarbon accumulation and trap identification. The interreservoir migration is a very slow process, even though the displacement pressures of bounding faults may be very low.

No description available.